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Summary 


The biodynamiecs of poultry health is a function of a complex set of 
events and interactions between the host (poultry), disease causing agents 
and environmental factors. In today's poultry industry, although the 
presence or level of many infectious agents have been reduced and birds are 
kept in relatively controlled environments, spatial and temporal 
differences in managerial and environmental factors affect the prevalance 
rates of diseases amongst farms. This also leads to subsequent differences 
in condemnation rates at the time of processing. Since the birds! genetic 
potential and environmental determinants are closely monitored to attain 
maximai growth with minimal losses due to diseases, the condition of birds 
in a flock tends to be homogeneous. It may then be possible to determine 
the health status of such a population of birds by examining the 
performance records and other possible indicators of health and by 
inspecting a representative sample from the flock before slaughtering the 
entire fiock. Amongst the important health status indicators are such 
factors as strain of birds, hatchery management, types of brooding ana 
housing, sanitation, vaccination, medications, feed utilization,morbidity 
and mortality, macro and micro climate and carcass inspection procedures, 
Each of these affect condemnation rates due to diseases either directly or 
indirectly. If the variables which are causally related to condemnation 
rate could be identified and their significance quantified, then one could 
device appropriate mechanisms to monitor the condemnation rate due to 
diseases and use it as an adjunct: to existing carcass inspection 
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In recent years, inflation and the poultry industry's tremendous 
growth has resulted in a substantial rise in the cost of inspection. 
Because of these costs, a closer scrutiny of existing inspection methods 
has become necessary so as to determine if more efficient and less costly 
Ways could be devised. The purpose of this study, therefore, was to 
explore the possibility of developing a predictive model for the prevalence 
of disease conditions in broiler chickens auring routine post-mortem 
inspection. Based on such a model, it was hoped that less costly 
inspection procedures, which would not compromise public health safety, 
could be developed. With this in mind and with the goal being the 
isolation of those variables which will be useful predictors of 
condemnation rates, the need to systematically examine the 
interrelationships and interactions of all potentially useful predictor 
variables was necessary. For such a task, systems analysis and 
multivariate analytic methods were selected. The systems concept is a well 
tested problem-solving tool widely accepted in engineering, economics, 
ecology and human medicine. In veterinary medicine, it is relatively new 
but is now being applied to solving epidemiologic problems. 

The application of multivariate techniques to develop predictive or 
discriminative models is also widely used in various scientific and social 
desciplines where the background in the use of quantitative methods is well 
established. It is a standard statistical tool, the basis of which rests 
on sound mathematics. Aithough not widely used in veterinary medicine, 
examples of the applications of such analytic methods are available. 

In the area of SBULEEY health, the application of the systems concept 
and multivariate analytic models appears limited. However, in one recent 


Study, the relationship between various turkey porduction factors, and 


condemnation and downgrading in turkeys were examined by using factor 
analysis. Using the 1977 live bird production and processing data from one 
processing plant in California, the study indicated that, week of year 
placed, climiate on day of processing and feed conversion may have 
predictive value. 

Since the application of systems analysis and multivariate models in 
veterinary medicine and more specifically in poultry health research is 
rather scanty, it is believed that the results of this study may have 
broader benefits in examining comparable scientific problems. 

The poultry systems approach: The rationale of the study relied on 
the fact that since it is highly unlikely that diseases suddenly arise in 
transit (although existing infections could be exascerbated as a result of 
transit stress) or while birds are under processing, it was rational to 
think that disease conditions that lead to condemnations during inspection 
must nave had some time to develop and be recognized as lesions. In that 
case, the infection must have taken place while the flock was still on the 
farm, This, necessitated a critical examination of the farm ecosystem to 
isolate and identify variables which would be useful in predicting 
condemnation rates; this was the task of the research design. Therefore, 
environmental, managemental, as well as host and agent factors which may 
affect the incidence of poultry diseases and the condemnation rate at 
Slaugher plants, were scrutinized via causal diagrams and multivariate 
epidemiologic models designed for possible predictive and/or discriminative 
purposes. Since poultry (broiler) production involves a sequentially 
predetermined set of factors which includes breeder flock sources, 
hatching, brooding, growing and processing, the identifications of 


potentially useful predictors of the condemnation rate at slaughter were 


iv 
designed to take the total poultry production ecosystem into consideration. 

Procedure for identification of variables: Three subsystems were 
considered; the hatchery,the broiler house and transit and processing. In 
the natchnery subsystem, the task was to identify the variables which had 
potential causal influence on condemnation rates. 

The infection in the hatchery may arise via genetic or egg-borne 
transmission (from breeder flock), or as a result of exposure to infectious 
agents in the hatchery. The subsequent stages which follow hatching are 
brooding and growing, which take place within the broiler house subsystem. 
Again, the birds are exposed to various risk factors in such a setting and 
epidemiologicaliy relevant factors were identified with the aid of the 
systems diagram. Tne last stage in the sequence of poultry production 
phases were factors related to transit or those that arose while birds were 
being processed, referred to as the transit and processing subsystem. 

Sampiing: Five cooperating poultry firms with 5-7 growers per firm 
were selected to form the study sample. To ensure a representative cross 
section of farms, the study units were selected from three strata of good, 
average and poor production performance based on criteria provided by the 
Parent firm. The sampling time frames were longitudinal during which time 
the health status of the population of cohorts and specified multivariable 
data which affects the condemnation rate were collected. By definition, a 
cohort of birds that were sequentially monitored during the three phases of 
broiler operation (hatchery, broiler farm, and processing plant), were 
referred to as. a brood cycle. It requires about 50 days to raise a bird.to 
the gesirable marketing weight as a broiler. The data collection time 
frame was one complete year during which time data on 5 brood cycles were 


coliected. After an initial data gathering visit to a hatchery where the 


























study farm chicks were hatched, each study farm was visited at two 
aifferent time periods; one when a given brood of birds were about four 
weeks old, and again when they were about six to seven weeks old. The last 
Stage of the sample data collection time frame involved obtaining 
condemnation records of study flocks the gathering of serum for serological 
analysis and whole bird samples for gross pathological studies in the 
laboratory. 

Types of data: Preexisting data (e.g. condemnation rate at 
processing, weekly mortality rate etc.) were collected from records kept by 
the firm, the farm or from the processing plants. Others were generated 
from field observations (e.g. management variables, weekly feed utilization 
Rave eu o)% 

Finally, laboratory data were generated from serological and necropsy 
profile studies where a predetermined smaplesize of 200 randomly selected 
birds were examined. Sera were tested for three indicator diseases using 
the HA/HI test. These dissease agents were Newcastle, disease virus 
‘Mycoplasma gallisepticum and Mycoplasma synoviae. The necropsy task 
involved gross pathological examination of the birds for any abnormal 
tissue changes which could be attributed to diseases. 

Data processing: Since the raw data collected from the various data 
generating centers had to be stored, retrieved and analyzed using a 
computer, a specific file handling software was developed for local use. 
Tne software used BASIC language and was designed in an interactive mode. 
The data were stored on floppy diskettes in a model II Radio Shack 
microcomputer. At the time of data analysis, the information was 
transferred to a larger central computing facility. The analysis of the 


data was performed using BMDP statistical software, on a VAX 11/750 


vi 
computer at Tuskegee Institute. The specific RMDP statistical programs 
used were: P2D for descriptive statistics, P6) for examining scatter plots 
of variables, P1R and P2R foi: multiple and stepwise regression analysis and 
P74 for conducting stepwise discriminant analysis to develop a 
classification function into high or low condemnation groups. The final 
screening of variables to arrive at the best set of predictor or 
discriminator variables was accomplished by a computerized stepwise 
Selection process which utilized forward and backward stepping techniques. 
In this procedure, the criterion for adding or deleting a predictor 
Variablie equivalently involved reduction in sum of squares, partial 
correlation coefficient, or an evaluation of the F statistic. Once a set 
of predictor variables were identified, the resulting model was examined 
further by performing residual analysis to evaluate the adequacy of the 
modei. A variety of residual plots were employed to identify the nature of 
fit, outliers, deviation from normality and other deficiencies. 

Predictor model: Tne data base contained a maximum of 154 cases for 
141 variables. These variables formed the initial data matrix upon which 
the various PUDP statistical programs were run. Since the objective was to 
develop a predictive model, only the multivariate data analyses are 
presented. The first step was to reduce the lengthy list of variables via 
systematic screening to a more manageable level, which was arbitrarily set 
at about 30-50 predictor variables. The final model for a subset of 
predictor variables was then selected from these. In such a variable 
reduction step, one consideration was to delete all variables which were of 
local significance to tne study region so that only variables which were of 
locai significance to the study region could be generalized to other areas 


were utilized in the analysis. Included in such omissions were the 


Significant variables representing hatchery source, vaccine source and 


location of air intakes or exhausts. 











Via regression analysis, the best predictor model was determined to 


be: 


19. 3543-0. 21023X +0, 10855X,),+0.20511% 


Wnere: 


63 109 

#0. 3227K 159-0. 4371% 596-0 .00702K,,, 40, 19356 Yor 

(Re= 0291" Ro = O.82;°F ratio = 19:07) (1) 

Y5 = predicted overall condemnation rate (%) due to 

diseases, 

X63 = the average hatchability rate (%) for the hatchery, 

Xoy = the total number of genetic strains of birds on the 
farm, 

X109 = the disease rate variable (%) derived from 
necropsy examination of poultry 

X409 = water spacing per bird in inches 

X505 =) total precipitation (inches) during a brood cycle, 

X503 = the distance from the poultry farm to the processing 
plant (miles) 

Yor = the overall condemnation rate due to disease (%) 


during one previous brood cycle. 


This equation was obtained at the end of the forward stepping 


regression analysis program. The finding was consistent with the 


underlying biological expectation in regards to the variables with a 


positive relationship to Ys. For example, if the previous 


condemnation rate (Y57) was high, one may expect the current 


condemnation rate or the response variable (Y,) to also be high. The 


positive relationship between Xan and Y5 says that if the number of 


genetic strains kept on a poultry farm are large, this leads to high 
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condemnation rates. This of course couid be a reflection of 
limitations in the management or that the chances of getting poor 
goers from among a larger genetic pool is greater than that from a 
few which would have been selected specifically for good performance. 
A similar management related argument could be made for X409° The 
variables with an inverse relationship, Xp3e Xo05 and % 503 are not so 
Straightforward to explain except for the first case. If the 
hatchability rate (X53) is high, this could be indicative of good 
Management and therefore the returns in lower condemnation rates at 
processing may be expected. The model indicate that the total 
precipitation during the brood cycle (X5 95? is inversely related to 
condemnation rate. Although this could be questionable, one reason 
here may be that significant precipitation is observed during the 
spring or early summer months which may be the most conducive tines 
for raising poultry. This may be partly due to the cleansing effect 
that rain may have on infectious agents in the air. The variable for 
aistance from the farm (%503) is somewhat peculiar in the sense that 
one would in fact expect to see a positive influence of X 503 on Y, 
i.e. as the distance increases, the birds are exposed to more stress 
due to travel which may then exascerbate existing abnormalities 
leading to higher condemnation rates. However, the relationship in 
the model was inverse. This can only be speculated here in that 
possibly due to the fact that transit stress is known to be harmful 
to animals in general, the management may nave provided above average 
care at the time of transporting of birds to the processing plants. 
An examination of the standardized regression coefficients given 


in equation 2 is useful in providing comparative information on the 


viii 


importance of each predictor variable. 
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In this case, X 109 Was the most valuable predictor followed by 


Xx and finally X,, and X For the predictor 


223° lor 1 129° 


2 
equation, the R and R- values were 0.91 and 0.82 respectively. This 


X505* *63: 


was a very high correlation coefficient and the equation successfully 


explained 82% of tne variability in Y It also had a significant F 


aye 


ratio of 19.07 which showed strong lineraity between Ys and the 
independent variables. 

An essential part of any regression analysis is a careful 
examination of residuals to ensure that the assumptions of least 
Squares theory are fulfilled. Tnerefore, residuals of the predictor 
model were scrutinized with the aid of computer plots. In the first 
plot, the objective was to examine the distribution of ri about the 

. +e A. 2 ° 
mean (zero). A plot of yi vs the predictea Yi indicated that the 
residuals appeared to be randomly cistributed around the zero mean in 
a balanced manner i.e. about half of the cases were positive (above 
0), while the rest were negative (below 0). The residuals were all 
within an acceptable standard deviation of + 2.0. In a second case, 
om Re 
the squared value of the residuals (ri) were plotted vs the 
: _ f ss 

predicted Yi. As expected, this snowed a ciustering of most of the 

aC : : 
Xi close te zero. Therefore, there were no major problems with 
outliers and in terms of model inadequacy. 
Discriminant models: To perform the analysis for developing the 
linear discriminant functions to classify a flock into high or low 


condemnation groups at tne time of processing, various 


decision-making options were tried. These included criteria for high 
or low condemnations set at 5%, 2%, 1%, 0.5% and 0.4%. Inese 
demarcations were set somewhat arbitrarily but with the idea of 
providing 2 variety of options for consideration. 

Tne dependent variable used to develop the discriminant models 
Was the overali condemnation rate due to diseases (Y,)- To perform 
the stepwise discriminant analysis, the F-to-enter and F-to-remove 
values were set at 4.0 and 3.996 respectively. For the two 
categories of condemnation rates (high or low), prior probabilities 


were set at 0.5. 


Two cases were selected and presented here for demonstration 


purposes. 
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Case i: Demarecation between high and low condemnation rate set at 1% 


LGE= 4091 .15234+114.77253%¢ 447 HMO3K, +4 89433X 


63 
HY IS 9BBX 6 4-18 .59022% 4411. 9N 752K pg thS. 3914X 50, 
4. 75752Y 574346 2482X, (3) 
HC = -3969.5529+107 671%, 0444. 3BHU5X, +4, 1A753K, , 
+H, 3899% 6-17. 56498X 412. 64814 9445,0957 1% 50. 
-0. 93103Y,7#31.05351X, (4) 


Percent correct classification overall = 91.9 
Percent correct classification for under 1™ = 935.4 
Percent correct ciassification for over 1% = 87.5 


Case 2: Demarcation set at 0.5% 


LC = —1863.71973+78.5039 1X¢q+77 .2542K 6 +28 6B584X, 

~85.92729X, ,9-8.73708X 59041, 32752%,,, (5) 
HC = ~4699 «33487476 B25NOK. +75. BNZ22K¢ 429 30284X, 

~83.62459K 1 59-7. 796 38k a 9541. 21017 X55 (6) 
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Fercent correct classification oversil = 9" .5 
Percent correct classification for under 0.5% = 88.2 
Percent correct classification for over 0.5% = 109 
Where: 

LC = equation for low condemnation rate group, 


HC = equation for high condemnation rate group, 


Xg = strain group 

Xan = total number of strains on farm 
Xe = incubator type 

X64 = incubator type 


Xoo = average fertility rate (%) for hatchery 


X63 = average natchability rate (%) for hatchery 

Xen = average fertility rate (%) for brood of birds 

X109 = average disease rate (%) based on necropsy at 7 weeks 
Xi09 = water spacing per bird in inches 


X205 = average total precipitation for study brood (inches) 


os 
[ 


= distance from farm to plant (miles) 
y, = previous condemnation rate (%) due to diseases for study 
farm 

In each cf the discriminant models, a classification function to 
specify whether a flock should be categorized into a high (HC) or low 
condemnation (LC) group are provided. In each case, a flock was assigned 
to the group (HC or LC) with the largest value of the classification 
function. Tnese equations were coupled with the respective correct 
classification probabilities. The probability of correctly classifying a 
flock into nigh condemnation category ranged from 37.5% to 100%. The range 


for classifying into low condemnation category was 88.2% to 93.1%, with an 
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Overail correct classification probability of 91.9 to 94.6%. This 
jnaicates tnat the accuracy of the model in providing decision-making 
information was, therefore, reliabie, 

It should be noted that the objective in this study was to develop a 
predictor model wnich could serve as an adjunct to the existing poultry 
inspection procedures and standards. Therefore, the results presented are 
visualized in terms of whether a decision pertaining to the current 
inspection procedure is to be made may involve either a scaled down 
procedure which in effect means that the processing line speed may increase 


Or to maintain the status quo currently operational in the various 


processing plants. 

In order to effectively utilize the models, two prerequisites should 
be fulfilled. Firstly, information for the variables identified in the 
equations should be made available by the firm/farm at the time of 
processing or immediately prior to processing of the birds. Secondly, the 
veterinarian, must take the data for the variables selected and generate a 
vaiue for Y, wnich is the predicted condemnation rate for the particular 
fiock to be processed. Finally, using the values computed, the 
veterinarian decides on whether to introduce scaled down inspection 
procedure or maintain the regular inspection routine. Based on current 
information, the veterinarian could for example use a wicrocomputer and 
obtain a set of instructions to facilitate the decisionmaking task for the 
type of inspection best suited to individual flocks. Such a computer 
program was prepared using BASIC. The program is interactive with a menu 
format where one persues available alternatives and responds to those 


questions until the objective has been achieved. 
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Note that, via such microcomputers located in the various poultry 
processing plants throughout the country, one could maintain a network of 
information gathering devices with connections to a central data base in 
Some specifiea place. Since such a data base would be inherently rich and 
reliable in the case of poultry, one maybe in a position, not only to use 
microcomputers for implementing the modeis in the field, but one could also 
utilize the aata base for other problem-solving tasks. Therefore, it is 
feit that one of the best alternatives for implementing the use of the 
predictor modaeis in the field could be via the Peed ea cn of 
microcomputer to facilitate data handling and decision-making in poultry 


inspection. This would be of great benefit to the poultry industry and to 


the veterinary profession itself, whereby, for example, some of the 
expertise of veterinarians and poultry inspectors may be diverted to expand 


the utilization of current information technology in their respective areas. 
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INTRODUCTION 


The precursor to a problem-solving task such as in the case 


at hand, involves firstly a definition and a detailed description 


of the problem under study so that the objectives to be attained 


could be better focussed. In this case, a brief review of poultry 


inspection and the problems associated with it as well as a 


description of the poultry industry, followed by a review of the 


biodynamics of poultry health provides an adequate background to 


comprehend and appreciate the complexities involved in this study. 


Uae 


The Problem: The poultry industry has seen tremendous 
growth since 1959, when inspection of poultry first became 
mandatory. The total volume of poultry products slaughtered 
under Federal inspection rose from 5 billion pounds (1bs) in 
1960 to 16.5 billion ibs in 1981. More specifically, there 
has been an increase in broiler production (1.5 billion 
broilers slaughtered in 1960 to 4.1 billion in 1981, (1-5). 
Such increases in production have caused higner inspection 
costs on one side, while also leading to higher speed 
processing lines (23 birds per minute, per inspector on the 
average) to keep up with higher processing rates. This, no 
doubt, has resulted in increased efficiency of processing 
plants, but it may have diminished the accuracy and 
consistency of inspection. Correspondingly, the cost of 
inspection has steadily risen and when compounded with the 
inflation rate, a closer scrutiny of existing inspection 


methods and other research into possible cost reducing 


inspection alternatives seemed timely. 

The responsibility of poultry inspection falls under 
the jurisdiction of the Food Safety Inspection Service 
(FSIS) of the United States Department of Agriculture 
(USDA). Annually, over 4 billion slaughtered birds to be 
inspected by FSIS personnel in accordance with the amended 
Poultry Product Inspection Act (PPIA) of 1968. This program 
is responsible for assuring poultry and poultry products 
slated for commerce in the United States (U.S.) are 
inspected to provide confidence to consumers that such 
products are wholesome, unadulterated and truthfully 


labelled. The inspection law covers raw carcass poultry 
(except for small farm flocks), as well as all processed 
poultry products such as frozen dinners and soups. About 
3000 processing plants which slaughter and/or further 
process poultry are under inspection by FSIS. Only 245 of 
these slaughter young chickens. Chickens comprise 96 
percent (%) of the slaughter and turkeys, ducks, geese and 
guinea fowl comprise the remaining 4 percent (16), 

The process of inspection, although standard and 
comparable from plant to plant, requires careful and 
diligent supervision and execution of tasks by well trained 
personnel. It consists of ante-mortem (before death) and 
post-mortem examinations. Ante-mortem inspection for 
disease involves making spot checks on the health of birds 
as they arrive at the processing plant so that those which 


are already dead, or those which appear unfit for slaughter, 


are automatically condemned and destroyed. The slaughtered 
bird carcass and its internal organs are individually 
examined on high speed processing lines for diseases or 
conditions which would cause all or part of the carcass 
unfit for human consumption. To assure uniformity and 
effectiveness of the inspection process, USDA veterinarians 
Supervise the post-mortem inspection procedures and other 
duties performed by food inspectors. 

In recent years, inflation and increased productivity 
in the poultry industry have been straining the inspection 
costs, which is borne by appropriated Federal funds except 
for inspectors' overtime work. To help control costs and 
ensure consumer protection, FSIS has been developing and 
reviewing inspection procedures to determine more efficient 
Ways to inspect poultry and poultry products. For example, 
most young chickens or fryers are now inspected under the 
procedure known as a modified traditional inspection which 
makes use of a mirror to reflect the back side of a bird, 
thus reducing the need to manually turn over the carcass as 
in the traditional method. This newly implemented procedure 
(April 1979) has been shown to save inspection time and the 
inspectors' task has also become less tiring (6). 

Other ways of lowering inspection costs are also being 
considered, and the one that is the central theme of this 
report is that of flock testing. This new approach to 
inspection may offer considerable savings in manpower. This 


potential could be realized if it is possible to predict the 


ies 


types and prevalence of disease conditions in poultry flocks 
before they are slaughtered; and then develop a scaled down 
inspection procedure which could be safely used to inspect 
fiocks where the predicted incidence of disease is low. The 
rationale for such an approach is quite simple. Most 
poultry, especially broilers, are raised and brought to 
market under highly controlled conditions by an industry 
which is vertically integrated such that one firm has 
control of breeding, hatching, feeding, growing and 
Slaughtering. Since the birds' genetic potential and 
environmental determinants are closely monitored to attain 
maximal growth with minimal losses due to diseases, the 
condition of birds in a flock tends to be quite uniform and 
homogeneous. It may then be possible to determine the 
health status of such a population of birds by examining the 
performance records and other possible indicators of health, 
and by inspecting a representative sample from the flock 
before slaughtering the entire flock. In light of the 
decline in condemnation rates due to diseases among poultry, 
flock testing which relies on population health indicators 
may prove to be useful and appropriate. 

The Objective: It is with this background that the contract 
to conduct a "Flock Testing Feasibility Study" was awarded 
to the School of Veterinary Medicine of Tuskegee Institute 
on September 30, 1980. The specifie objective of this 
contract was to determine whether it is possible to predict, 


accurately and by practical means, the presence and 


ancidence of those conditions in broiler chickens observed 
during routine post-mortem inspection and which may be a 
cause for condemnation of all or part of the carcass. 

The purpose of the research design was directed at 
exploring the possibility of developing a predictive model 
to determine prevalence of disease conditions in broiler 
chickens during routine post-mortem inspection. Based on 
such a model, it is hoped that less costly inspection 
procedures, which would not compromise public health safety 
could be developed. 

The contractual agreement for this project was signed 
on September 30, 1980, and the date of approval by the 
Office of Grants Management of Tuskegee Institute to proceed 
with the project was made on November 1, 1980. The actual 
implementation of the project officially began on February 
15, 1981, with its termination date slated for 
June 30, 1982. Initially, activities and estimated times 
for completion of the tasks were outlined in a graphical 
procedure for analytical and synthetical evaluation and 
program review (GASP) chart. A modified version of that 
chart is indicated in figure 1. Since the study was 
centered around the poultry industry in Alabama, a brief 


review of the study area appears pertinent at this time. 
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Table 1 


KEY FOR GASP CHART 


TASK - FLOCK TESTING FEASIBILITY STUDY 








REFERENCE ACTIVITIES DESCRIPTION ESTIMATED COMPLETION 
DATE 
A Proposal preparation September, 1980 
B Finalize contract November, 1980 
C ~Develop an Epidemiologic format 
for flock testing study May, 1981 


-Initiate contact of farms and 
processing plants 


-Prepare data collection and coding 


format 
D Finalize arrangements with farms 

and plants March, 1981 
E -Hire personnel April, 1981 


-Purchase equipment/supplies 


F Conduct pilot studies March, 1981 
G Data collection March, 1982 
H -Develop file handling software February, 1982 


-Code data and enter on computer 
A! Data analysis and interpretation June, 1982 
J Prepare final report June, 1982 


K Begin Phase II of study Undetermined 


ee 


ae THE POULTRY INDUSTRY IN ALABAMA 
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Background of Industry: In broiler poultry production, 
Alabama is ranked third in the U.S. preceded by Arkansas and 
Georgia (Table 2). The poultry industry in the state is 
sizeable and diverse in nature, ranging from small 
independently operated farms to large integrated firms with 
10,000 chickens or more per grower. There are 11 major 
integrated broiler companies in Alabama (7). The number of 
growers per integrated firm is variable as is the size of 
the poultry population per grower. Most of these growers 
are iocated in the northern and southeastern parts of the 
State. (Figure 2). 

In 1980 and 1981, 39,309,000 and 42,149,000 birds 
respectively were inspected by FSIS personnels in Alabama. 
During this period, the condemnation rate in the state 
varied from 1.32% to 2.19% on a monthly basis; the national 
figure for condemnation varied from 1.57% to 2.25% for the 
Same period. The average annual condemnation rate for 
Alabama in 1979 and 1980 was 1.93% and 1.72% respectively. 
The figure for the same period for the U.S. was 2.01% and 
1.84% respectively, indicating that the rate was declining 
(Table 2) in the state as well in the country as a whole. 
Tne disease specific condemnation rates for the five major 
poultry producing states for the period 1976-1982 is given 
in Table 3. Tne values for 1982 are incomplete. 

In Alabama, as well as in other states, the vertically 


integrated firms operate their own breeder flocks, hatching, 
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Table 2. 


COMPARISON OF YOUNG CHICKENS SLAUGHTERED UNDER FEDERAL INSPECTION IN 
THE TOP FIVE (5) BROILER PRODUCING STATES (1975-1982) 


(Number of Chickens In 1000's) 
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No. Inspected 49,857 43,711 STL L 
Post. mACOnd et.) 1266 eos basse) 










21,060 
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No. Inspected 36,838 34 20 Bi 935 22), DOO He Ar) 
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Fig.2 MAJOR BROILER PRODUCING COUNTIES 
IN ALABAMA 
Flock Testing Feasibility Study 
Tuskegee Institute (1981/82) 
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growing, feeding and slaughtering. These operations are 
closely supervised. However, there are managemental 
variabilities amongst firms. Each firm has growers which 
could be categorized as above average, average and below 
average based on performance characteristics which are 
further defined under methodology for stratified sampling. 
The major poultry firms in the state are either national 
or independents companies. Such a combination of 
independent and national firms provide the state with a 
representative blend of the poultry industry in the U.S. 
Typical Poultry Production Process: A typical broiler 
raising enterprise consists of a source of breeder flocks 
which provide the high quality eggs for hatching. The 
genetic make up of breeder flocks in Alabama are mainly of 
the Corbett, Arbor Acres, and Peterson strains. The eggs 
from such breeder flocks are taken to hatcheries where under 
optimally controlled temperature and relative humidity 
conditions, the eggs are set and hatched after a period of 
21 days. An important factor in the management of the 
hatcheries is the control of incubator temperature and 
sanitation of the environment in the hatchery. Hatching egg 
quality, fertility, hatchability and chick quality are an 
interrelated set of events that ultimately affect the cost 
of production. So, the sanitary management of the hatchery 
becomes crucial. 

Sanitation in a hatchery typically involves washing, 


fumigating and disinfecting. The disinfectants utilized are 
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quaternary ammonium compounds, organic iodine combinations 
and phenol. In the application of these chemical agents, 
manufacturer's recommendations are followed, 

The cleaning and disinfecting procedures vary somewhat 

depending on the parts of hatchery as described below: 

- Setters- in some hatcheries, the setters are sprayed 
three times per week with a solution of quaternary 
ammonia. In addition, broken eggs and debris are 
removed from the bottom portion of the setter on a 
weekly basis. Every four months, the setters are 
scrubbed with a commerical detergent and rinsed with 
hot water using a high pressure sprayer, followed by 
a disinfectant such as quaternary ammonia compotmd, 
In others, the setters are washed every six months 
with a commerical detergent and rinsed with cold 
water under high pressure. 

- Hatchers- after the chicks are removed, the hatcher 
is rinsed with water and drained, Scrubbing with a 
detergent solution and rinsing with hot water under 
high pressure ensues, This is followed by the 
application of a disinfectant solution such as a 
quaternary ammonia compound or iodine solution. 

- Vaccinating machines— are hand cleaned daily with a 
detergent solution. Rinsing with cold water follows. 
Needles are autoclaved every day and changed twice a 
day. 


~ Floors and hallways— are washed with hot water under 
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high pressure and sprayed with a disinfectant such as 
a quaternary ammonia compound every hatch day. In 
some hatcheries, the work rooms are fumigated three 
times per week with a quaternary ammonia solution. 

~ Chick buses- are washed with hot water under high 
pressure two times per week in some hatcheries and 
once a week in others. A solution containing 
quaternary ammonia is sprayed on the inside of the 
buses. 

~ Hatching crates and egg flats—- are washed with hot 
water and disinfectant on a tray washing machine at 
the end of every hatch. 

- Chick delivery crates— are washed with a disinfectant 
and hot water twice a week on a crate washing 
machine. 

The above brief glimpse of sanitation practices in a 
hatchery is one of many managerial responsibilities. After 
hatching, baby chicks are kept in the hatcher for about 24 
hours. During this time, the chicks are not fed. In many 
of the hatcheries studied, chicks were vaccinated for 
Newcastle, infectious bronchitis and Marek's disease before 
they were moved to the broiler houses. The baby chicks were 
then transported, via chick buses, to the broiler house. In 
other hatcheries, the chicks were vaccinated within 10-14 
days of hatching. Again, good sanitation practices and 
a warm and comfortable environment during the transition 


period is important. The chicks are either placed in 
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partial or whole house brooding. Partial house brooding is 
an attempt to save cost (energy/heat, space) and maintain 
the chicks in a smaller confined space, where more birds can 
huddle, to keep warm, under a few brooders. 

The responsibility of brooding and growing birds to 
marketable weight, at about 50 days of age, rests on each 
grower. The birds are monitored for feeding, water 
consumption, disease problems and variations in microclimate 
Within tne house. Growers usually cull a proportion of 
their undersized birds during the first few weeks, Through 
this practice, the cost of feed can be reduced since the 
undersized and inefficient birds are depopulated. During 
the growout period, good farm management practices coupled 
with a planned disease control program, proper feeding and 
Sanitation, are crucial factors that determine the health 
Status of the population. At the end of the 50 day growing 
period, a catching crew is brought in and the birds are 
loaded onto trucks and transported to the respective company 
Slaughter plants. The grower then has about 7-10 days 
during which time the litter may be changed and the houses 
cleaned and disinfected in preparation for the subsequent 


growout or brood cycle. 


Processing plants in Alabama: The processing plants in 


Alabama appear to be representative of the FSIS supervised 
plants in other parts of the country. There were 25 
processing plants operating in Alabama in 1981-82 (8). Some 


of the plants operate a two shift system, doubling up on 
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night and day processing. A majority of the plants operate 
on a regular eight hours a day basis. In terms of 
inspection, both the traditional and modified traditional 
procedures are utilized. In the former case, an individual 
inspector examines both the external aspects of a bird and 
its internal organs, to determine if it is fit for food, In 
the newer approach, a team of three inspectors are placed on 
each line with responsibilities split into an "outside" and 
"inside-viscera" inspection of each bird. Inspection of the 
external parts of the birds is facilitated with the aid of a 
mirror so that the outside "parts" inspector does not have 
to handle the bird; a mirror behind the carcass makes all 
surfaces visible. Two inspectors handle the "inside- 
viscera" or internal organs of the birds. The "inside- 
viscera" examination has been redesigned so that fewer hand 
motions are used to check the inside of the carcass and time 
can be saved, making the process less tiresome. Upon 
completion of inspection, the birds are chilled and packaged 
for distribution to consumer outlets. 


Is Alabama Representative of the Poultry Industry in the 


a ee 


South and in the United States? Among the eleven (11) major 


integrated broiler companies in Alabama, five are national 
companies that also have broiler operations in other states; 
the other six are independent companies owned by individuals 
or families in the state (Table 4). Among the five firms 
which cooperated in this study, three were national and the 


other two were independent companies. Their poultry flocks 


Table 4 
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Broiler Processing Plants and major companies operating in Alabama 


(1981-1982) 








Name of Firm Location 
City County 
Cagie's Incorporated Collinsville DeKalb 
Colonial Poultry Co. Albertville Marshall 
ConAgra Poultry Co. Athens Limestone 
Ww " " Decatur (office) 
He " rt Enterprise Coffee 
Gold Kist Poultry Co. Cullman (office) Cullman 
_ ay MW Boaz Marshall 
Ws " " Li Guntersville Marshall 
Ms " " uN Trussville Jefferson 
Golden-Rod Broilers Cullman Cullman 
" " n Greensboro Hale 
Marshall Durbin Co., Inc. Hayieville Jefferson 
ui " " " Jasper Walker 
Peco Foods, Inc. Gordo (office) Tuscaloosa 
i a “ Tuscaloosa Tuscaloosa 
Poultry Products Co., Inc. Montgomery Mongtomery 
Southland Broilers, Inc. Enterprise Coffee 
Spring Valley Foods, Inc. Oxford (office) Etowah 
" ut ML " Blountsville Blount 
i " u u Heflin Cleburne 
u , Ww Us Gadsden Etowah 
“3 i: " i Ashland Clay 
Wayne Poultry Union Springs Bullock 
" " Decatur Morgan 
ys Wy Albertville Marshall 


SOURCE: Food Safety and Quality Service, Poultry and Dairy 


— 


Quality Division, List of Plants Operating Under USDA 


Poultry and Egg Grading and Egg Products Inspection 
Program (Washington, D.C., United States Department of 


Agriculture. 31979).en.mo. 





























were located in the major broiler producing areas of the 
state. 

Based on the brief presentation given above, Alabama 
appears to represent the diversity and complexity of the 
poultry industry in the U.S. Tne state is in an area where 
more than 60% of the nation's broilers are produced 
(Arkansas, Georgia, Alabama, North Carolina and 
Mississippi). Thus, it was logical to conduct this study in 
Alabama where the poultry industry is well established and 
occupies a significant economic position in the hierarchy of 
agricultural industries both in the state as well as in the 


entire U.S. 
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Sie THE BIODYNAMICS OF POULTRY HEALTH - A SYSTEMS VIEW 


salle 


Poultry Population Health Dynamics: Poultry health is a 
function of a complex set of events and interactions between 
the host (poultry), disease causing agents and environmental 
factors. In today's poultry industry, although the presence 
or level of many of the infectious agents have been reduced 
and birds are kept in relatively controlled environments, 
differences in managerial expertise and diligence, as well 
as variabilities in environmental factors, still remain 
Significant. Such variabilities eventually lead to 
differences in prevalance of diseases amongst farms and 
subsequent differences in condemnation rates at the time of 
processing. In order to reduce the cost of production in 
terms of factors which have influence on poultry health, a 
detailed review of the literature and an understanding of 
biomedical processes as these relate to poultry health is 
necessary. 

As presented earlier, since one deals with relatively 
uniform flocks on a given farm and the major variability is 
from farm to farm, such a review of the literature 
emphasizes population health status indicators of individual 
flocks of broiler chickens. To be consistent and systematic 
in such a review, one should preferably start from the 
breeder flock source and proceed through hatching, brooding/ 
growing and finally, inspection of bird careasses. In light 
of this premise and to correlate the subsequent 


identification and specification of variables which were 
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used as health status indicators in this study, a sequenced 
review consisting of strain of birds, hatchery management, 
brooding, housing, sanitation, vaccination, medications, 
feed utilization, morbidity and mortality, climate and 
inspection procedures will. be examined; as these may affect 
ultimate condemnation rates due to diseases. Along with a 
literature review of the biodynamics of poultry health, and 
Since we utilize a systems approach and multivariate 
analytic models to analyze the data, a brief review of the 
literature on the application of such techniques in 
veterinary medicine may be informative and subsequently 


follows. 


The multiple determinants of poultry health: 


- Strains of birds- total performance and resistance or 
susceptibility to disease varies amongst strains or 
breeds of birds. Some are resistant to one disease 
but not to another. For example, several lines of 
chickens with genetic resistance to Marek's disease 
have been selected and maintained in the laboratory 
(9). Thus, Knowing which strain of birds is more 
susceptible to a specific disease that may result in 
high condemnation, could be helpful in formulating a 
method to predict condemnation. There are many 
varieties of poultry, but the more popular meat 
production male lines are: Hubbard, Cobb, Peterson 
and Indian River. The popular female lines are: 


Hubbard, Arbor Acres, Cobb, Indian River, Pilch and 


Shaver. Therefore, the genetic potential of a strain 
of birds, or the performance of the breeder flock as 
measured by such indices as genetic quality, 
fertility and hatchability of hatching eggs and chick 
quality are very critical in determining the health 


Status of a flock of birds. 


- Hatchery management- due to a less competent immune 


System newly hatched baby chicks are extremely 
susceptible to various diseases. Thus, hatchery 
Sanitation and management are critical in poultry 
health. Because of the continuous hatching program 
an the poultry industry, the air in the hatchery may 
be contaminated with various organisms, although most 
of these are non-pathogenic, a few pathogenic species 
can cause serious health problems for the broilers. 
This is particularly true for some respiratory 
diseases. For instance, the time between hatching 
and transporting of chicks out of the natchery may 
not seem important, but in reality it is critical. 
The longer these birds stay in a hatchery, the higher 
the possibility of their exposure to pathogens such 
aS Aspergillus species. Since the birds do not 
receive water and feed in the hatcher, to hold the 
birds in the hatchery for an extended period of time 
Will cause dehydration. Low temperatures or drafts 
in the hatchery can also add stress to these birds, 


Those stressed usually have more violent reactions to 


Ze 


| vaccinations and are more susceptible to disease 

agents at the onset or in their later life. In the 

case of Aspergillus spp., a few spores inhaled by the 

baby chicks, while they are extremely susceptible, 

i will most likely result in high condemnation at 
processing because of airsacculitis (9, 10). 

~ Brooding- in order to reduce the high cost of 
heating, many poultry farmers resort to partial house 


brooding instead of the traditional method of 





whole-house brooding. Some farmers prefer to block 
off a large segment of the poultry house and use the 


brooders in a small restricted area. In doing so, 


large numbers of baby chicks are crowded together, 





thus, increasing the density of the population (10). 
Other farmers prefer to turn on every other brooder 
in the poultry house; though the birds have the same 
amount of space, they are usually crowded together 
under the brooder. This in turn also increases the 
density of the population. Therefore, the type of 
brooding (whole-house or partial house) may prove to 
be a potential indicator in predicting condemnation, 
since crowded conditions created by partial—-house 
brooding may result in higher mortality, increased 


feed cost and respiratory disease conditions that 























lead to higher condemnation rates. 
- Housing and ventilation- in recent years, houses 


equipped with 3 or 5 cubic feet per minute (CFM) fans 
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have become fairly popular in this part of the 
country and it is anticipated that many more will be 
built in the near future. Because these houses are 
closed, with continuous rapid air movement, once a 
disease agent is present in a house, it will likely 
Spread rapidly. It is thought that birds reared in 
this type of house may have higher respiratory 
infection rates than those reared in conventional 
poultry houses. Thus, comparing condemnation records 
of flocks reared in conventional houses to those 
reared in the mechanically ventilated houses with a 
Similar method of management, may prove to be useful 
in predicting condemnations in the processing plant. 
Although not correlated witn condemnations due to 
diseases, housing environments and seasons have been 
Shown to lead to increased down-grading of broilers 
(12, 13, 14). Another vital aspect is that of the 
density of birds per given area which may affect the 
performance of broilers. Studies have shown that 
there is greater bacterial build up when birds are 
housed at a density of 0.42m> as compared to 
0.84m2/bird (15, 16, 17). Therefore, densities of 
bird population could be a useful indicator of 
population health. 

Feed consumption and feed conversion- a steady 
increase in feed consumption in young broilers is an 


indication of a nealthy flock. A sudden change in 














feed consumption may indicate disease problems or 
environmental changes. Thus, monitoring feed 
consumption of broilers throughout the growing period 
is a potentially useful indicator for predicting the 
condemnation rate of the flock. Other factors such 
as mold contaminated feed and their adverse effects 
on the performance of birds, are of particular 
concern since moldy corn and other concentrates, 
including ground feeds often carry mycotoxins. The 
effect of mycotoxins have been documented previously 
(18-23). Tnerefore, the type and source of feed, and 
its utilization could be useful indicators. Poultry 
are among the most efficient food producing animals. 
Flocks with an excellent health status usually attain 
feed conversion rates of two pounds of feed to one 
pound of meat. However, during stress and disease 
conditions, feed conversions usually are much poorer. 
Since feed efficiency is directly relatcd to the 
health status of flocks, its correlation with 
condemnation rates may be crucial. 

Morbidity and Mortality- although these may not 
necessarily reflect the rate of condemnation, the 
type and severity of the disease and the time between 
occurence of disease condition and processing, may be 
useful indicators of population health, In some 
eases, a flock of broilers may be affected with 


coccidiosis at an early age, and exhibit morbidity 
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and some mortality, but such an infection usually 
does not affect condemnation rate. On the other 
hand, morbidity and mortality caused by 
airsacculitis, Marek's disease, synovitis and 
Newcastle disease and infectious bronchitis 
particularly at 5-6 weeks of age, could lead to an 
increased condemnation rate. The relationship 
bewteen mortality and other conditions causing 
economic losses and types of housing have been 
studied previously. The results indicated that an 
automatic curtain house had the least mortality rate 
and gave the best returns to land, labor and 
management (24). Therefore, utilizing morbidity and 
mortality records in conjunction with necropsy 
records and other data, may serve as a reliable 
indicator of population health. 

Vaccination- all commerical poultry firms have their 
own vaccination programs for broilers, some of which 
are more effective in disease prevention than others. 
Variations in vaccination programs could make a 
difference in condemnation. The vaccines commonly 
employed are for the prevention of Newcastle disease, 
infectious bronchitis, infectious laryngotracheitis, 
fowl pox, Marek's disease, infectious bursal disease, 
ete. The age of birds at vaccination, use of one or 
multiple vaccines concurrently, use of booster 


vaccinations, or vaccines made by different 


26 



































27 


companies, may have important effects on the health 
and eventually affecting condemnation rate of the 
flock (25-31). Therefore, the history of 
vaccinations may prove to be a useful indicator. 
Medications— most commerical poultry firms institute 
their own health programs, to provide immediate 
treatment, when disease conditions are diagnosed. 

For example, there is definitely a difference in 
condemnation due to airsacculitis between birds 
treated with medication and those that were not 
treated during a mycoplasm infection. A prophylactic 
medication program specifically aimed at coccidiosis 
is routinely administered to broilers in the U.S. 
(11, 31). Use of various coccidiostats and other 
prophylactic and/or therapeutic medications could be 
useful indicators of flock health. 

Weather and seasonal differences— most commerical 
poultry houses are well equipped with heating units, 
so that the temperature inside poultry houses remains 
constant during winter months. During summer months, 
temperatures in chicken houses may increase to 110°F 
or even higher, particularly in the Southeast and 
Southwest. In the Fall and Spring, day and night 
time temperatures fluctuate drastically, (difference 
of 30-40°F) as a result of these changes birds are 
under stress and respiratory problems may become 


serious. Temperature and humidity are important 


factors which affect hatchability and chick quality 
(11), and could eventually reflect in increased 
condemnations. In order to study the weather and 
seasonal effect on condemnation, flock performance 
data collected during different months of the year 


could be useful. 


= Inspection procedures and condemnations —- the need to 


reexamine current poultry inspection procedures with 
the goal of seeking more efficient and less costly 
alternatives has recently been addressed (33-40). 
Currently condemnations of whole birds in FSIS 
inspected poultry slaughter plants are grouped by 
cause into eleven categories viz. Tuberculosis, 
leukosis, septicemia-toxemia, airsacculitis, 
Synovitis, tumors, bruises, cadavers, contamination, 
overscald and others. The first six of these are 
disease categories and amongst these conditions, 
incidence of condemnation due to septicemia- 
toxemia has been the highest in the past few years. 
For example, in 1978, a total of 38,799,191 young 
chickens were condemned for various reasons and of 
this total, 13,589,943 or 35% were condemned for 
Septicemia-toxemia, and in 1979, about 41% of the 
condemnations were due to the same category (41). 
The increase in numbers of birds condemned for 
septicemia-toxemia is a significant point and may 


very well mean that this category of production loss 
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needs further study (42). 

- Sanitation- is the mark of good management and is an 
area that requires close scrutiny. The economics of 
today's poultry production mandates high bird 
densities in commerical operations, which have led to 
other problems, High bird densities increase the 
probability of disease spreading via contact and 
airborne organisms from contaminated dust and air. 
These basically are sanitation problems. 

It was reported that coliforms represented aobut 10% of 
the microbial population of poultry litter (43). Reports 
also indicate that the housing of birds at 0.42m>, as 
compared to the density of 0.84m>/bird, led to significantly 
greater bacterial build-up of both aerobic and anaerobic 
bacteria (15), while dust was shown to be responsbile for 
establishment of the microflora in fresh litter in poultry 
houses (44). Toxigenic molds, such as Asperigillus, 
Penicillium, Petriella and Scopulariopsis, have also been 
isolated from poultry litter (18). 

The concern for sanitation in poultry production and 
processing is two sided. Firstly, to prevent the spread of 
infections in poultry establishments, and ultimately to 
protect humans against such organisms as Salmonella, 
Staphylococcus and Clostridium botulinum (45). Since the 
air and dust, birds, insects and rodents, equipments and 
personnel could serve as sources for microbial flora, proper 


sanitation (washing, disinfecting, etc.) and close 


333% 


Supervision of personnel are crucial in determining poultry 
health. 

With respect to sanitation, one very useful method of 
surveillance in use in hatcheries and broiler houses is air 
exposure samplings of these sites to monitor the population 
size and types of organisms found in the environment. Its 
usefulness has been documented. In one study, using a New 
Brunswick STA 200 microbiological air sampler in poultry 
houses, Bacillus, Micrococcus, Proteus, Pseudomonas, 
Staphylococcus, four species of Clostridia and nine genera 
of molds were identified. The author emphasized the need 
for vigorous sanitation practices, especially in houses with 
high bird densities, since the numbers of airborne 
microorganisms in such houses were high (46). The inter- 
relationship between dust and bacterial contamination and 
type of houses as well as the air contamination in 
hatcheries being associated with early chick mortality has 
been documented (47, 48). 

Interrelationships and Interactions Between Factors 
Affecting Poultry Health: Although the major determinants 
of poultry health nave been presented briefly, the list of 
the multiple variables for consideration is rather extensive 
as will be shown later. The essential point is to recognize 
that each of these variables are interrelated biologically 
to directly or indirectly affect poultry condemnations in 
processing plants. With this in mind and with the goal 


being the isolation of those variables which will be useful 
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predictors of condemnation rates, the need to systematically 
examine the inter- 

relationships and interactions of all potentially useful 
predictor variables, using the systems approach, was an 


obvious and logical choice. 


FLOCK TESTING — IS IT A RATIONAL ADJUNCT TO POULTRY INSPECTION? 


In closing this section on the background preparation 


required to grasp the problem solving task specified in the 


objectives of this contract, it is important to interrelate the 


literature review presented so far, to the overall goal of 


predicing condemnation rates. 


els 


The Assumptions: It is vital to briefly touch on the idea 
of why flock testing is a rational adjunct (conceptually to 
poultry inspection by referring to figure 3 . The diagram 
presents a summarized list of factors which could affect 
condemnation rates in processing plants. Some of these 
factors are related to processing errors e.g. scalding, 
contamination, improper processing such as bruising, missing 
viscera, etc. Other factors are related to possible errors 
during inspection and record keeping, which are present in 
any human task e.g. experience, fatigue, rapid line speed, 
CLG. 

One major assumption is that due to the uniformity of 
basic training for all FSIS inspectors, the errors that are 
human related would be comparable among processing plants so 
that if sample data are collected from various plants, this 


should not introduce major problems. If so, other than 
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processing plant errors, the major determinant of 


condemnations is disease (figure 3). 


The Rationaie: Since it is highly unlikely that diseases 


Suddenly arise in transit (although existing infections 
could be exacerbated as a result of transit stress) or while 
birds are being processed, it is only logical to think that 
disease conditions that lead to condemnations during 
inspection must have had some time to develop and be 
recognized as lesions, and in that case, the infection must 
nave taken place while the flock was still on the farm. 
This, of course, necessitated a critical examination of the 


farm ecosystem to isolate and identify variables which would 


be useful in predicting condemnation rates; this was the 
task of this research design. 

The above reasoning was supported in another 
preliminary diagrammatic model (figure 4). For a bird to be 
condemned due to diseases, in general, the following must 
hold: 

a. The bird was already diseased (D) on farm or at 
time ofevshipping. let us caliethis: P(D) = 
prevalence rate or proportion infected on farm, 

b. The bird is showing lesions (L) for the disease. 
P(L/D) = probability of finding lesions given that 
the bird is diseased. 

e¢. Tne inspector (I) picks out the lesion correctly. 
P(I/L) = probability that the inspector correctly 


picks out the diseased bird, given that the bird 
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exhibits lesions. 
Tne figure and the above listed criteria showed that the 
probility of condemnation was mainly a funetion of whether a 
bird was diseased P(D), the inspectors’ experience accuracy 
in picking out birds that were diseased with high 
Sensitivity (finding and condemning birds that are diseased) 
and high specificity (not condemning birds that are 
healthily. If we assume that the inspector has the above 
characteristics due to comparable USDA training, then the 
emphasis will be on identifying diseased birds P(D) on the 
farm. Note that the assumption that inspectors have 
comparable knowledge may be questionable. 

The function for the condemnation rate (CR) or 

probability of condemnation P(C) is: 

GReemftecUjeer tL), ue Ul/L). 
This determines the joint probability of condemnation at the 
plant. Note that if the P(D) is low, the CR will be low and 
this is the area of emphasis in flock testing. 


Die THE RESEARCH DESIGN FOR FLOCK TESTING 





5.1. Study Format-The Poultry Systems Approach: The study design 
utilized a systems concept, whereby a detailed evaluation of 
the determinants of population health were performed by 
systematically identifying and interrelating all variables 
which nave biologically explainable influence on 
condemnation rates due to diseases. Environmental, manage- 
mental, as well as host and agent factors which may affect 


the incidence of poultry diseases and the condemnation rate 


Dees 


Deceu 


at slaughter plants, were scrutinized via casual diagrams and 
multivariate epidemiologic models designed for possible 
predictive and/or discriminative purposes. Since poultry 
(broiler) production involves a sequentially predetermined 
set of factors which include breeder flock sources, 

hatching, brooding, growing and processing, the 
identifications of all potentially useful predictors of the 
condemnation rate at slaughter were integral parts of this 
research; designed to take the total poultry production 


ecosystem into consideration. 


Procedure for Identification of Variables: The complete 


list of original variables for which data were sought are 
provided in Appendicies 1 and 2. 

The hatchery subsystem: The first task was to identify the 
variables which operate in the hatcheries and had a 
potential causal influence on condemnation rates (figure 5). 
In this and subsequent diagrams, dotted lines indicate the 
paths of causal factors which may influence the population 
(state variables) to become infected, e.g. sanitation 
affects the population and type of microflora in a hatcher, 
Which may in turn affect the rate of infection in the 
hatched chicks. The bold lines indicate material transfer 
from one state to another, e.g. eggs obtained from breeder 
flock sources are set in incubators. The fertile ones are 
hatched, then transferred to brooder houses. The circles 
represent auxillary variables which have an influence on 


State variables represented by rectangles, e.g. in the 
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hatchery there are two state variables representing fertile 

and infertile eggs in the setter. From this point on, those 

fertile eggs which hatched transit to two states, viz 
infected or uninfected chicks. The infertile and unhatched 
ones are discarded. The infection in the hatchery may arise 
via genetic or egg—borne transmission (from breeder flock), 
or as a result of exposure to various disease factors and 

infectious agents in the hatchery as indicated in figure 5. 

Note that the transitions from fertile to hatched are 

sequenctial events that take place in the hatchery, and 

these transitions could have rates such as fertility, 
hatchability. 

The most important areas to examine are the variables 
(circles and the dotted lines) that originate from or lead to such 
areas; these provide the clues to the causal influences on 
population health which eventually influence condemnation rates 
further down the line. The variables identified via this figure 
were: 

- Hatchery Source (Variable codes Xooe Xege Kegs Xo60° X56! Xogo8 
X563)° There were 8 hatcheries which provide baby chicks to the 
Study farms in this project. Each of these were separate 
entities with different management practice so that 
Variabilities do exist in the source of hatched chicks for 
broiler production. Since the hatchery source was a categorical 
or nominal variable, it was not possible to quantitate it 
directly. There will be others like this, so it is important to 


note how one establishes an indicator or dummy variable for such 
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cases. Using hatchery 1 (Xe) as an example: 
Xoo = 1 if chicks came from hatchery 1, 

0 if chicks did not come from hatchery 1, 
Thus any time the source for baby chicks was hatchery 1, a value 
of 1 is given to variable Xe7 and O assigned to others, In 
general, when establishing such indicator variables, the rule is 
to use one less than the actual number of nominal variables 
(49). In this case then, since baby chicks were obtained from 
any one or two of the hatcheries and since there were eight 
hatcheries, the total number of dummy variables established were 
seven. 
Type of incubator (Xe os Xe4- There were two major types of 
incubators in use, viz, Chickmaster and Big J, with a few 
"other" types. Again, since this was nominal data, three 
categories (Chickmaster, Big J., others) were set up and two 
indictor variables were created. 
Ventilation is an important area, and what is taken into 
consideration is the type of ventilators Xogs X545 Zy) and the 
number of ventilators (X55). Three categories of ventilator 
types (3 CFM, 5.CFM, and others) were established and two dummy 
variables were created. Location of air intakes (Xne6s X53 Xo» 
Z.) and exhausts (Xoeas XogXogs Ze) were all nominal variables 
with 5 levels (directions) for each. The directions considered 
were: roof, east northeast, north northwest, west southwest, 
and south southeast. An appropriate number of dummy variables 


were then established. Note that ventilation is important 


either in introducting microbial agents from outside, or 
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spreading the ones already in the hatchery or poultry house. 

The information on ventilation was obtained from observation of 
hatcheries. 

Sanitation assessment of hatchery (Xo4)> This was another type 
of variable which was difficult to quantify directly. It was a 
composite index made up of various categories. To quantitate 
this data in a somewhat objective manner, and to make it useful, 
Six criteria were established. Each hatchery was graded on a 
seale of 1-3 with one being "good", 2 representing “average” and 
3 for "poor" or "unsatisfactory". The criteria were: offensive 
odor, dustiness, cobwebs on walls/ceilings, freedom from garbage 
and refuse, availability of waste disposal within easy reach or 
Sight, and unwashed (unsanitary) equipments/areas. Based on a 
grade for each of the six criteria, a final “average” grade 
referred to as the sanitation idex was computed. Thus, this 
Variable was one which was ranked based on specified criteria. 
Note that ranking is another established quantitative technique 
in handling nonnumerical data as in the case of nominal 
variables (49). 

A set of directly quantifiable variables for which numerical 
vaiues were available for, were also identified. Only the 
variable codes were given here and one could easily refer to 
data collection form 1 and obtain a description of what these 
Variables were. These were X60 = Xoo and Kage The information 
about these variables were gathered from regular visits to 


hatcheries. 


- Another set of variables closely associated with hatchery data 
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were those pertaining to breeder flock sources and major strain 
of birds. Four levels for the major strain of birds (and the 
breeder source) in Alabama were assigned for Corbett, Arbor 
Acres, Peterson, and "others". This was based on prior 
information on the major groups of breeding stock in the state. 
Again, breed/strain was nominal data and appropriate dummy 
variables (X., Xg» X19) were set up according to the method 
described earlier. By establishing these variables for the 
strain of birds, it also automatically took care of the breeder 
flock sources since these sources basically specialized in 
producing only certain strains of birds. Other quantifiable 
data, such as Xay and X45 pertaining to the number of strains 
and breeder flock sources, were handled directly. Such data 
were obtained from visits of hatcheries. 

The broiler house subsystem-brooding and growing- are the 
subsequent stages which follow hatching. Again, the birds are 
exposed to various risk factors in such a setting and 
epidemiologically relevant factors were identified with the aid 
of a systems diagram (Figure 6). The keys to the diagram were 
presented under hatching data. In this case, there were two 
major state variables, an uninfected and an infected sub- 
population constituting the total broiler population in a given 
house or farm. The infected subpopulation arises as a result of 
exposure to risk factors while the birds are in the broiler 
houses. Of course the converse, that is transition from 
infected to recovered (uninfected), occurs as a result of 


therapeutic medications, while vaccinations, prophylactic 


4] 


42 





O31 
Ost 

$40}284 
WWD 


SyRW}}>- 
weds uss} 





L $201922 {SUSI} 









F wayzepndod 
i, susiuetsg 47 





; : “ - a 
Bujssasasd jo Buyoejd joy % / 


2eak 40 Arg 





Ge 
























j NOLIWWidOd C3LINIAN 
) 
i, ne / \ \ 
Y RUPUD VO} WINGO 39 #RJ- i 
4 i | ‘ 
I ! aSeyasnou Ag azys uorRin dod 3SHOH INOUE ' \ 
| — 
4? 1 H ‘ 
at ‘ 
j 4 
4 ’ j \ 
/ ! ‘ 4 ' 
pet \ d t 
sf ‘ j : 
/ ' * gia 
t] \ $9jpNOI- °939 @IUnO3- 
' SebEsop- here pdky> | 
4 yoyyesrip- voya2sn9- bujpo0sg 
j \ acy 3 
O\ SUD YOUJIWA / 
: alo neet . wat \sap2yhyzosg 
(43103403) N ? 
91 199/p390 S : 7 
NS pid 





$y 0° 

buj2e0se 
B1NpB4IS- 
BIAQ- 


Gu] 0}8m/5U]p984 


(woqy>OyU}) WAHLSASGNS ASNOH UATIOUG AHL 9 °3Ta 


Aypiqow 





medications and other disease preventive programs maintain the 
portion of the uninfected subpopulation in that state, Again, 
the significant point was the identification of all potential 
variables which may infiuence the subpopulation of uninfected, 
to transit to the state of being infected. The variables so 
identified were: 

- Population descriptive variables - included in this category 
were initial flock size (X3) and the flock size at time of 
processing (X55) and estimated proportion of under sized birds 
in the population (X7). These variables were directly 
quantifiable and were gathered during visits to the growers, 

- Housing - the type and area (density) available for the flock 
Was considered. In terms of house type, three categories were 
established for type 1 (conventional/fan = Xig)s type 2 
(conventional/no fan = X45 and type 3 (others). Since these 
were nominal variables, two dummy variables were set up 
accordingly. The second aspect of housing was population 
density per square foor (X19), which was directly quantifiable 
from available data. 


- Ventilation — type and the direction that the intake and 





exhausts face is important in introducing microbial agents from 
external sources and/or spreading them within poultry houses. 
The air itake variables considered were (Xone XagXaqs Zoo: 
The corresponding air exhaust direction variables were (Xo65° 
Xe47Xg29 Z,,)- The respective directions were roof, east north 
east, north north west, west south west, and south south east. 


The appropriate data were obtained from field observation of the 


broiler hourses. 

- Sanitation index for broiler farms (X4 96) was computed in a 
comparable manner as discussed under hatchery sanitation index. 
Tne assessment of sanitation was performed by the same two 
individuals who visisted and graded the hatcheries. This way, a 
comparable type of data with less variability amongst observers 
was generated. 

- Litter is related to housing, and inciuded here were variables 


such as the type of litter (X for shavings or sawdust), 


189° 
litter change for subsequent flocks (Xigoe partial or complete), 
depth of litter (Xi gq) interval between litter changes (X194) 
and number of batches of birds on the same litter (Xy9y)+ The 
last three variables were directly quantifiable, while the first 
two were dichotomous and one dummy variable represented each of 
those. Other variables, such as procedures for house cleaning 
after litter removal and methods of removing litter, were also 
considered. Such data were obtained from growers. 

- Management assessment criteria - as in the case for sanitation 
index, relatively objective criteria for ranking purposes were 
established for farm record keeping, cooperative attitude of 
farm personnel and overall farm organization and management. By 
grading on a scale of 1-3, an index of record keeping (X5Q9) 


cooperation index (X544) and management index (X,,.) were 


ele 
established. These criteria are indirect measures and are 
somewhat subjective and thus limited, but data were gathered for 


possible descriptive use. Other management concerns such as 


farm security, and rodent control programs were also obtained 
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when available. 


- Brooding variables - directly quantifiable data such as day of 
year of placing with day one for the year beginning January 1 
(Xee). number of chicks dead on arrival (Xpe). length of time 
(hours) for preheating of broiler house (Xe). partial vs whole 
house brooding (Xgo)s age (days) when brooder guards were 
removed (Xoo) and flock size per house (X94) were obtained from 
growers by requesting such information during scheduled farm 
visits. 

~ Macro and microclimatic data — more specifically, this referred 
to the macro temperatures (average maximum, average, average 
minimum temperatures) -— (X5957Xogy) total precipitation (X5q5? 
and number of rainy days (Zr) per month, these affect the micro~ 
climate (Xg3-Xg¢) in the broiler house. This information was 
obtained from the National Climatological Data Center (ref.). 
The data for the nearest weather station where study farms were 
located, was used for this purpose. It was not possible to 
obtain microclimatic data from broiler houses, in a regular 
basis, since growers neither keep such records, nor would they 
cooperate in keeping daily records on such data. 

- ason - although season is a function of the climatic variables 
listed above and the study design takes this into consideration, 
season in itself was also used as an independent index with 
possible predictive power. Four levels of season Cfalaewinter,: 
spring and summer) were designated by three dummy variables. 
They represented fall (X56) winter (X597) and spring (X59): 


Thus, the influence of season on a given brood of birds could be 


captured in an approach which has been utilized in similar 
fashions in other fields (50). 

Morbidity indicators —- as given in figure 5, the number of 
infected birds divided by the total population on a given farm 
or house represents the morbidity rate, which directly leads to 
infected birds being slaughtered. However, the morbdidity rate 
of the subpopulation of birds that may have been infected were 
difficult to identify. Data of morbidity were unavailable from 
growers and company records. Thus, various indirect measures of 
morbidity were devised to assess this vital information, The 


indirect indicators used in this study were: 


mortality data obtained from existing farm records 


~- serological profiling 


necropsy profiling and 


feed utilization 

Except for mortality, data of the rest of the above indicators 
were gathered from field for laboratory studies by the research 
team, and these represented a major portion of the field work. 
These last three items will be discussed in greater detail 
further into the report. With respect to farm mortality data, 
reasonably accurate weekly mortality records were available 
although the cause of mortality by specific disease were 
unavailable. Mortality rate on farm (X29) and more 
Specifically, age specific mortality rate at 4 weeks of age 
(X29) and at 7 weeks of age (Xy9) were computed using the 


formula: 


X = total number of birds dead by 4 or 7 weeks on farm (house) 
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average population at risk on farm (house) during same 
period 


Where X is X29 or Xno- 

In addition, specific mortality rate in a house (X347X39) X59) 
was also utilized as a predictive variable. It is important to 
note that mortality data in itself is not a very useful 
indicator for condemnation rate since the birds are already dead 
and never reach the processing stage, but the assumption is that 
other than mortality due to non-disease causes, those deaths 
that occur as a result of infection, may be correlated with 
morbidity in the population. The relationship is not so much 
with morbidity, but with case fatality. This is so because 
there could be a chronic infection which may not result in the 
death of the bird but may eventually be the cause for 
condemnations at the time of processing. 

Microbial profile data - It is vital to see how the research 
design utilized indirect measurements of the population 
morbidity rate to correlate with condemnation rates. If one 
examined both areas carefully, it can be recognized that many 
variables, directly or indirectly, affect either the population 
and/or the types of microflora in the hatchery or broiler house 
subsystems which pose risks of infection in susceptible birds. 
For example, in the hatchery, factors pertaining to sanitation, 
ventilation, climate and type of incubator allaffect the 
mircobial flora, which in turn poses the risk of infection for 
newly hatched chicks. In the broiler sybsystem of the diagram, 


such factors as rodent control, sanitation, litter, ventilation, 
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housing and climate determine the type and population of micro- 
organisms in broiler houses. 

The high density of birds in large poultry production units 
increases the possibility of disease spreading through airborne 
microroganisms, dust and air contaminants. Therefore environ- 
mental quality control and the collection of information 
concerning airborne bacteria and fungi in poultry establishments 
becomes crucial. Data concerning the overall population of 
organisms and of their type, specifically on their pathogenicity 


and/or toxin producing potential, are useful (18, 46). 





The research design emphasized the need for 
obtaining a microbial profile of the flora in 
hatcheries and poultry houses; these were 
respectively referred to as variables X00 to Xiao 
for the population size of bacteria and fungus in 
randomly selected air samples from poultry 
establishments. Additionally, because of the need 
to identify the types of microbial organisms in the 
air samples, the variables designated X 103 - X07? 
were created. Although it was recognized that 
microbial profiling (quantitating the population and 
identifying the types of organisms) was essential in 
establishing the link of factors that affect poultry 
health (Figures 5, 6) due to limited financial 
supports and lack of semiautomated equipment, to 
handie the large column of samples collected, 
bacterial and fungal agents picked up from 
hatcheries and poultry farms were identified only 
for a few study farms. 

Medication ~ Veterinary health care, (therapeutic 
and/or prophylactic medications), plays an 

important role in determining morbidity and 
mortality rates in poultry (Figure 6). The 
specific variables identified with potential 
influence on morbidity and mortality were types of 


medication, duration of administration of 


medication, dosages and routes of administration. 
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Both therapeutic and prophylactic medications 
contain the variables named X438 - X63 and Zan - 
Lage 
Vaccination — in the systems diagram (figure 6), 
prophylactic medications and vaccinations were 
placed together in one category, because both are 
disease preventive programs intended to minimize 
the transition path from an uninfected state to an 
infected state. In terms of vaccination data, the 
variables for this study include vaccine 
manufacturers (X69 _ Xaqys Zig)» route of 
vaccination (Xizoe X73" Zog)s vaccination schedule 


(X Z,,) type of vaccine (Xigye Xigos Zz 
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and degree of adverse reaction to the vaccine (Xigy 
- X1g5)° Most of these variables were categorical 
and appropriate indicator variables were set up to 
quantitate them. 

- Feed utilization - aS one of the indirect health 
status indicators in a population, feed utilization 
occupied an important part in this study (figure 
6). Feed utilization of birds is affected by such 
variables as ambient temperature, feeding/watering 
practices, (type, spacing, scheduled vs ad libitum 
feeding) administration of additives and above all, 
morbidity or prevalence of infections in the flock. 
Conversely, feed utilization affects rate of weight 


gain, which in turn reflects the health of the 


50 





flock, including morbidity. 

Because of the importance of feed utilization data 
in providing clues to the morbidity rate on a farm, 
a detailed listing of variables from type of feed 
(Xi 46° Xa47s feeding Oi 49? or watering (X49) 
Space to a procedure for collecting data 
establishing weekly feed utilization and the rate 
of weight gain, was designed (Xj 35 - X435)+ Some 
of the variables were nominal in nature, while 
others were directly quantifiable. 

5.2.3 The transit and processing system: The third stage in the 
sequence of poultry production phases which were evaluated 
Via a systems approach were factors affecting condemnation 
of poultry due to transit or those that arose while birds 
were being processed (figure 7). This is similar to 
figure 4, except that in-transit factors have been added 
and the figure expanded further to focus on condemnations. 
Although it has been pointed out earlier that an infection 
must have necessarily arisen while birds were on the farm 
for it to have enough time to reach the stage of lesion 
development, other factors which may have some influence 
on condemnations such as transit stress including: 
distance from farm to plant (X53) number of hours of 
driving from farm to plant (Kooy) drivers and trucks 
involved (X05 - Xoo) types of coops used (X59) and 
population density per coop (X520) and per truck (X534)) 


method of cooling birds in processing plant (X53) and 
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macroclimate on date of processing (X539 - X 40> were also 


2 
considered. Some of these data were nominal, while others 
were directly quantifiable and the methods used to 
establish values for the nominal data were the same as 
given in earlier sections. The most important aspect of 
this last segment was that data on condemnations due to 
diseases and other causes (Y,, - Yuq)» which served as the 
dependent variables, were also collected at this time, 
Such data were provided by management of the plant upon 
completion of the condemnation report by the inspector in 
charge. The same categories for causes of condemnation 
due to diseases (FSIS reports) were used in the data 
collection forms (appendix 2). 

In addition to condemnations due to diseases, data on 
total number of pounds (lbs) of parts trimmed and 
non-disease condemnations were also obtained. Additional 
data on condemnation for birds on all farms processed in 
that plant on the day that birds on study farm were 
processed, obtained for comparative purposes. from these 
data, the condemnation rate due to disease for the study 


farm (Y.) was computed using the formula: 
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Y, = Total number of study farm birds (whole carcass) 
condemned due to diseases 


Total number of study farm birds inspected 
Disease specific condemnation rates of birds on study 
farms were also computed for the following diseases: 


leucosis (YQ). septicemia—-toxemia (Y.,), airsacculitis 


10 


(Yow) synovitis (Y-.), tumors (ls.) ,eand souners (Y41)- 
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These disease specific rates were computed from formula: 


Y = Number of study farm birds condemned due to specific 
disease 


Total number of study farm birds inspected 


Where Y = disease specific rate. 

Note that since these were the dependent variables, the 
code names begin with Y while all other possible predictor 
or explanatory variable names were coded as X's, One 
other variable of importance was the condemnation rate for 
the study farm at one previous brood, referred to as Toes 
(Y at t-1 period, or condemnation rate at one previous 
brood), which was used as a historical reference variable 
for predictive purposes. This is based on the rationale 
that the performance of a poultry grower may in fact be 
judged by how one performed or managed at previous times, 
i.e. based on experience or historical performance. 

Since the systematic examination and identification of 


variables with potential predictive value have now reached 


the end point of the poultry production process, a more 
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holistic systems diagram of the poultry exosystem and its 
determinants as there affecting condemnation of poultry is 
presented (figure 8). The various segments of this 
diagram have already been discussed, so what is shown here 
Should be self explanatory in providing one with the total 
conceptual and design framework for analyzing the data 
generated from such a model. 

NOTE: The complete list of original variables identified 
from the hatchery, broiler house, transit, and processing 
subsystems are in appendices 1 and 2. 

Study Unit: Based on pilot studies to determine where 
adequate data could be collected from an individual 
poultry farm was initially considered the study unit. 
However, the original design specified data collection 
from a farm on a house by house basis too whenever 
possible. The major problem to collecting the data on a 
house basis was that data for the dependent variable 
(condemnation) on a house basis at the slaughtering plant 
Was unavailable. Later, by prior arrangements with the 
Management of each cooperating firm and inspector in 
charge of the processing plants, birds were processed and 
data recorded on a house basis. Therefore, both the house 
and farm served as two experimental units. 

Sample selection: Having identified what the experimental 
unit was, five poultry companies in Alabama was selected 
for the study. Originally, a letter was sent out to nine 


firms in Alabama, Only four firms responded to the 
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request for cooperation in this study, and personal visits 
and discussions held with officials of these four firms. 
Only three firms were willing to cooperate in this study. 
However, with more assistance from Dr. Edgar and FSIS, two 
more firms eventually participated in the project. The 
selection criteria for experimental units (farm/house) 
from each of the five firms consisted of identifying 
cooperative growers within reasonable driving distance of 
each other. Thus, the question of cooperation and 
logistics was vital as opposed to the originally planned 
method of random sampling of farms. 

Stratification of sampling units: To ensure a 
representative cross section of farms, the study farms 
were stratified across three strata, of good, average and 
poor production performance based on the criteria for 
classifying a grower in one of these categories. The 
criteria for ranking growers varied from firm to firm as 
presented below: 

Procedures for ranking grower performance for firm 1: The 
growers were ranked both weekly and annually. For weekly 
ranking, all growers that had birds processed in a 
specific week were compared to one another via an index of 
performance criteria to be presented below. Annually, the 
farms were ranked on production costs. 

Tne fixed cost factors considered were: cost per baby 
chick, cost per ton of feed, cost of vaccination, de- 


beaking cost, and overhead and service costs. 
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The variable cost factors were: sanitation costs (wash 
water, cleaning area around chicken house), and medication 
costs (depends on type of drug and duration of usage). 

The total cost of producing the flock divided by the total 
pounds of birds processed resulted in the cost for a pound 
of meat which was the criteria for ranking growers as; 
good, average or poor. Growers with the highest 
production cost were ranked, "poor". Those with the least 
production cost were ranked, "good", and the rest were 
rated as average. 

Procedures for ranking grower performance for firm 2: The 


growers were ranked weekly and annually in a comparable 
manner as described for firm 1. 

The fixed cost factors considered were: cost per baby 
chick, and cost per ton of feed. 

The variable cost involves medication and other costs. 
Again, the cost items were summed and cost per pound of 
meat computed as in firm 1, Ranking was assigned 
accordingly. 

Procedures for ranking grower performance for firm 3: The 
growers in this frim were ranked weekly through a 
computerized system of accounting. There were three 
categories viz. above average, average, below average, 
The fixed cost factors considered were: baby chick cost, 
feed cost, and service and overhead cost. 

Tne variabie cost factors were: medication cost, and 


miscellaneous cost. 











Tne above cost items were summed to obtain the total cost. 
The firm then considered three indices: the average live 
weight that birds should attain was set at 3.85 pounds, 
attain the least amount of feed consumed to produce a 
pound of meat, and cost due to condemnation loss computed 
from the report of IIC. 

The farm production cost plus the loss due to 
condemnations were combined to derive the grower's cost. 
Tne grower's cost for a certain week were summed and this 
total cost was then divided by the number of growers to 
arrive at the average growers cost. If an individual 
grower's cost was higher than the average grower's cost, 
he/she would be ranked, “below average". If the cost was 
lower, the individual grower was ranked, “above average”. 
Procedures for ranking grower performance for firm 4: The 
growers were placed in one of three categories monthly 
through a computerized system. The formula for ranking 
growers is an adjusted feed conversion index which was 
derived from the average live weight of birds + a standard 
weight of 3.67 lbs. For example, if a grower produced an 
average weight of 4.12 ibs broilers and had a feed 
conversion index of 2.07, therefore, the adjusted feed 
conversion would be: 4.12 — 3.67 = 0.45 lbs. 

The above value was then divided by 5, an arbitrary fixed 
value of 0.09 lbs was otained. One then substracted 0.09 
lbs from the original feed conversion index, of 2.07 lbs., 


to arrive at the final adjusted feed conversion index for 
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the grower, which in this case would be 1.98 lbs. of feed 
for one 1b. of meat. In this ranking procedure, the 
lowest adjusted feed conversion index is the best while 
the highest corresponded to a poor rank. 

Procedures for ranking grower performance for firm 5: 
Growers were placed in one of the three categories in 
association with production costs. The fixed cost factors 
were: cost per chick and cost per ton of feed. 

The only variable cost factor was medication. 

The total production costs of per pound of marketable meat 
was computed for each grower after their birds were 
processed. The average production cost for all farms with 
birds processed in a certain week, was obtained by 
dividing the total production cost by total pounds of 
marketable meat after condemnation. If an individual 
grower's production cost was higher than average, a minus 
(-) factor would be assigned to him; if the production 
cost was lower than average, a plus (+) factor would be 
assigned. These cost figures were calecuated for each 
grower after each processing of each brood of birds and 
also annually. 

The annual production cost was derived by dividing the 
total production cost for the year by the number of 
broods. To encourage lower production cost, a prize was 
Ziven to the grower with the lowest production cost for 
the year. Using the criteria for the selection of study 


units from each firm, at least two growers were selected 
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from each strata. This decision was made to ensure that 

at least one grower would represent each strata in case 

some of them decided to withdraw from the study for 

reasons beyond the control of the research team. If all 
selected growers continued to stay in the project, this 
would expand and strengthen the data base. 

One major adjustment was made at the end of August, 1981- 
to reduce the total number of study farms per firm from 6~7 down 
to 3-4. This adjustment was necessitated as a result of a 
reduction in the number of student research assistants who were 
working on full time basis during the summer, but had to go back 
to school in early September. Thus, the practice of fielding 
two or three teams concurrently to collect data was not possible. 

Although the number of study farms were reduced to three 
or four, the criteria for selection of which farm to drop and 
which to retain was based on: 

—- maintaining representation of good, average and poor 

farm rankings as described earlier 

~ keeping the farms where the growers were most cooperative 

in providing detailed data regularly collected at four or 
seven weeks 

The final make up of the study farms consisted of four 
farms from firm 1(A), four farms from firm 2(B), three farms 
from firm 3(D), four farms from firm 4(E), and five farms from 
firm 5(F). Therefore, there were a total of 20 study farms 
which consisted of 40 houses. Although there were potentially 
40 sets of observations, from each experimental unit (house) 


during any one brood cycle, one firm (D) did not provide 
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condemnation data on a house basis. 


5.3.3 Sampling time frames: Two type were followed: 


a. Sequential or longitudinal follow-up which involved: 

-examining the health status of the population and its 
different determinants at various cross-sections in 
time for a group of cohorts. 

-collecting data on the appropriately identified 
multivariables from farms during early, middle and 
later stages of a growout. 

-collecting corresponding condemnation data on the 
Same cohorts from processing plants. 

b. Cross-sectional sampling: 

-conducted at one cross-section in time to provide a 
measure of prevalance of the specific determinants 
within the specified farm ecosystem, 

-multivariable data collected from the farm before 
birds were shipped to slaughter and then obtaining 
the corresponding condemnation data from the 
processing plant. 

In essence, these were specified sample data callection 
time frames from a hatchery, broiler operation and 
processing plant (figure 9). Using these sampling time 
frames, data were collected during one complete year (5 
brood cycles) to accomodate possible seasonal and other 
time related variables. 


5.3.4 Sample data collection: by definition, a cohort of birds 
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that were sequentially monitored to collect data from each 
of the three phases of broiler operation (hatchery, 
broiler farm, processing plant), was referred to as one 
brood cycle. It requires approximately 50 days to raise a 
bird to the desirable marketing weight as a broiler. Such 
a data collection time frame was projected to take place 
for a period of one year. The first visit involved the 
hatchery, with specific goals to: 

~- assess the sanitation index for the hatchery 

- obtain air samples from randomly selected sites 

in the hatchery as well as from the hatcher 
where chicks were scheduled to be sent to one of 
the study farms. 

- collect preexisting hatchery data 

Subsequent to collecting hatchery data, each study 
farm was visited during the first week that birds were 
placed on the farm. During this time, air samples were 
obtained for microbial profiling to provide information on 
the potential infection causing factors in the brooder 
houses. During the first week of life, baby chicks are 
most susceptible to infection, therefore the microbial 
sampling of farms was performed at that time. 

Each farm was then visited at two different times; 
one when a given brood of birds were about four weeks old 
(visit 1) and again when they were about six to seven 
weeks old (visit 2). Note that, as given in the sampling 


time frames and as indicated in figure 9, a cohort of 





birds were followed up from the hatchery through the 
growout period. 

4-Week-Data (Visit 1): The data concerning flock history 
and different variables were obtained on about the 28th 
day, following placement of the chicks in the broiler 


houses, by questioning the grower and studying the 


performance charts in the chicken houses. A questionnaire 
(Form 1) concerning the variables mentioned previously was 
filled out by the visiting research assistants. At the 
same time, ten chickens per house were randomly selected 
and bled. The sera were transported to Tuskegee Institute 


for serological tests. 


7-Week-Data (Visit 2): These data were usually obtained 


on about the 49th day, following placement of chicks in 
the broiler house. Form 1 was completed for this time 
period with all available information. 

The final phase of the sample data collection time 
frame was that of obtaining condemnation data of study 
flocks from the processing plants, and the information was 
stored in Form 2. Later on in the report, this task will 
be discussed in greater detail. In addition to the 
sampling time frames, figure 9 presents the type of data 
collection forms to use and where. For example, all farm 
data were filled out on form 1 while processing data were 
filled out on form 2, The design of the data capturing 
format in this study has been an important aspect of the 


research development phase and will be presented in the 
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pages to come. 
Types of Data 
Tnere were three main types of data in this study. 
Preexisting data - e.g. condemnation rate at processing 
(Y5), weekly mortality rate from records kept by the 
grower (X29)s etc. In such cases, the data were collected 
from existing records either on the farm or in the 
processing plants, and the quantitative vaiues of the 
Variables were incorporated with the corresponding set of 
data in a multivariate setting, both in terms of time and 


place. 


Field generated data -— two subcategories were considered 


here. 

a. Data generated from field observations, e.g. 
sanitation index (%186), management variables 
(X212), etc. Research assistants of this project 
assessed the respective cases and provided 
quantitative estimates (numerical values) for the 
variables under study. 


b. Data generated from field measurements —- e.g, th 


(99) 


weight of birds (X2), the proportion of undersized 
birds 0X7) ,..the-weekly, feed utulization, etc. thn 
tnese cases, actual measurements, e.g. counting 
undersized sbinds, in japsample,of 100 birds wel ening 
a few randomly selected birds, or measuring feed 
bins to estimate weekly feed utilization, were 


performed to generate data. An important part of 
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this task was data collection on weekly feed 

utilization; a detailed description of which follows: 
Feed utilization data collection procedure: Three farms 
from each cooperating firm (one good, one average and one 
poor) were monitored for feed utilization. These farms 
were visited once a week when possible. The level of the 
feea in each feed bin was estimated, compared to the level 
of the previous week. This provided the volume of 
displacement, which was converted into pounds of feed, If 
the level was higher than the previous week, it was 
assumed that additional feed must have been delivered. In 
some cases the growers kept the feed delivery tickets, 
from which the amount of feed delivered and the date of 
delivery were recorded for calculations at a later date. 
On one farm, the feed that was actually weighed to measure 
the amount of feed was dispensed at a particular time. 
Tnis again, varied with the size of the birds and the 
temperature inside the chicken house. 
Tne method for estimating the amount of feed in feed bins 
was indirect; the details of which are provided in 
appendix 3. Based on the method described, an estimate of 
the amount of feed consumed per bird per week was computed 


for a limited number of flocks with adequate information. 


5.4.3 Laboratory generated data: 


a. .oanple coltection for serological and» necrops 





profiling: An important aspect of the study design 


was the sampling of a representative number of birds 
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from the study farms for serological and necropsy 
examination. As referred to earlier, the study farms 
were selected by stratifying representation from good, 
average and poor production performance farm groups 
based on information provided to us by the repective 
cooperating firms. Although there was some stratifi- 
cation across three production strata, there was no 
random sampling of farms from each strata, 
The importance of these parameters is primarily, for 
providing indirect measures of morbidity. As pointed out 
earlier (figure 6, 7), since morbidity data are difficult 
to obtain, and since most infected birds are eventually 
condemned upon processing, some mechanism of obtaining 
prevaliance data on a given population of birds immediately 
before slaughtering could prove to be the most important 
piece of information in predicing condemnation. 

It was with this background that serological and 
necropsy profiles of a flock of birds was studied since 
both parameters could possibly be adequate indicators of 
morbiaity. 

Serum Samples Collection Procedure: From each of the 
Study units, blood samples were obtained at two stages 
during a brood eycle. The first stage samples were 
collected at four week of age from ten randomly selected 
birds per house, 

Bleeding Procedure: A lice hypodermic syringe with a 26 


guage 1/2" needle was used to collect blood from the 
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subclavian vein under the wing of the birds. The blood 
was then slowly forced out of the syringe into a 5ce test 
tube. A cap was placed on the tube to prevent 
contamination and spillage. The sampies were transported 
to Tuskegee Institute, where the serum was separated from 
blood celis by centrifugation (2,000 rpm for 10 minutes) 
anda stored at ait in a serum bank until serological 
analysis were performed. 

The second stage sampling of sera was done at the 
slaughter plant when the birds were being processed at 
approximately seven weeks of age. Blood samples were 
obtained from birds after the jugular veins were cut and 
were being bled out. Approximately four ml of blood was 
obtained from randomly selected birds, as described under 
the sampling procedure at processing plant. The blood 
samples were treated in a similar manner as in the first 
Stage sampling. Sera from additional birds of the same 
flock were stored in the serum bank of the poultry Disease 
Research Laboratory, and will serve es future references. 
Three indicator diseases were tested in our laboratories 
by using the hemagglutination/hemaggiutination inhibition 


(HA/HI) tests. These disease agents were Newcastle, 


Mycoplasma gallisepticm and Mycoplasma synoviae. A titer 


Ofsdc4O waseusedaas the cut off point.for, positivity— 
hefgativitysin the easeyof thermycoplasma, while a,titer jof 
1:80 was used as a cutoff point for Newcastle disease. 


Since broilers were vaccinated with Newcastle disease 
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vaccine at an early age, a low titer of antibodies usually 
persist for some time and therefore a higher level of 
titer was used as the cutoff point as an indicator of 
infection. 

During the first visit to the farm when the birds 
were about 4 weeks old, blood samples were collected from 
a randomly selected sample of 10 birds per house. From 
the processing plants, a predetermined sample size of 200 
blood samples from each farm were collected for 
serological testing. Concurrently, 290 uneviscerated 
poultry carcass samples were also obtained randomly for 
necropsy profile studies. 

Necropsy sampliing: In order to randomly select 200 birds 
for necropsy sampling and blood testing, informations 
concerning the time required to process the study flocks, 
the line speed and number of birds being processed were 
acquired before processing was initiated. For example, if 
it required one hour (3600 seconds) to process a flock, to 
obtain 200 birds during that 1 hour period, the research 
assistants must pick up a bird every 18 seconds. 

All necessary preparations were completed before the 
birds from study farms were processed. Insulated boxes 
With ice were placed at the sampling location. Blood 
sample collecting required one person, bird sample 
collecting required two people. A fourth person served as 
messanger to communicate between the back dock, sampling 


Station and tne inspection stations. 





In collecting the sample birds, instead of using 18 
second intervals, one used 15 seconds as the time 
interval, and picked the third bird that passed down the 
Chute. The sampling was conducted on line right after the 
hock joints were cut. Approximately 20 birds were placed 
in an insulated cardboard imbedded with ice. Each box was 
properly marked to identify the grower's name, house 
number and processing date. At the end of processing, the 
boxes containing sample birds were re-iced and loaded onto 
a truck and taken back to Tuskegee Institute, 

The necropsy profile (Xi99): The birds were usually 
necropsied and inspected a short time after arrival by the 
poultry pathologists associated with this project, Gross 
lesions were recorded in terms of specific conditions such 
as leukosis, septicemia-toxemia, airsacculitis, synovitis, 
tumors or others. Tne mere fact that the bird exhibitied 
lesion of disease, whether mild or severe, was sufficient 
evidence for that carcass to be considred a positive case. 
The extent of lesions and parts affected are vital to 
decision making as to whether or not to condemn a carcass, 
In this case though, our necropsy efforts were not 
intended to make that decision, but to find out the 
prevalance of lesions which may or may not lead to whole 
bird condemnation in slaughtered birds. 

Specific rates such as overall disease rate and 
specific disease rates were computed, using the total 


number of birds necropsied and inspected (approximately 
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200) as the denominator. The proportions ( or 
percentages) of diseased birds from each study farm were 
then utilized as one of the independent variables 


(specifically X ) to predict the overali condemnation 
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rate for that study farm at processing. Also, more 
specifically, specific disease proportions (X49 -Xy 45) 
were computed and utilized to predict the corresponding 
specific disease rates (Yg -Y 41) at processing. (Appendic 
1, 2). Data collection for X109 progressed through four 


brood cycles and had to be terminated thereafter. 
b. Microbial profiling: In order to determine: the effect 
microbial flora on condemnation of a specific flock, a 
limited microgial profiling study was conducted. Two 
aspects were involved in this process, Firstly, 
quantitative estimation of the microbial population was 
sought. The second aspect involved isolation and 
identification of microogranisms which occur in the 
hatchery and broiler house. Primary culture plates were 
prepared using standard bacteriologic technics, Initially 
three types of primary media were used; blood agar base 
for general bacterial culturing, MeConkey agar for 
differentiating gram negative bacteria, and corn meal agar 
for fungi culturing. Since some bacteria do not grow in 
blood agar base alone, blood was added to the blood agar 
base in subsequent culture attempts. 

The culture plates were placed in a bacterial 


incubator fork24 hourse= to! test formsteriiity, wands only 














sterile plates were used for bacterial culture, 

Microbial sampling of the hatcheries: Culture plates were 
exposed to the air at different sites in the hatchery. In 
order to determine an optimal air exposure time, intervals 
of one, two, five, ten, and fifteen minutes were tried. 

It was found that ten minutes exposure time was 
acceptable. After exposure, the culture plates were taped 
to avoid accidental exposure in transit and transported to 
the microbiology laboratory at Tuskegee Institute. 
Microbial sampling of broiler houses: Upon placement of 
baby chicks in the broiler house, microbial samplings were 
performed within the first week of their life, Three sets 
of culture plates were exposed in three locations (at both 
ends and in the middle) of the chicken house, for one 
minute each (the exposure time was predetermined as 
acceptable, by a pilot study). 

Bacterial and fungal colony counting: The samples from 
the hatchery and chicken house were incubated at 37 C for 


bacterial growth or were kept at room temperature for 


fungal growth. The number of bacterial and fungal colonies 


were counted, by using a Fisher Accu-Lite Colony Counter 
and recorded at 24, 48, and 72 hours. Since fungal 
colonies grow slowly, corn meal agar plates were kept for 
one week before discarding. 

identi ticavion.of abacter va: identification of -hacteria 
eollected from the study farms and hatcheries were not 


Piece avedemtitathe lacter part of this study due to lack 
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Ofsfunds:. 

Standard laboratory procedures were employed to 
iaentify the cultured bacterial colonies as deseribed 1 
"Diagnostic Procedures in Veterinary Microbiology" by 
Carter and Bergey's Manual of Determinative Bacteriology, 
8th edition (52). The steps involved in this aspect of 
the work are presented in greater detail in Appendix 5. 

In closing the section on research design, one should 
refer back to figure 8 to interrelate the significance of 
each variable, and grasp the holistic picture of its 
influence on the end result of predicting condemnation 
rates at processing. In that light, all: existing data 
collected from hatcheries, farms and processing plants, as 
well as newly generated data from our laboratories and 
field work (microbial, serological, necropsy profiles, and 
feed utilization), should be viewed as occupying a 
specified causally identifiable position vis-a-vis 


condemnation rate at processing. 
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Data Collection Format: The design of a poultry health 
data collection format has been a much more complex and 
time consuming process than originally anticipated. 
Basic requirements: The intention of developing such a 
form was: 

- to be comprehensive enough to accommodate all 
epidemiologically relevant data pertaining to poultry 
health and management within the framework of the 
poultry farm ecosystem, 

~ to enable the collection of quantitative or 
numerical data for analysis, using multivariate 
techniques, 


- to develop forms with the potential for serving as a 


Standardized poultry health/management reporting format 


So that the data gathered could be computer coded, 


stored, retrieved and analyzed if a national or regional 


program were made operational, 

- to provide relatively objective criteria to assess 
Sanitation, management, farm security, rodent control, 
and record keeping. 

The data base also allowed a retrospective review of data 

already collected, if such a need arose. For example, a 


format for disease outbreak investigation could be 


followed, and detailed data by house, age, strain, breeder 


source, hatchery, etc. could be examined to provide 


feedback to the grower and poultry firm if high 
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condemnation rates were encountered. Other additional 


data on health care and cost of the same have also been 


included, 


5.5.2 Types of forms: With these objectives in mind, and after 


5.6. 


five lengthy revisions, the final two types of data 


collection forms were developed, (Appendices 1, 2). 


Farm Data (form 1) — for collecting data from the farm 


environment, which included the hatchery and broiler 
houses. Such data were collected at least at two 
different times sequentially for the same brood of birds; 
when they were about four weeks of age and right before 


tney were processed at about six—seven weeks of age. 


Processing Plant Data (form 2) — for collecting 


corresponding processing and condemnation data from flocks 
for whicn farm data had been obtained. 

Confidential Information Form — on this form, only the 
names of cooperating firms and farms were kept under 
Strict confidentiality. The other forms bear only the 
codes of these farms and firms. 

With the aid of these predesigned forms, all visits 
made by research assistants of this project and data 
generated from hatcheries, broiler houses and processing 
plants (figure 8) were filled in. Tnis also meant that 
the logistics and scheduling of the field teams to gather 
such data at the correctly specified time was vital in 


obtaining valid and usable data. 


Protocol for data collection: A descriptive protocol of 


cm = 








the research was prepared and distributed to all members 
of the research team. The protocol was self explanatory, 
with detailed descriptions, directions and examples where 
necessary. Tasks to be performed were planned and 
Schedulea by the research associate, based on the 
implemented study design and routine managemental 
requirements. Major tasks to be performed during a given 
month were printed, and persons responsible for carrying 
out the tasks were assigned. Prescribed tasks were then 
carried out, following detailed weekiy task schedules. A 
Status report was filled out and followed by the research 
associate, as needed. In some cases, the tasks to be 
completed depended on the convenience of the study farm, 
and this was taken into consideration in planning and 
scheduling. Project staff meetings were held on Fridays 
to discuss and solve problems that arose. Changes of 
Schedules, travel plans, expense accounting and 
communication were directed by the research associate 
before reaching the principal investigators. 

SN 6 Data handling: 

Beal eoottware for handiing poultry health data base: Since the 
raw data collected from the various data generating 
centers (figure 10) had to be stored, retrieved and 
analyzed by using computers, a specific file handling 
software was developed for this project. The software 
used BASIC language and was designed in an interactive 


mode. It is simple and can be used to enter new data, and 
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Figure io DATA GENERATING CENTERS 


update or change existing ones in the data base, With a 
Short introduction to its use, persons with little 
training will be able to apply it. The data were stored 
on floppy diskettes in a model II Radio Shack 
Microcomputer. There are three main data files which 
-nclude: 
~ file for entering variable data for which no 
computation from tables were required. 
- file for handling data to perform analysis, 
using the biomedical data program (BMDP) 
Statistical package. 
- file for storing confidential information on 
the firms and farms involved in this study. 
To facilitate the analysis of data for descriptive and 
multivariate statistical methods using biomedical data 
program (BMDP package), FORTRAN files were created. Jn 
the early stages of the project since a large number of 
predictor variables (Xi) were involved, the variables were 
partitioned into: 
- hatch data 
- broiler data 
- biological and feed data, 
- transit and processing data 
Tnerefore, there were four sets of predictor variable 
groups, along with the respective condemnation rate data 
(Yi). Such groupings appeared to be satisfactory from the 
natural delineations that exist in the sequential events 


in poultry production, and the systems diagrams (figure 


(pe) 


Sf oo) 


Shes 


5-8) were useful in making such a decision. 

Vaizaation procedures: The data generating centers for 
this study and the need for coordination of these data 
inputs were indicated in figure 11. In order to manage 
the incoming data, proper planning, scheduling, evaluation 
of data for validity and completeness was required. The 
steps involved and method of handling this aspect is 
described in figure 11. Note that the descriptive 
protocol for management of the study corresponded closely 
with the mechanism laid out for data handling in this 
figure. The incoming data were cross checked manually by 
the research associate and the principal investigators at 
least three times before it was finally stored in the 
computer data base. 

It is significant to stress that careful data 
management is a crucial factor to the success of any 
project, and especially when a large volume of data is 
involved. Such as in this project, this point becomes 
even more critical. Thus, the data processing and 
management practices have been streamlined to ensure 
accuracy, validity and completeness throughout the life of 
the project. 

Data analysis: With a background on the essence of the 
problem that we seek solutions for, a brief review of the 
literature pertaining to the systems concept and the use 
of multivariate models to develop the predictive equations 


in this task is also timely. 
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5.8.1 Review of the literature on data analysis: The systems 


concept is a well tested problem solving scientific method 
which is widely accepted in engineering, economics, 
ecology, biological sciences, biometerology and human 
medicine (58-66). Although, it is relatively young in the 
area of veterinary medicine, it is now beginning to be 
applied to solving epidemiologic problems. Its future in 
veterinary medicine appears bright, and its application to 
examining the problem defined in this project could be one 
good example of such trends. The application of 
multivariate techniques to develop predictive or 
discriminative models is also widely used in various 
scientific and social disciplines where the background in 
the use of quantitative methods is well established. It 
is a standard statistical tool, the basis of which rests 
on sound mathematics (49, 67). In veterinary medicine, 
the application of multivariate techniques has been 
somewhat limited so far. In one such application, the 
best method for the control of African trypanosomiasis was 
investigated, and based on this result, rational disease 
control recommendations were obtained (68). Others have 
aiso used it to evaluate disease control methods (659, 70). 
In a recent study, discriminant analysis was utilized to 
classify a group of feed lot calves into poor or good 
doers (71). In the area of poultry health, the 
application of the systems concept and multivariate 


analytic models appears to be unavailable except for two 


eee 


Deosa.| 


Cases, both yet unpublished reports (72, 73). In those 
two reports, the relationship between various turkey pro- 
duction factors, condemnation and downgrading in turkeys 
were examined by using factor analysis. The data base for 
this study was obtained from the 1977 live bird production 
and processing data from one processing plant in 
California. Data from 48 flocks, representing 19 growers, 
were available for analysis. The results of that study 
indicated that "week of year placed", climate on day of 
processing and feed conversion may nave some predictive 
Value, Since the literature on the application of systems 
analysis and muitivariate models in veterinary medicine 
and more specifically in poultry health is rather scanty, 
it is believed that the results of this project would add 
a little more to our present knowledge. 

As mentioned earlier, relevent variables that affect 
the health of a population of birds at given sites and/or 
times, served as potential predictors for the prevalence 
of poultry diseases on farms, and the condemnation rate at 
slaughter plants. Data gathered for such a set of 
explanatory or predictor variables (Xi) were examined for 
significant effects on the dependent variable (Y) viz. 
condemnation rates at slaughter plants, by using the 


biomedical data program (BMDP) statistical package (57). 


Data Analysis Methods: 


Descriptive Statistics: Detailed data descriptive 


Statistics were conducted by using PMDPeD. This program 
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enables counting and listing of the distinct values for 


each variable in the analysis. It computes univariate 
statistics including the mean, median, standard deviation, 
skewness and kurtosis. It then plots the positions of 
several estimates of location on a line. From the 
histogram and other descriptive statistics, the type of 
distribution and variabilities present in the data were 
assessed. P2D was followed by a scatter (bivariate) plot, 
whereby each of the predictive varaibles (Xi) were plotted 
against the dependent variable i.e. condemnation rate 
(Y,), by using the BMDP6D program. This program computes 
and prints equations of simple linear regression and 
indicates the intersections of the regression lines with 
the frame of the plot. From such a seatter diagram, it 
was possible to obtain information on linearity or on the 
type of transformation which may be needed to lead to 
linearity. 

5.8.2.2 Comparative (analysis of variance-BMDP2V): This program 
performs an analysis of variance (or covariance) for a 
Wide variety of fixed effects models. In this analysis, 
the levels for a factor were specified by grouping it in a 
manner that would allow each factor to be evaluated, 
reither individually (one way analysis of variance) or 
cross classified with another factor (two way analysis of 
variance). The dependent variable with which each of the 
factors were evaluated was Y,, the condemnation rate at 


2 


processing. 


As an example, flock size X3 could be grouped into 
three levels: flock size of 20,000 or less, 20,000 to 
40,000 and over 40,000. The other variable to be 
considered could be Xi09° the disease rate at necropsy. 


Three levels may be established for X as 0.2 or less, 
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0.2 to 1.0 and greater than 1.0. An analysis of variance 
using BMDP2V could be conducted to test the hypothesis: 
a) test of equality of row means: 
b) test of equality of column means: 
ec) test of no interaction: 

In this study, for example, the hypothesis tested 
were: 

There is no difference in the response variable 
for the three categories of flock size Ox) which 
were laid out across rows. There is no difference 
in the response variable for the three categories 
of aisease rate (x 199) which were across columns: and 


that there is no interaction between the row and 


column factors. 


One should note that in creating 3 or 4 levels of Xa 
(flock size) to set it in an ANOVA format, what one does 
is really group data in summary form rather than examine 
the true spread of that variable by looking at its actual 
Values in a more continuous setting in the form of 
groupings or levels. More importantly, the basic 


objective was to devise a predictive model; not to assess 


Deueces 


differences between variables. 

Multivariate linear model layout: As indicated in the 
research proposal, the major thrust of the analysis was 
based on linear multivariate models. It was proposed that: 

Vee! HOO) 

and the linear multivariate model: 

Yar= eGo °61X ii cone ta BKAK: hues 

Where, 

Y = condemnation rate at the processsing plant, 

X, -- Xe = determinants of condemnation rate, variable 
data collected from farms. 

Bo -- Bk = coefficients to be estimated, which 
indicate the relative strength or 
importance of each variable (X) in 
predicting Y. 

ek= random error term — unexplained variability. 

A diagrammatic view of the linear causal model which 
is an abridged and simplified form of figure 8, but one 
which is more straight forward to see, is given in figure 
12. Note that figure 8 was a dynamic representation of 
the problem, while figure 12 was the summarized static 
(still picture) form of the same problem. The array of 
data layout for the flock testing project compatible with 
the diagrammatic view of figure 12 is presented in 


figure 13. 
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The specific linear statistical models employed to 
develop the predictive model for condemnation rate were: 
a. Multiple linear regression — the objective was to 
establish a predictor equation for condemnation 
rate from a set of predictor variables, The 
assumptions in using multiple regression models, 
are standard as presented in regression texts 
(49, 50, 67) 
Each of the assumptions were evaluated by performing 
residual analysis to see if any violations were made. 
Such a discussion would be found under the results and 
discussion section. For performing the regression 
analysis, the BMDPIR and RMDP2R programs were utilized. 
The BMHDP2R was especially important since it is a stepwise 
procedure which facilitated the selection of the best 
subset of predictors for condemnation rates. This program 
(P2R) computes and estimetes parameters of a multiple 
linear regression equation in a stepwise manner, 
Variables are allowed to enter into an equation (forward 
stepping) or are removed (backward stepping) from the 
equation, one at a time, according to any of four possible 
criteria to be selected. In this analysis, the F method 
was used, whereby the variable with the smallest 
F-to-remove was removed if its F-vaiue was less than the 
Preassiened value. If no’ Variable meets this criterion, 


the variable with the largest F-to-enter was entered if 
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the F-to-enter exceeded the F-to-enter limit (57). 


Assuming that some of the predictor variables may require 
transformations e.g. logarithmic, such options were 
exercised as needed. P2R computes univariate statistics, 
covariance and correlation matrices, analysis of variance 
tables at each step, F-to-remove and F-—to-enter values, 
and other information. A summary of variables entered, 
the multiple correlation coefficient (R) and the 
coefficient of determination (Ro) were also provided along 
with residuals. A zero intercept equation useful for 
establishing the descriptive and comparative value of a 
variable, was provided. Since forward stepping was 
followed by backward stepwise regression, the potential 
predictive value of the variables were not lost. 

Aithough the decision making criteria for variable 
elimination or selection from the data base will be 
discussed later, it is important to point out that the 
results from multiple regression would be evaluated by 
using F values, R and Ro levels to determine whether a 
selected variable(s) has biologically explainable 
correlation with poultry condemnation, 

Ro measures the proportionate reduction of total 
Variation of the dependent variable (Y,) associated with 
the use of a set of independent (Xi) variables. 
Generally, the higher the value of nae the better the 
usefulness of the model. Correspondingly, R measures the 


Strength of correlation between the dependent and 





independent variables. Both of these values are useful 
indicators of the predictive power of the model. 

b. Discriminant analysis is the other linear multi- 
variate model used in this study. Tne objective 
in this case was to obtain a subset of 
discriminator variables which would enable one to 
classify a flock of birds from a given farm 
and/or house into high “or low condemnation 
groups or more subgroups, based on predetermined 
criteria for such a decision. 

To accomplish this, BMDP7M, a stepwise discriminant 
analysis procedure was employed. In this method, 
variables used in computing the linear classification 
functions were chosen in a stepwise manner with both 
forward and backward stepping options. At each step the 
variable that adds the most to the separation of the 
groups was entered into, or the variable that adds the 
least was removed from the discriminant function. A Jack 
Knife-validation procedure option was used to reduce bias 
in group classifications. In all statistical analysis to 
be performed, the level of significance was evaluated at 
either *the0.05,0or,0.01.levels. 

It should be reiterated that the main objective of 
this research was not to simply establish differences or 
associations but to develop a predictive model. 
Tnerefore, the two linear multivariate models presented 


above were selected as the most appropriate ones. 
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5.8.5 Procedure for selecting the final predictor model(s): 


After the total number of predictor variables have been 
reduced, the final aspects of data analysis to establish 
the predictive or discriminative models, was conducted in 
the following manner. Using an individual house on a farm 
as the experimental unit, cases for which a complete set 
of data were available for, were analyzed first. Secondly, 
using a farm as the experimental unit, comparable data 
analysis were performed. A decision as to whether to use 
the house or the farm as the experimental unit of choice 
was made based on detailed analysis of the data. Finally, 
to perform the validation of the model, the total number 
of cases were divided into two groups. The first set of 
cases formed the study data set. Based on these 
observations, the predictive/discriminative models were 
developed. The second group of cases served as the test 
data set. Once the multivariate models have been 
establisned, the validity or performance of the models 
were assessed using the cases from the test data set. In 
developing the predictive/discriminant models, BMDP2R and 
BMDP7M programs were used. For the discriminant model, 
various cut off points were used and the options discussed 
accordingly, including the probabilities and the 
consequences of misclassification. Results from the 
detailed descriptive statistics provided via PMDP2D and 
BHDP5D, and other statistical tests, specifically, 


analysis of variance were utilized where needed to reveal 





additional useful information. 


5.8.6 Assessment of adequacy of model and validation: Once the 


predictive/discriminant models have been established, 
using study data, the validity of the models were 
evaluated by entering these data into the predictor 
equation and predicting the condemnation rate before 
inspection, and then comparing the prediction with the 
actual condemnation rate. A very important step at this 
stage was the assessment of the adequacy of the model by 
performing detailed residual analysis. A variety of 
computer generated residual plots were employed to 
identify the nature of lack of fit, outliers, deviation 
from normality and other deficiencies. This aspect of the 
analysis addressed the question of the basic assumptions 
involved in using linear models. Such assumption as 
linearity, homogeneity of variance and normality of 
distribution were evalauted. Where inadequacies existed, 
the variables of interest were transformed, for example to 
respective exponential or quadratic forms to improve the 
adequacy of fit. Correspondingly, confidence intervals 
were computed for the respective regression coefficients 
in the predictive model using Bonferroni's method (49), 
These techniques are all related to evaluate the validity 
of the predictive model and should provide useful 


information. 


5.8.7 Weighted least squares: A weighted least squares 


procedure was used in the analysis of the data via BMDPOR. 
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Weighted least squares estimation is appropriate when the 
error variance is not homogeneous but varies from case to 
case. In the data base used in this analysis, the 
estimated standard deviations of the regression 
coefficients and the coefficient of variation of some 
variables indicated that weighted least squares could be 
useful to improve the predictor equation. 

A two stage least squares procedure was used here as 
described below. 
Weighted least squares: in BMDP2R weighted least squares 
procedure, the goal is for the estimates of the 


coefficients to minimize: 


2 
al 2 


Where Wi is the weight for case i. The problem though is 


EWE CY yan) 
i 


in obtaining a value for Wi for each of the cases. 
However, we know that: 

3 Gel) atk QO hee ORD) 
and an appropriate weight estimate could be obtained by: 

Wiles j 
LES CYS) ere eC 

But, 

RWI) e=36o. + bie 
indicating that the above weighting procedure involved 
unknown parameters. One way to solve this problem was by 
using a two stage least squares procedure, i.e. in stage 
1, using BMDP2R, stepwise forward and backward stepping, a 


regression model was obtained using a specifie predictor 
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equation. 


obtained 


In stage 2, estimates of the weights were 


from the results in stage 1 above where: 


Wi = j 
Y2iIgG-y2i4) 


Tnese estimated weights were then used to obtain the 


weighted 


least squares regression model. To enable this, 


a small Fortran subroutine was included to calcuate Wi as 


given above and in the BMDP2R program, the variable which 


contained case weights was indicated appropriately. The 


adequacy 


residual 


of the model via weighted least squares and its 


analysis were performed as described earlier. 


Hardware and Software for Data Analysis (BMDP Package): 


The anaiysis of the data collected was performed using 


BADP statistical software on a VAX 11/750 computer at 


Tuskegee 


Institute. The sequence of statistical analysis 


consisted of: 


BMDP2D — 


BMDP6D — 


BADP2V — 


BMDP1R — 


BMDP2R — 


descriptive statistics for each variable; 

scatter plots of each predictor variable vs. the 
dependent variable; 

An analysis of variance program for specified 
factors; 

multiple regression program whereby all variables 
were entered for a specified dependent variable; 
a stepwise forward selection regression procedure 
followed by backward stepping. This enabled the 


evaluation of each predictor variable so as to 


5.9 
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eventualiy select the best subset of predictor 


variables; and 

BMDP7M - a stepwise discriminant analysis program to 
develop the classification function into high/low 
condemnation groups. 

Decision making criteria for deleting variables from data 


I 


base: From an understanding of the nature of the problem 





and based on a detailed causal analysis of variables which 
may influence condemnation rate of poultry (figure 8), 286 
potential predictors (X's) were identified. This was a 
lengthiy list of variables and systematic screening was 
needed to reduce the number of variables to a manageabie 
level. The ultimate objective was to select only a very 
few, possibly less than 10 subset of predictor variables 
which could then be used for predictive or discriminative 
purposes. The reason for shortening the list of predictor 
Variables was mainly because a multivariate model with a 
limited number of independent variables was easier to 
analyze and understand. But one has to be careful not to 
compromise predictive power. The search for the "best" 
set of independent variables is complex. ‘The decision 
making criteria devised for this study were; 
Phase I preliminary data’ sereening variables deleted if: 
a. data on specific variables were missing and were 
unavailable, 
b. variable was redundant 


c. variable was not of fundamental importance to the 








problem, 

computations of specific rates which enlarged the 
raw data for a variable e.g. disease specific/age 
specific mortality/morbidity rates, In some cases 
such specific rates were unavailable e.g, morbidity 


rates, 


5.9.2 Phase II second stage screening: 


ae 


Based on descriptive statistics, variables with large 
variability especially those with large coefficient of 
variation were considered for deletion. 

From scatter plot diagrams, variables were examined 
visually as well as by their simple linear regression 
results to see if transformations, such as logarithmic, 
quadratic forms or interactions could improve their 
linear form. If not, these variables, were considered 
as candidates for deletion. 

From a printout of correlation matrices of each group 
of variables, highly intercorrelated variables (high 
multicollinearity) were considered for reduction, 
After fitting a regression model containing an entire 
set of potential independent variables thich may 
predict condemnations, those variabies wich had small 
absolute value of the t statistic were also considered 


candidates for screening. 


Meth ose stage wrdaitnomen anyvof the above criteria did not 


automatically lead to deletion of the variables, those 


that failed any two or more of the above criteria were 


2M) 
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dropped from further consideration. It was believed that 


at the end of the second stage of screening, the final 
potential predictor variables could be narrowed down to 
between 30 to 50. 

Phase III final stage of selection for the "best" set of 
independent variables for a predictive and/or 
discriminative model: The final screening of variables to 
arrive at the "best" set of predictor variables was 
basically accomplished by a computerized stepwise 
selection process which utilized forward and backward 
Stepping techniques. The stepwise regression search 
method was widely used in many fields. In this procedure, 
the basic criterion for adding or deleting a predictor 


variable equivalently involved: 


reduction in sum of squares, 


partial correlation coefficient, or 
hee Siete st 1c 

Quite simply, at each step, the F statistic was 
computed by the computer program, and the predictor 
Variable with highest F value was the candidate for 
addition. Conversely, in backward selection, the variable 
with the smallest F vaiue was identified and if it was 
less than some predetermined limit, the predictor variable 
Was dropped. More details on the involment of such a 
selection process, specifically, when using BMDP2R. and 
BMDP7M are available elsewhere (57). Once a set of 


predictor variables were identified, the final model was 
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examined further by performing a residual analysis to 
evaluate the adequacy of the model. A variety of residual 
plots were employed to identify the nature of lack of fit, 
outliers, deviation from normality and other deficiencies. 
Tne model building process did not stop at the conclusion 
of the selection of the "best" subset of predictor 
Variables discussed above; until the final predictor 
equation was validated. This was accomplished by testing 
the predictor power of the model on a new set of data. 
Strategy for data analysis: The first decision to be made 
in the data analysis task was to determine whether to use 
a house or the farm as the experimental unit. Data have 
been collected using both of these as study units and the 
analysis using two sets of data would have been 
cumbersome, expensive and time consuming. To settle this 
question, the following plan was devised. Initially to 
facilitate analysis, the data set were subdivided into 4 
partitions or submatrices. This was necessitated since the 
number of variables in use at one time was large and 
beyond the limits of the BMDP package. However, during 
the last part of the project the BMDP package had already 
been expanded to handle up to 200 variables at the same 
time. Therefore, the partitioning of the matrices was not 
necessary in the final analysis of the data. On the other 
hand, for the decision of whether to use the house or farm 
as the study unit, the partitioned matrix approach was 


utilized since such an analysis was performed a few months 
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earlier. The partitioned matrices were labeled: 
a. Hatchery data-which contained data of hatchery 
related variables 
b. Broiler data-which contained variabies related 
to broiler house environment 
ec. Feed and biologicals data-which contained data 
pertaining to the use of feed and biologicais in 
poultry raising 
d. ‘Transit and processing data-which contained 
variables dealing with the final phase of 
poultry production and processing 
Such partitions appeared satisfactory from the 
delineations that are inherent in the sequenced events in 
poultry proauction. The systems diagram (Figure 8) was 
used for making such categorizations. ®y using the 
dependent variable (Y,) with each of the above partitioned 
matrices, the. previously indicated sequence of BYNP 
programs were conducted. The goal was to reduce the 
number of variables from each of the respective matrices 
so as to end up with one final data base of about 30-50 
predictor variables, thus compressing the 4 partitioned 
matrices into one data matrix for final analysis. After 
the decision of which experimental unit to use in the 
analysis was made, the final aspects of data analysis to 
develop the required models was initiated. The final data 
analysis were conductea after: 


a. Data for the last brood cycle (brood 5) was 
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entered on flock testing computer data base and 


all other information was updated, corrected and 
verified. 

b. Tne respective disease specific condemnation rates 
in addition to the condemnation rate due to all 
diseases (Y,) were entered into the data base. 
These were: Yg (leucosis), ‘40 (septicemia), Yay 
(airsacculitis), Y42 (synovitis), Y43 (tumors), 
Yay (others). 

c. The respective Ye-4 values (condemnation rate at 
one previous brood for each study farm) was 
entered along with the disease specific 
condemnation rates into the data base as one 
additional potential predictor variable. 


Gea Data Orne. the disease rate variables based on 
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the necropsy profile performed at Tuskegee 
Institute, was entered along with the 
corresponding respective necropsy based disease 
specific rates for leucosis (Xa49)s septicemia 
(Xi 44) airsacculitis (Kayo) synovitis (X4 43) 
tumors (Xo 44) and others (X445)- 

Once the decision of which study unit to use was made, all 

the predictor variables in the poultry data base were 

nandlea simultaneously in one sweep without resorting to 

partitioned matrices. This way, the influence of each 


Variable in condemnation rate evaluated simultaneoulsy 


along with all other potential predictor variables. 
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Therefore, a sequence of BMDP of analysis were performed 
with the specific objectives of: 
a. Developing a predictor/discriminator model for 


Y5 by using X as one of the many independent 


109 
Variables. Each of the disease specific condemna— 


tion rates (Y, = Yay) were handled in the same 


9 
manner as Y5 but using their respective necropsy 
based disease rate variables (X)497% 445) 

b. Developing a predictor model for Xr 09 itself as 
the dependent variable using the remaining 
independent variables. Each of the necropsy 
based disease specific rates were then evaluated 
individually by replacing X109 as the dependent 
variable. 

c. By replacing the value of Xi09 in the predictor/ 
discriminator model for Y, using the result of 


the prediction of X from section b, the 
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predicted value for X109 was used in the model 
thus bypassing the need to perform necropsies 
before the birds were processed. However, one 
has to accept the level of error present when 


using the predictor model for X,.,. as an input to 


another predictor equation. 


ape) Limitations faced in the project 
2.10.1 Research design related: A problem inherent in the 


research design itself was the lack of a random sampling 
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of study firms and farms to serve as the experimental 
unit. Because of this drawback, the inferences to be 
made from the results will have to be guarded. It is 
important to note that although such a sampling 
procedure was desirable from the statistical standpoint, 
it was neither practical nor feasible. From a total of 
11 major poultry firms in Alabama, nine of them were 
requested to cooperate in the study and only three 
accepted to do so initially. Eventually two others 
joined the project. The firms studied were satisfactory 
in terms of the representations of the poultry firms in 
Alabama. On the other hand, the random samplings of 
growers would have been desirable; however, the question 
of cooperation and logistics were overriding factors in 


the decisions made in this study. 


Field data collection related problems: Problems 


encountered in dealing with growers as well as those 
wnich arose in the processing plants are presented 


below, in summary form: 


a. Lack (or loss) or cooperation of study units. 
Some growers restricted the visits made by the 
research team, and data gathering was hampered. 

b. Unavailability of farm data-records on such 
variables as morbidity, medications, feed 
utilization and microclimatic data 


(temperature and humidity) were not available. 
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When requested to keep some of these records, 
most of the growers were unwilling to do so. 
ec. Some of the study farms ceased to raise 
chickens, and data on these farms could not be 
followed up during the whole year. 
d. The logistics of gathering data from 
hatcheries, farms and processing plants was 
beyond the original contractual agreement of 2 
farms per firm. 
Originally, quite often, 2 or 3 teams were sent to the 
field at the same time performing microbial samplings of 
hatcheries and farms, collecting weekly feed utilization 
data, visiting each brood of broilers on different farms 
to collect blood samples and existing data, and 
collecting blood and necropsy samples from processing 
plants. This usually involved long hours of driving at 
night and day. Because of the exhaustive schedule as 
well as shortage of manpower, the number of study farms 
were later reduced to 3 or 4 per firm while still 
Maintaining representation across the three strata of 
production performance rankings. 
e. Weekly feed utilization data remained elusive 
Since such data were not kept by the poultry 
growers, An estimation procedure was 
developed toa accomplish this phase of the 
Study. Since the procedure involved rough 


estimations, the accuracy and usefulness of 
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such data was limited, it did however, provide 

some preliminary indication of patterns of 

feed utilization among farms. Due to the 
demand on time and manpower, the data 
collection of weekly feed utilization from 
study farms was conducted in the following 
manner: 

- Data were collected from north Alabama study 
farms only if other tasks such as 4 or 7 week 
data collections or microsampling were 
performed, thus limiting the amount of 
travel to a manageable level. This also 
meant that feed utilization data was un- 
avaiiable on regularly scheduled weekly 
visits. 

- Data from study farms in south Alabama 
were collected on a weekly basis on week- 
ends since these farms were relatively 
close to Tuskegee Institute. 

The possibility of developing a flow meter 
type of device which could automatically 
record the volume of feed being utilized 
on chart paper, on a daily basis was 
explored, Although the idea appeared 
feasible, it was not pursued further due 
to lack of time for development of such a 


device, This is an engineering problem, 
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and one has be to intrigued why the industry 
has not already developed such a feed 
utilization monitoring device since feed cast 
is one of the major expenses that the grower 
has to contend with. ur discussion with a 
few growers on this subject and on the 
possibility of designing such a device 
confirmed the need and the demand for such an 
instrument. 
Tne collection of condemnation data by house 
was another problem area. One company 
declined to assist: us«in thercollection of 
such data. Occassionally, even in other 
processing plants who have volunteered to keep 
track of slaughter birds on a house by house 
basis, flocks from two houses were mixed 
together during shipping of birds from the 
farms to slaughter plants. Although every 
effort was made by the research team to 
maintain the identity of lots on a house 
basis, the problem appeared to be ever present 
and because of that, some condemnation data on 
house basis were unobtainable. 
Unavailability of trim data due to diseases: 
condemnation due to diseases involved either 
whole carcass or parts trimmed. Records 


furnished by IIC of FSIS on whole carcass 








condemnation were accurate. However, data for 
trimmings due to diseases was unavailable. 

The processing plants kept all trimmed parts, 
regardless of causes for condemnation, in one 
container, i.e. trimmings due to diseases as 
well as those due to plant errors were all 
placed in one container and then weighed. 
Diseased parts trimmed and condemned were for 
such categories as necrotic dermatitis, patchy 


skin disease, breast blisters, and slip 


tendons. Conversely, trimmed parts due to 
plant errors include such items as mangled 
birds, bruises, broken hones or skin, 
hemmorrhage in legs, wings, thighs or backs, 
etc. All these condemned parts were placed 
together, since FSIS does not require 
processing plants to acquire and maintain 
records on trimmings. 

Yet, a considerable amount of time and energy of the 
inspector on the line were actually devoted to identifying 
and indicating what parts were to be trimmed and how much. 
Trimmings due to diseases appear to contribute 
substantially to the condemnation rate of a flock. The 
keeping of such data should be standard practice by FSIS, 
Since both in terms of disease surveillance, and 
inspection manpower requirements, the taks of trimming 


would have to be considered emphatically. Additionally, 
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the expertise required in identifying trimmed parts due to 
Such diseases as tendonitis, abdominal flap sore, patchy 
skin disease, hock deformation, breast blisters trims, 
etc. may not be standard from plant to plant. As the 
disease condemnation of whole carcass decreases, 
condemnation due to parts trimmed may become significant 
and yet as the system operates currently, such data are 
unavailable. Therefore, relevant comments on such trends 


could not be made. 


Laboratory related problems: 


a. Necropsy profile, (X109)-described in great detail 
earlier was an attempt to obtain an estimate of disease 
rate on a study farm from 200 randomly selected birds 
at processing plants for necropsy examination at 
Tuskegee. Since the purchase of birds was a 
rather costly one, data for this variable was 
somewhat limited. 

b. Microbial sampling: The need for an expanded micro— 
bial profiling work has been presented in our systems 
diagrams (Figure 8). The necessity for performing 
microbial population estimates by colony counting as 
well as the identification of the types of organisms 
in the hatchery and poultry house ecosystem has been 
described. Tne identification of microbial flora 
using semiautomated equipment to handle the large 
volume of samples in this project was essential. 


Unfortunately semiautomated microbial identification 
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equipment and supplies were not included in the 
original budget. Because of the strong conviction of 
the investigators on the need for this aspect of the 
study, a limited number of samples were identified by 
using manual procedures during brood cycles 4 and 5. 
The limited nature of this part of the study provided 
only descriptive information. The number of cases 
generated were not sufficient enough to be evaluated 
along with the other predictor variables. One cannot 
overemphasize the potential importance of microbial 
profile data in affecting condemnations (especially 
the types of organism which predominate the poultry 
ecosystem) since these are central in affecting 


the biodynamics of poultry health. 


5.10.4 Data analysis related constraint: One of the crucial 


techniques for data analysis calls for the use of a 
discriminant analysis model. The objective was to devise 
a discriminating mathematical equation to separate 
condemnation rates into at least two groups. This was 
accomplished only if prior decisions are specified as to 
how to separate condemnation rates into high or low; high, 
medium or low categories. This meant specific numerical 
values or cut off points had to be established prior to 
analysis of the data. In the progress report of August 
1981, FSIS was requested to provide such a differential 
eriteria so that the discriminant classification functions 


could be developed. However, FSIS declined to provide 
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such guidelines since they did not have numbers on what 
constitutes "high" or "low" condemnation rates. It should 
be noted that even for the predictive models being 
developed, prior demarcation lines have to be established 
for condemnation rates in order to provide guidelines as 
to when to implement possible scaled down inspection 
procedures, as opposed to continuing the current practice 
of inspection. Therefore, in this study arbitrary cut off 
points based on the existing data base were utilized in 
the analysis. Since a practical decision making criteria 
was unknown, therefore, time was spent trying various 


options which could have been avoided otherwise. 


6. RESULTS AND DISCUSSION 


onl 


The data base: The research design whereby the 


determinants of poultry population health were 
systematically identified and scrutinized via causal 
diagrams and predictive multivariate epidemiologic models 
have been described earlier. It was also noted that 
poultry (broiler) production involves a sequentially 
determined set of events which include breeder flock 
sources, the hatching, brooding and growing ecosystem and 
Variables which result from poultry processing efforts. 
All variables in these various sites were identified, and 
appropriate quantitative data were collected, forming the 
data base. 

A tabular summary of the study subjects (firms, 
farms, and the individuals houses of cooperating farms) 
for the respective brood cycles is provided in table 5. 
The complete set of data on both the independent and 
dependent variables is in Appendix 4; the complete list of 
the variables for which adequate data was collected is in 
appendix 9. The data base contained a maximum of 154 
cases for 141 variables. Of these only variables Yos 


Y -Y 4 and X 


= t 
ot 1H? “109 11077 115° served as the dependen 


variables, while those with subscript Xr (113 variable) 

formed the independent variables. The other Y subscripted 
gud Fa = 4 i: 

variables Yas Y3 YQ. and Yo5 Ya were included for 


comparative purposes only. Although there were a maximum 


id 
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TABLE 5 Study Units For Flock Testing Feasibility Study, Tuskegee Institute (1981/ 


ck BROODS TOTAL # OF 
STUDY UNITS _ Sl et 5 OBSERVATIONS 

No. of Firms 3 iS paps 

No. of Farms 2A 28 138 

No. of Houses 455 (4). 1585026) 16/027) 689012) 7 50 RC on 299 

No. of Houses with Condem- 

nation data 6 16 87 
No. of Houses with Necropsy 

and Condemnation data 6 0 47 


No. of Farms With Condem- 
nation data 2A: 





28 138 








of 154 cases, all cases did not contain the complete set 
of predictor variabe data corresponding with each 
agependent variable. Therefore, the number of cases with a 
complete set cof data varied, depending on which 
multivariate variable were used. 
Preexisting data ~- a major facet of this study was the 
collection of data from existing records kept by a broiler 
farm, processing plant or hatchery. Such information was 
acquired at specific times and for specified variables 
(Appendix 1, 2). Selected examples of those which 
comprise part of the data base are given in tables 6 - 8. 
The first table provides a segment of the poultry 
population profile and significant paramets in the flock 
testing study. In table 7a, specific preexisting data for 
condemnation rates of chickens from the five study firms, 
which spanned a year of study, are presented. Table 7b 
presents related information on a disease specific basis; 
the specific disease being leucosis, septicemia, air-— 
sacculitis, synovitis, tumors, and others, as provided by 


Polo inspectors Operaingy im “the five, study ‘firms. 


Field generated data ~ the data collected in this manner 


formed part of the flock testing study data file, a 
computer printout of which’is tin appendix 3.¢ Of the field 
generated data, the one with special significance is feed 
utilization. It should be noted that a good portion of 
feed utilization data at 4 weeks and 7 weeks of age was 


available from the records of the firm or poultry farm. 
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TABLE 6: Population Profile ~ Flock Testing Feasibility Study, 
Tuskegee Institute 1981-82) 


STUDY | UNIT FLOCK SIZE AT | MORTALITY RATE (%) _| CONDEMNATION RATE (%) 
Wk. #1 : PROCESSING 





TABLE 6: Population Profile - Flock Testing Feasiblity Study 
Tuskegee Institute (1981-82) 


STUDY UNIT FLOCK SIZE AT MORTALITY RATE (%) | CONDEMNATION RATE (%) 
FARM Wk. #1 ease A Wie x PROCESSING 


1.417 
0.450 
0,747 
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0,288 
OFZ 3) 
0.034 
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0,662 
dee 27 
26,103 
207 
0.401 
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0.816 


0.686 
The Si Re) 
0.940 
isle: 
Uo 
0.380 


pe03o 
0.472 
0.330 
0.412 
0.414 
Oe321 


e720 
0.464 
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TABLE 6: 
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STUDY UNIT _ FLOCK SIZE AT 


FIRM 
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Population Profile - Flock Testing Feasibility Study, 
Tuskegee Institute (1981-28) 


#1 


MORTALITY RATE (2%) 





CONDMENATION RATE (%) 


PROCESSING 


0.344 
Ont 2 
Ones) 
02226 


5.800 
d is Pa 
dag 65 
L002 
0.865 
0.982 
0.947 
0.314 
Aig die ys 


0.296 
0.754 
0.790 


02232 
0.203 
0.205 
0.833 
Wee eI 


0.938 
Gras 
Coats 
0.494 
O.. 529 
0.742 
0.700 


0.367 
0.243 


0.248 
0.114 





TABLE 6; Population Profile - Flock Testing Feasibility Study 
Tuskegee Institute (1981-82) 






STUDY UNIT FLOCK SIZE AT MORTALITY RATE (%) | CONDEMNATION RATE (%) 


ate eee ewe 4 We. [7 Wk. LS ee PROCESSINGD ss = 

48310 1.822 2.593 
48255 1.980 2.829 
43790 1.644 M 

53155 1.813 2.876 
64411 2.595 3.657 
24774 Sw) M 

11352 1.378 2.463 
17750 | 2.979 5.847 
30260 1 pi 4.102 
36276 | 3.692 4.928 
27046 | M M 

15504 2.167 2.816 
16133 1.911 2.843 
49193 iL opis 2.379 
68944 jargeyAl 2.999 
15320 3.092 3.991 
20774 2.542 3.960 
32262 1.875 2.442 
26913 1.826 ae we 
22266 ST 4.109 
42324 2.407 3.015 
31862 1.503 aloes 
25222 2.027 2.667 
47463 3.435 4.729 
47982 1.937 2.860 
42910 1.941 2.622 
53285 1.911 2.929 
29183 2.860 | 3.229 
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Table 7a: 
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Condemnation Rate Due to Diseases for Study Farms 


Of 


oO oO @>2 © OC © O0Oro0o09cjn0eco 


OO OF NK O 


(e} (@) (2) (= 


o747 


APav | 


Zoe 


662 


pe Joy 


444 


5092 
e30l 


Brood 


eooooocko 
> 
N 
wn 


ocOoOrrF OW (2) (2) ore 
(on) me Ww 
oO Ww > 
ie) ~N > 


- Ore © 
Ww 
— 
No) 


(eee) (eo) (e)(2) 
No 
io) 
Go 


oO 


(=) (=) 


@) (=) (e) (e) te) (oJ (=) 


oOrore 


eo) (oe) @) (&) (2) SS =) 


Orr Oo 


(ey le) e) (ey ley iy’) 


118 


119 


Table 7b: 


Condemnation Rates and Disease Specific Rates by Study Firms for 5 Brood Cycles 
(1981/82) 

















Firm 1 Brood 1 Le O23 0: 
Firm 2 Brood | He. 0.005 CF 
Firm 3 Brood 1] 0. 0.510 OF 
Firm 4 Brood 1 0: 0. Po O- 
Firm 5 Brood 1 0. O20 i). 0. 
Firm 1 Brood 2 0. 0.118 0.2 
Firm 2 Brood 2 0. 0.01 cs 
Firm 3 Brood 2 ie 0.704 0. 
Firm 4 Brood 2 0: 0.134 0. 
Firm 5 Brood 2 OF 0.091 GO; 
Firm 1 Brood 3 0. 0.060 Ge 
Firm 2 Brood 3 0. 0.008 O7 
Firm 3 Brood 3 1s Laan OF 
Firm 4 Brood 3 QO. 0.087 0. 
‘Firm 5 Brood 3 O¢ 0.094 0. 
Firm 1 Brood 4 0. 0.109 0. 
‘Firm 2 Brood 4 Oe 0.014 0. 
Firm 3 Brood 4 0. U7225 QO. 
Firm 4 Brood 4 Ox 0.05 0. 
Firm 1 Brood 5 0; 0.148 Oy 
Firm 2 Brood 5 QO. 0.022 0 
3 5 0. 0. 0. 
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However, the one of interest from the standpoint of 


population health characterization, was weekly feed 

utliization rate. As described under methodology, this 

task was rather difficult and expensive and only a limited 
amount of information was available for it. An example of 
one such data is presented in figure 14. Since the number 
of study farms for which such data could be collected were 
few, the feed utilization on a weekly basis has not been 
in the analysis. Other variables for which field 
measurements or evaluations were conducted included weight 
of birds at 4 and 7 weeks, (Xp) and proportion of 
undersized birds (X7). In such cases, 50-100 dirds were 

Selected arbitrarily and the necessary procedures 

performed to generate the data of interest. 

Laboratory generated data — three specific types of 

laboratory related profiles are considered here. 

a. Microbial profiling -— two aspects were involved. 
Firstly, the quantitative extimation of the microbial 
population size in hatcheries and in the poultry 
houses as described earlier. Secondly, the isolation 
and identification of the types of microorganisms 
which predominate in the broiler house and hatchery 
environment. The results of this work are provided 
in Appendix 5. The limited studies on 
microbial profile in hatcheries and broiler farms did 
not offer clear indicators of whether bacterial and 


fungal populations in hatcheries and chicken houses 





























LBS. OF FEED CONSUMED PER BIRD PER WEEK 


Figure 14 WEEKLY FEED UTILIZATION OF BROILER BIRDS 


IN SELECTED FARMS (example 3) 


Flock Testing Feasibility Study, Tuskegee Institute (1981/82) 
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would directly relate to condemnation. There were 
few bacterial and fungal colonies isolated from the 
hatcheries, however, the bacterial and fungal 
populations in broiler houses were much larger. 
Non-pathogenic or low pathogenic bacteria do not 
pose serious health problems for poultry flocks. 
However, the pathogenic ones cause diseases of 
Significance to the poultry industry. Therefore, 
the identification and characterization of the 
bacteria isolated from these study farms were 
required. Colonies from twenty-six cases were 
Studied, the data of this study is summarized 

in Appendix 5. Most bacteria isolated from 

the hatcheries and broiler houses were non- 
pathogenic, however, several potential pathogens 
(bacteria and fungi) were also isolated. These 
were Staphylococcus, Moraxella spp. 

Pasteurella, Pseudomonas, Proteus, Listeria, 
Corynebacterium, E. coli, Kiebsiella, 
Streptococcus, Aspergillus, Rhyzopus, Mucor, 

and Penicillium. One of the potential causes 

of airsacculitis, Aspergillus spp. were isolated from 
11 cases out of 26 studied. The result of this 
Study indicated that there numerous potentially 
pathogenic organisms in the hatcheries and 

broiler houses examined during 1981-82. These 


organisms may cause serious health problems 


under favorable environment, Therefore to 


monitor bacterial and fungal population in 














| poultry houses and hatcheries may be necessary 
for a well managed poultry enterprize. Due 

to the limited nature of the data studied in 

| comparison to the rest of the multivariate 

| data collected during the life of the project, 
the data obtained from the microbial profile 

were not analyzed with the rest of the data. 

b. Serology prota ter= therresui tot sunis work is 

| provided in Appendix 6, for the test sera which 
were collected from birds at 4 weeks of age. 

The titers for Mycoplasma synoviae were negative 
except for 5 study farms. These had a titer of 

| 1.:10. In the case of Mycoplasma gallisepticum, 
only one study farm had a titer of 1:10, the 

rest were all negative. Tne titers for Newcastle 
| disease varied from 1:10 to 1:80. Only three farms 
had a higher titer of 1:80, 10 farms had a titer 
Ofel:H0 15> farms nad a titer of 1520 and" 20 


farms had@ a titer or, mw t0Y in Newcastle disease, 


ae 


most of these are vaccination related titers. The 
: question was what level of titer could be used as 


the best predictor of the condemnation rate (Yo). 











For such a decision, the Newcastle data was 
analyzed using P2R. However, none of the titers 


showed a good correlation with Yo. Because of 





this, the cutoff for a, positivertitersusoius 
for predictive purposes could not be determined. 
Tnerefore, the serology profille data provided 
only descriptive information and was not used 
with the rest of the multivariate data. The 
hemagglutination inhibition test for the sera 
collected from birds at the time of processing 
Was conducted only for Mycoplasm galiisepticun. 
As in the case of the 4 week-serology, the few 
Samples analyzed had negative titers. Additionally, 
due to the fact that information generated from 
Sampling of birds at the time of processing was 
determined to be of limited practical value, 

it was felt that the serology data at processing 
was not critical. Therefore, it was not carried 
outset. 

c. Necropsy profile (ig 927 in an attempt to obtain 
an extimate of the disease rate on the study 
farms, 200 randomly selected birds, at the time 
of slaugher, were necropsied and examined for the 
presence of gross pathological changes. The 
disease and the respective disease rates are 
provided in Appendices 7 and 8. The conditions 
considered were septicemia—-toxemia, airsacculitis, 
synovitis, salpingitis, squamous cell carcinoma 
and necrotic dermatitis, The latter condition, 


was Significant an terms. ofy trimmings for 
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disease parts, which was unavailable from the 





processing plants. Therefore, the necrotic 
| dermatitis rate could be indicative of the 


| importance of trimmed parts in poultry 





condemnations. For comparative purposes, 
the disease rate, lesions, regardless of degree 
| of seriousness which may result in whole bird 
condemnation as well as trimming of parts) based 
on necropsy (X99) and the condemnation rate 
| reported by FSIS inspectors (Ys), were also 
| provided. The rate for X09 varied from 0.114% 
to 5.8% 
6.2 Data analysis to determine the study Mot Chol ce; 
The first major decision to be resolved involved the 
question of which experimental unit to use for the 
lengthy series of data analysis. This was crucial from 
the standpoint of saving time and expenses in computing, 
| as well as to simplify the complexity involved in 
performing analysis on two sets of data - one on a house 
basis and the other on a farm basis. Since comparable 
data, both on a house and farm basis, were collected 
| only up to the end of the 4th brood cycle, this data 


base was used to resolve the question of which study 




















unit (house or farm) to use in the final model 
development. The complete set of data on a house or 
farm basis are given in Appendix 4, Basic descriptive 


Statistics of these are also provided in Appendix 4. 
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6.2.1 Variable Reduction: Based on detailed causal analysis 

of variables which may influence condemnation rate of 

poultry (Y5), systematic screening was needed to reduce 

the potential predictor variables to a more manageable 

level, of about 30-40 predictor variables. 

The variable reduction decision criteria have been 

presented earlier under the research design methodology. 

Following those guidelines, the variable reduction steps 


were divided into 3 phases. 





Ley 


eee ei Phase al, preliminary data screening-using the criteria for 


selection as described under methodology; a large portion of 
the variables were screened out. 
The variables selected at the end of phase 1 were: 

RD (Cee 0 week ose 542 X00. X61, X62-X72, X18 _X83. 


X84, and X100 for hatchery data; 


Neco tee? 2 eX 72) 100 1X20 X30-X33. X39, 
X40, X85-X91, X93, X94, X97, X98, X102, and X109 for 


broiler data; 


X116-X125, X127-135, X138-X148, X164, X165, X169-X174, 
Alii vex tole locees1G0, 41995 X190. X1925 X194, 
X202-X209, and X212 for feed and biologicals data; 


and 


NCPR), (eet, VORLEV Gch SeRG ey -yellne eL am 


X243-X255 for transit and processing data. 


6.2.1.2 Pnase 2 - second stage variable-there were 159 potential 
predictor variables at the beginning of phase 2. Since the 
experimental unit was an individual house, as well as the 
farm, only one of the data base was used to proceed through 
the steps of variable reduction. As it is, a large number 
of computations were involved and performing the analysis, 


both on a house and farm basis, was unnecessary and 


Step 1: 


expensive. The larger data base with more cases and a more 
complete set of information was the farm data and therefore, 
it was used to conduct the variable reduction process. 

To facilitate further data analysis, the data set were 
divided into 4 partitions or submatrices. These partitioned 
matrices were labelled: 

Hatchery data which contained 34 predictor variables, 

Broiler data which contained 42 predictor variabies, 

Feed and biological data with 56 predictor variables, 

and Processing data with 27 predictor variables. 

By using the respective Y5 value with each of the above 
partitioned matrices, the previously indicated sequence of 
BMDP programs which included BMDP2D, 6D, 1R and 2R were run 
and are reported below in the form of steps. 

Tne final set of variables which were selected from the 4 
partitions, combined into one data base to facilitate the 
final stages or predictor model development are given in 
table 8. Note that in table 8, some variables had to be 
inciuded in the final data base simply because these were 
needed to make the retention of other variables meaningful. 
For example, house type 2 (Ky) was selected based on 
selection criteria however, in order to retain Xa7 in the 
final data base, Xags house type 1 also had to be included, 
Je eS ce a7 was an indicator variable, it can only be 
used meaningfully if the associated indicator variables X46 
was also included. Such was the case for hatchery source 
Variables, house type 1, air intake direction, exhaust 


direction, feed source, and vaccine source. At the end of 


step 2, the number of predictor variables were reduced to 45. 
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_ Table 8 
Variables Selected at End of Step 1 
Flock Testing Feasibility Study, 1981-82 
_ Variable Code Variable Name Reason For Inclusion aun 
| Xoo Hatchery Source 1 b 
| Xe 8 Hatchery Source 2 b 
: Xe9 Hatchery Source 3 a 
X Hatchery Source 4 a 
x4 Hatchery Source 5 b 
X53 Hatchery Source 6 a 
Xe Hatchery Source 7 b 
Xe Fertility Rate a 
X63 Hatchability Rate (Hatchery) a 
Kee Hatchability Rate (Brood) a 
X50 Slaughter Weight a 
Xe Strain Type 1 a 
Xo Strain Type 2 a 
/ Xho Strain Type 3 a 
X House Type 1 b 
| 16 
Xa7 House Type 2 a 
) Xig Population Density a 
Xi9 Distance Retween Houses a 
| Xue Air Intake (Roof) b 
| X53 Air Intake ENE Qa 
X54 Air Intake NNW b 
X55 Air Intake WS a 
| Xn Exhaust (Roof) a 
X56 Exhaust ENE b 
X57 Exhaust NNW a 
X59 Exhaust WSW a 
X Mortality Rate at 4 weeks a 
Aaa Mortality Rate at 4-7 weeks a 
X30 Date of Placing a 
Xg7 Length of Time for Preheating a 
X No. of Birds Shipped a 
221 a ; 
Xan No. of Strains on Farm a 
Xi55 Feed Source 1 b 
X 153 Feed Source 2 A 
X oy Feed Source 3 b 
Xn Feed Source 4 a 
X Vaccine Source 1 b 
169 ; 
X70 Vaccine Source 2 a 
X74 Vaccine Source 3 a 
X49¢6 Sanitation Index a 
X139 Litter Type 1 a 
565 Litter Type 2 a 
XS 6 No. Dead on Arrival a 
Gi No. of Inspectors a 
RSS No. of HourssDriving a 


a variable included because it fulfilled selection criteria 
“ variable included to make selected variable under (a) meaningful, 
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Step 2: Subsequently, using the farm data, with the 45 
predictor variables, BMDP1R and BMDP2R programs were 
run. ‘This time the selection criteria for retention 
were: 
- if selected by forward stepping BMDP2R computer 
program 
- if the correlation coefficient of the predictor 
variable vs. Y2 was > 0.1 or < -0,1 
- if the absolute numerical value of the 
"tC" statistic pertaining to the predictor 
variable was > 1.25 
The satisfaction of any one of these criteria allowed 
the predictor variable to be retained for the 
subsequent steps. Those variables selected are shown 
in table 9. At this stage only, 30 variables 
remained. Four predictor variables were deleted for 
the following reasons: 

X45 = air intake direction via roof — there was 
only one case with such a value and so the 
variable was omitted. 

X123, X124 = feed sources 2 and 3 actually 
represent the poultry firms themselves. 
Since the vaiues for these variables were 
too broad and many factors affect the 
performance of the company and also since 
the result could not be generalized to 


other areas, these variables were dropped. 











a ee 








Variables Selected At End of Step 2 


Flock Testing Feasibility Study, 


Var. Name 


Hatchery Source 2 
Hatchery Source 5 
Average fert. rate 
No. Birds Shipped 
Strain Group 2 

House Type 17 

House Type 2 

Air intake-roof 

Air intake—NNW 

Air exhaust ENE 

Air exhaust NNW 

Date of placing 

Feea Source 2 

Feed Source 3 

Vaccine Manuf. 3 

No. of Inspectors 
Hatchery Source 3 
Hatchery Source 4 
Hatchery Source 6 
Hatchability Rate 
Hatchability Rate (brood) 
Now Ole otrasns 
Exhaust WSYW 

Mort. rate at 7 weeks 
Length of Prenheating 
Vaccine Manuf. 2 

NO wot Nes. driving 
No. dead on arrival 
Mort. at 4 weeks 
Farm Sanitation Index 


Selection 
by ©P2k 


yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 
no 


~ ”~ 
CouRreiws 


eae, 


0.2468 


0.3334 


0.3650 
0.3658 
0.1124 
0.360 
-0. 3301 
=(2 1005 
0.3433 
0.1565 
-0.2485 
0.2490 
0.2480 
0.1233 
-0.1570 
-0,2225 
-0. 3615 
0, 1202 


1981-82 


Coeff. 


2 eee 


t{ value 


-1.29 
~1.489 


—2 445 
=1.552 


Be el 


1.640 


1.570 


Ale 
-1.415 
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X186 = farm sanitation index was deleted because 

it was based on partially subjective 

criteria. 
Thus, at the end of step 2 of phase II, 26 predictor 
variables were retained to be used in the final 
analysis to determine which experimental unit to use, 
At this point, the disease rate based on necropsy 


(X49) was added to the data base. 
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Step 3. The final 27 potential predictor variables 
| established at the end of Step 2 were further 
Scrutinized and the following steps were taken. 
Although hatchery sources No. 5 (X5) and No. 6 (Xe 3) 
were useful predictors of condemnation rates in the 
Study area since these sources could not be 
Beneralized to other areas, the indicator variables 
for hatchery sources were deleted. Similarly, the 
indicator variables for air intake and exhaust 
directions (X56 0 Koos X59) were also omitted. Thus 
| only 18 predictor variables were retained. From an 
examination of the residual plots in steps 1 and 2, 
two outliers were identified and those 2 cases were 
: removed, Outiiers are extreme observations whose 
points ile beyond the scatter of the remining 
residuals usually four or more standard deviations 
from zero. At this point though, two cases with 
outliers which showed standard deviations values of 
only 2 or 3 were removed from the data base to see if 
a better fitting model could be devised. For the 


final decision making process using the 18 selected 





predictor variables, BMDP2R was run both for the farm 
and house data. 
Opese) Equatzons for determining the study unit of “cholce: 
Gaseql-_ocudy Unit = Farm 


The preferred predictive model was obtained at the 





end of Step No. 12 at the termination of the forward 


stepping procedure, the predictor model being: 


A 
Nees -2.753+0.038X¢ 540. 009K, -0.256X 3440. 157%Xy 0 


>= 
-0,,003X 9, -0.011X%9740.052X 49-06 275X47 4-0. 125X554 
-0.001X5354+0.095X5y, (a) 
F Value = 7.44 
R = 0.8652 n = 43 
R° = 0.7485 x, = 19 


ve = predicted condemnation rate 
Xoo = avg. fertility rate 
Xay = number of strains on the farm 
X = age specific mortality rate at 4 weeks 
Xn = age specific mortality rate at 7 weeks 
X90 = date of placing of birds in brooder house 
Kea = length of preheating of brooder house 
X109 = disease rate based on necropsy 
X74 = vaccine manufacturer No. 3 
Kooy = number of hours driving from farm to plant 
X535 = number of birds dead on arrival 
Xan 4 = number of inspectors in plant 
n = number of cases 
Xy = number of predictor variables to select from 
R = multiple correlation coefficient 
R = multiple coefficient of determination 
Tne predictive value in equation (a) was high and 74.85% 
of the total variability was explained by the model. Of 


course the predictive variables (X) are indicated in the 
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equation, 


Case 2 - Study Unit = house 


The predictor model was: 


A 
= Q = = 5 
Y5 0.0223+0.281X24 0.081%) +0. 152X449 0.003%536 
Rees Oscc05 nese co 
2 
R- = 0.6998 X = 18 (b) 
Where: 


ve = predicted condemnation rate 
X39 = age specific mortality rate at 4 weeks 
Xno = age specific mortality rate at 7 weeks 
X09 = disease rate based on necropsy 
X5365 = number of birds dead on arrival at plant 
n = number of cases 
x = number of predictor variables to select from 
The respective R and Ro values were 0.8365 and 0.6998. In 
this case, 69.98% of the variability was explained by this 
predictor model. 
Case 3 — Study Unit = farm (via weighted least squares): 
In this third case, a weighted least squares 
procedure was used in the analysis of the data via BMDP2R 
as presented under research methodology. Weighted least 
Squares is a two stage procedure whereby in stage 1, using 
a predictor equation as obtained (e.g. equation a) and 
from this equation, estimates of the weights were obtained. 
The weighted least squares predictor model was: 


a 
Y2 = -0.264 + 0.027 X63 -0.015 X65 -0.142X14-0.267 X16 


la 
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-0.052 X39 + 0.001 X85 + 0.136 X170 ~-0.096 X224 + 
0.051 X24u1 (c) 
(R = 0.8935, R© = 0.7983) 
where: 
A 
Y5 = predicted condemnation rate 
X... = average hatchability rates (%) for hatchery 
X65, = average hatchability rate (%) for brood 
Xay = number of strains on farm 
X = house type 1 
X39 = age specific mortality rate at 4 weeks 
Xgc = age specific mortality rate at 7 weeks 
x = vaccine manufacturer no. 2 
X = number of hours driving from farm to plant 
Xony = number of inspectors in plant 
R = multiple correlation coefficient 
R” = multiple coefficient of determination 
Case 4 — Study Unit = house (via weighted least squares) 
The weighted least squares procedure via BMDP2R for 
the case where the house was the experimental unit was 
handied in the same manner as described for the farm case. 
The result of the weighted least squares was: 
{2 = 17.202-0.067 X62 -0.052 X63 -0.079 ¥55 + 0.184 X17 
-0.123 X39 + 0.003 X85 + 0.953 X87 (d) 


(Res Ose ties Ro = 0.7590) 


<> 
1] 


2 predicted condemnation rate 


Xoo = average fertility rate (%) for hatchery 


























63 = average hatchability rate (%) for hatchery 
X65 = average hatchability rate (%) for brood 


X = house type 2 
X29 = age specific mortality rate at 4 weeks 


Xo0 = date of placing of chicks in brooder house 


X eq = length of time for preheating (hrs) 
R = multiple correlation coefficient 
2 
R= multiple coefficient of determination 


The signigicant point to be made now when examing these 
four final equations is that it is to be remembered that 
there was interest right from the start to ensure the 
utilization of poultry flock in a house as the study unit 
as compared to the flock on the farm, After analysis of 
the data, using both study units, there are some 
differences in the type of predictor variables selected. 
However, when one examines the R° value or the amount of 
variability explained by the respective predictor 
equation, along with the selected predictor variables, the 
gain in collecting data on a house basis has not heen 
beneficial. We concluded therefore, that the data on a 
farm basis provide just as must useful predictive 
information as those based on a house basis and the 
additional effort and difficulties encountered in 
gathering data on a house basis did not justify its use as 
the preferred experimental unit. Instead, from here on, 
the rest of the data analysis was based using a broiler 


farm as the study unit. 
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6.3 


65361 


Data analysis for developing the predictor model: With 


the data base given under 6.1, a sequence of BMDP programs 
were executed in order to develop the final predictor 
model. The complete set of the data used in the analysis 
and the printouts from BMDP analytic steps are given in 


Appendix 4. 


Descriptive statistics of variables used in the study: 


Detailed basic descriptive statistics were obtained using 
BMDP2R. 

For the initial 141 variables, a summary of basic 
descriptive statistics is provided in Table 10. Two 


variables of interest, X and Y the disease rate 


109 apht 
based on sample necropsy and average condemnation rate due 
to diseases observed at processing, served as the response 
or dependent variables in the predictor model, and were 
therefore examined closely. 

For Xi99° there were 75 cases (n=75) available for 
analysis. The mean value of the disease rate for all 
farms and broods was 2.7%, with a range of 0-11.1% and a 
variance of 4.9. The histrogram of the variable was 
skewed to the right. 

In the case of Yos where n=137, the mean condemnation 
rate at processing was 0.68%, with a range of 0.08% to 
5.8%. The variance and standard deviations were 0.38 and 
0.62 respectively. The histogram of this variable was 
also skedwed to the right, indicating a non-normal 


distribution curve: 
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table 10 Descriptive Statistics for Variables Used In Flock Testing Feasibility Study 





Variable Sample Mean, (3) Standard Minimum Maximum 
Size (n) Deviation 
(SDS) 
X57 138 0.34 0.48 0.00 1.00 
X58 138 0.19 0.40 0.00 1.00 
X59 138 0.19 0.40 0.00 1.00 
X260 138 0207 Ot27 0.00 1.00 
X261 138 0.23 Or47 0.00 1.00 
X262 138 0.23 yey 0.00 1.00 
X263 138 0.07 Utz) 0.00 1.00 
X60 138 OF57 0.50 0.00 1.00 
X61 138 O.00 0.49 0.00 1.00 
X62 jx) 9077 2 Pia. 89.00 94.00 
X63 123 83.82 1.82 80.00 87.60 
X64 wen 90.91 bisa 89.00 94.00 
X65 T3532 83.93 2.48 79.80 88.60 
X66 139 3.34 0.84 3.00 4.00 
X67 139 3330 0.83 2.00 5.00 
X69 139 23:38 byes ie! 13300 25.00 
X264 139 0.96 O19 0.00 1.00 
X78 139 O25 0.42 0.00 1.00 
X79 139 0.46 0.50 0.00 1.00 
X80 139 O22 3 Oy 0.00 1.00 
X265 139 0.42 0.50 0.00 1.00 
X81 139 O723 0.42 0.00 1.00 
X82 139 0.84 0.36 0.00 1.00 
X83 139 O19 0.40 0.00 1.00 
X84 ESY 92.50 46.49 33.00 230.00 
x3 138 35145296 1579 ea? 14858 .00 63818.00 
221 137 34335.69 15480072 14760.00 63028.00 
X222 sl Se ey 0.46 7770 4.30 
X7 as ee 6.61 6.82 0.00 23.00 
X8 138 OnAZ 0.50 0.00 1.00 
X9 138 0.50 0.50 0.00 L206 
X10 138 0.57 Grou 0.00 POG 
X14 136 Wa Pe 1.09 } .00 4.00 
X15 t37 3.65 1.64 1.00 8.00 
X16 138 0.84 O20 0.00 1.00 
X17 138 0.19 0.40 0.00 1.00 
X18 138 1.36 0.19 O796 1.96 
X19 139 44.69 Zo e039 0.00 70.00 
X20 138 0.84 0.36 0.00 1.00 
X22 137 edd) ead iS} 0.00 38.00 
oe 139 0.88 Cr32 0.00 1.00 
X24 139 0.96 dedi 0.00 £200 
X25 139 0.88 O 232 0.00 1.00 
X26 138 0.26 0.45 0.00 1.00 
X27 138 0.34 0.48 0.00 1E0C 
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Table 10 (cont.) 


Descriptive Statistics for Variables Used In Flock Testing Feasibility Study 


cont. 
Variable Sample Mean, (x) Standard Minimum Maximum 
Size7(ny Deviation 

eee. oe pee (cee) ee 
X28 138 0.19 0.46 0.00 1.00 
X39 136 2.54 0.87 0.95 4.23 
X40 128 3.35 1.08 1.80 5.63 
X85 137 112.34 47.42 53.00 251.00 
X87 133 13.50 7.06 T2200 48.00 
X88 139 0.38 0.49 0.00 1.00 
X89 138 50.69 Beil 34.00 53.00 
X91 136 16131. 46 4497.68 9791.00 31660.00 
X109 75 phek ss 2.1/0 0,00 11.05 
X116 138 0.80 0.40 0.00 1.00 
X117 138 0.84 0.36 0.00 1.00 
X119 138 0.82 Oeat 0.33 1.50 
X122 139 0.42 0.50 0.00 1.00 
X123 139 O223 0.42 0.00 1.00 
X124 139 0.26 0.45 0.00 1.00 
X266 139 0.07 Of27 0.00 1.00 
X125 139 4.80 0.80 4.00 8.00 
QU ra7) 2 317 376 Li2202 70.00 600.00 
X129 134 0.54 0.33 (ae 5) Lsic 
X130 139 8.42 1.41 4.00 10.00 
X148 138 40.73 9.99 2100 53.00 
X169 138 0.92 O27, 0.00 1.00 
X170 Sw 0.50 0.50 0.00 1.00 
X171 137 0.53 0.50 0.00 1.00 
X186 110 Lhe 3! Oral 1.00 Ze 
X189 134 0.76 O42 0.00 1.00 
X267 Nels} 0.19 0.40 0.00 1.00 
X190 116 2205 0.92 1.00 4.50 
X202 136 77.39 8.99 62.60 90.80 
X203 136 Oeau 8.96 ot ba) 28) 80.90 
X204 L35 53 a0 9.26 39.20 71.20 
X205 136 405 1.09 1.60 a Fee, 
X206 139 0.03 0.19 0.00 1.00 
X207 139 0.00 0.00 0.00 0.00 
X208 139 0.88 0232 0.00 1.00 
X223 129 32564 26 4314 4.00 80.00 
X224 129 0.87 0.61 0525 2.00 
X231 7S 745 tas 2333 eL0 5600.00 15850.00 
X235 131 81.69 69.80 10.00 336.00 
X237 120 LaO7 0227 1.00 oe) 
X24] 134 10.26 2056 6.00 12.00 
X245 135 Omg2 0827 0.00 1.00 
X247 [ahs 0.42 0.50 0.00 1.00 
X248 135 0.23 0.42 0.00 F.00 





Table 10 (cont.) 
Descriptive Statistics for Variables Used In Flock Testing Feasibility Study 


cont. 


Variable 


Sample 
Size (n) 


Mean, (x) 


° 


MWoOonnownwwor dO hd 


NODCGCOWDOONTOCO 
OnNUUNnNF ODOAN OD 


13.13 
135453.03 
123698.18 

368.47 
11.62 
1.45 
0.00 
0.00 
0202 
0.00 
0.00 
0.00 
1.06 
0.19 
0.68 
Geol 
0.93 
0.01 
0.32 
OeZ2L 
0.13 
0.03 
0.06 
0.01 
O235 
UsZZ 
0.01 
0.04 
0.06 
1200 
Oe 21 


Standard 
Deviation 


Cons) 


0.45 
0.27 
2.20 


W NN = 
© Uo G1 c= 
me 


40365.80 
43530.83 
£00259 
3.95 
0.50 
0.00 
0.00 
0.02 
0.02 
0.00 
0.00 
O72 

0 ams 
0.42 
Leu? 
0.78 
0.01 
Use 
Onsz 
Oe U2 
Uru? 
0.06 
0.05 
0.25 
0.50 
0.02 
0.03 
On0F 
bale 
0.42 


Minimum 


0.00 
0.00 
Oc55 
0.00 
0.00 
4250 
~0.00 
0.00 
0.00 
0.00 
1.00 
1.00 
1.00 
1.00 
4.00 
81000.00 
47166.00 
250.00 
8.00 
1.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.02 
0.19 
Or22 
0.18 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.07 
0.00 
0.00 
0.00 
0.00 
0.22 
0.00 


Maximum 
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Table 160 (cont.) 


Descriptive Statistics for Variables Used In Flock Testing Feasibility Study 
cont. 


Variable Sample Mean, (x) Standard Minimum Maximum 
Size (n) Deviation 
CSiDe) 
Y9OT 134 0.02 0.06 0.00 0.57 
Yior 134 0.39 0.29 0.00 1.65 
YuEE 134 0.27 0.56 0.00 5.02 
YA 134 0.01 0.03 0.00 0.38 
Y 137, 134 0.03 0.03 0.00 O.15 
Y14T 134 0.08 Osb2 0.00 0.80 








Among the disease specific condemnation rates 
(YoY ay) the mean values for all study farms during the 
1981-82 calendar year were 0.01% for leucosis (Yq)s Ge354 
for septicemia (Yoo). Geeep for airsacculitis (Ya4)s 
0.0095% for synovitis (Yyo)> 0.036% for tumors (Y43)s and 
0.056% for others (Yip de This showed that in terms of 
importance, as observed in other national figures, 
condemnation due to septicemia/toxemia and airsacculitis 
accounted for 83.3% of a1] whole carcass condemnation due 
to diseases. In terms of emphasis of developing disease 
Specific condemnation rate models, these two categories 
will require emphasis as opposed to the others which, at 
this point, do not appear to be as important. 

In an effort to characterize the relationship of the 
independent and dependent variables, each of the 
independent variables were plotted against Yo and Xi09° 
using BMDPSD. A scatter diagram of each independent/ 


dependent pair is provided in Appendix 4. Amongst the 


pairs evaluated, X VS Y., showed the strongest linearity 


109 2 


with a correlation coefficient of 0.5888, followed by X502 


VS Y, and Kooy VS Ys which respectively had correlation 
coefficients of -0.3868 and -0.3720. 

Variable X503 represented the distance from the 
poultry farm to the processing plant, while Kooy 
represented the number of nours for driving from the farm 


to the processing plant. The rest of the variables did 


not show strong linearity. 
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De SIOe 


When X 109 was the dependent variable and the other 
independents were plotted against it, the variables with 
adequate correlation coefficients were the condemnation 
rate at one previous time period (Y5.) with a value of 
-0.2239, the average fertility rate for the hatchery (X65) 
with -0.2185, average hatchability rate for hatchery (Xe) 
with -0.2239, and average hatchability for the brood under 
study (X-.) with -0.2182, The rest of the variables did 
not show satisfactory linearity on the scatter plots. 
Comparative statistics to assess differences in 
condemnation rates among study units: To address the 
question of whether there were differences in the 
population means of whole carcass condemnation rates at 
processing, a number of fixed factor effects were 
considered. These included mainly, an assessment of 
differences (by place) among study firms, farms, 
hatcheries, or (by time) amongst the five brood cycles, 
Such assessments were performed graphically (figures 15) 
Via tables (11-14) and finally, using analysis of variance 
(two way layouts). 

Figure 15 provides an illustration of the profiles of 
average mortality rate of poultry and the corresponding 
average condemnation rate during the five brood cycles and 
the respective flock sizes. The average mortality rate in 
the study farms/firms appeared comparative from brood to 
brood. This also was the case with the average 


condemnation rate, except during the fourth brood cycle, 
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Table 11 


Comparison Of Condemnation Rates And Disease Specific Rates 
For Study Firms 


Youu | peaay9 Mei UC aul Vi? are Y14 
Firm 1 
(All Broods) 0.723 0.029 0.368 0.132 0.01 0.049 0.136 
Firm 2 0.259 | G200 5m eee 4 0.01 20a alomeeaod 0.008 0.003 
Firm 3 1.014 | 0.014 | 0.287 0.599 0.028 0.058 0.047 
Firm 4 10.728 0.008 0.540 0.112 001 0.022 0.053 
Firm 5 0.499 0.004 0.370 0.085 005 0.038 0.0 
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Table 14 
Condemnation Rate (Y7) For The Strata Of Study Farms Used In 
Flock Testing Study 
| Farm/Firm Good Piverace Welw .WihdRéoray 
| A 0.649 0.280 1.948 0.377 1.302 0.460 
| 0.631 0.938 teOLor "0.722 0.672 0.494 
| 0.648 0.529 0.962 0.936 0.463 0.700 
| Peee7e 09747 0.508 0.691 0.747 1.020 
| 0.720 0.603 0.542 0.888 0.425 0.701 
| 0.450 
| 0.464 
| 0.433 
| 0.621 
| CAG OgS76E 9 Pe 
B 0.138 0.243 0.255 01226 O220LER0 21) 
0.212 0.365 0.232 0.367 Or 1820 70"312 
Ovi 0.146 0.248 0.288 0.249 
0.543 O22 219 OFTT4 0.137 
0.344 0.259 0.368 
0.084 0.559 
| 0.231 0-218 
0.239 
D 0.934 0.674 1,037. 10023 2.044 0.331 
0.714 0.579 0.658 1.002 0.488 0.868 
0.662 0.697 0.896 0.982 OF40 18202478 
17 7en0 730 24/039 00 4797 Onsi7 
5.800 12207" 0: 618 0.865 
0.902 0.477. 0.449 0.960 
0.617 
1222 
| ern ron Pee en ee ee 00905 
E G51 ie 01214 0.601 0.380 
0.444 0.423 0.780 1.477 
0.686 17692 TAIN734 None Available 
| 1.115 1.319 0.754 
| 0.314 0.940 0.326 
RT 1.209 0.232 
es E 0m5503 0" 0.210 0.220 
/ 0.368 0.263 0.382 
2.033 0.330 0.414 
0.324 Oa? 0.312 
0.470 0.295 0.836 
0.283 0.598 
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when it was lower. Note that the comparisons could be 
made along the norizontal axis, whereby the respective 
values for each of the five firms and five brood cycle are 
presented. Tne five sets of graphs could also be examined 
vertically, whereby possible differences during the 
different brood cycles could be examined for eacn firm, 
one at a time. 

In figure 15, firms A and D showed a higher 
condemnation rate than the others. The condemnation 
rates, however, varied from brood to brood and from one 
study farm or firm to another. In Table 11, the 
respective percent (%) condemnation rate due to specific 
diseases is provided for each firm for comparative 
purposes, 


Table 12 presents the values for Y,, the condemnation 


>? 
rate for each of the study hatcheries and the 
corresponding brood cycles which also showed variabilites 
both by time and place. Table 13 provides comparative 
Values for selected hatchery parameters. Other than 
providing somewhat superficial profiles of the response 
variable (Y.)s indepth analysis and assessment of the mean 
differences for this variable could not be performed, 
Therefore, such variabilities were examined 
quantitatively using a fixed effects, two way ANOVA 
models. Tne tested hypothesis to examine differences in 


condemnation rates were: 


- test of equality of row means of condemnation 











eee el — 
SASS A 


rates (Y,) 


- test of equality of column means (Y,) 


- test of no interaction between column and row 


factors, 


The following cases were considered: 


a) 


to assess if differences existed by place (firms) 
and by time (brood cycle) the average 
condemnation rate (Y.,) was classified for the 
two factors, i.e. five levels of firms (A, B, D, 
E, F) and brood cycles, which had three levels, 
(Cycle 3, 4 and 5). An analysis of variance was 
performed (Table 15) and an F ratio of 4.13677 
was obtained for the firm factor indicating that 
there was a statistically significant difference 
in condemnation rates amongst the five firms at 
the 95% level of significance (F 0.95, 4, 8 = 
3.84). Conversely, there was no statistically 
significant difference in terms of brood cycles 
or in interaction between the two factors. 

study farm (4 levels) by study firms (5) levels: 
The condemnation rate (Y,) was cross~classified 
by study firm and farms. The rows represented 
firms and columns represented farms. Using a 
two way ANOVA approach, an F ratio of 3.7007 was 
obtained for the row factor (F 0.95, 3, 27 = 
2.96), while the column factor and interaction 


Showed a non-significant F value. 
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Table 15 Anova Summary Table for Comparing Performance of Firms 


Source 


Brood Cycle 
Firms 


Brood x Firm 
Interaction 


Error 


TOTAL 


76 


90 


° 107542 
Heo L253 


2 993236 
- 365631 


-4311 


F-Ratio 


> 294126 
4.13677 


beO225 


























{eS 


c. A similar analysis was conducted for firms A, B 
and D, but, with the brood cycle factor occuring 
at 5 levels. For this case an F ratio of 
13.3521 was obtained. This was also 
statistically significant (F 0.95, 2, 8 = 4,26). 

d, Flock size (X3) vs disease rate on necropsy 
(Xi99)5 the variable x. was divided into three 
categories viz. farms with (X3) flock size of: 
20,000 or less, 20,900 to 40,000 and those with 
over 40,000 birds. The disease rate (X19 9) was 
GividedeintonOve vom less, (Ove: to 1,0 and greater 
than 1.0. The result indicated that the F value 


for row (X.) and a column (X ) and 


3 


interaction, were all insignificant, 


109 


It is coneluded then that, for the firms studied in 
this project, both from the graphical analysis as well as 
the results of ANOVA, there were significant differences 
in the means of the samples for whole carcass condemnation 
rates among firms/farms. However, there were no such 
observable or quantifiable differences over time, f,e, 
(five brood cycles). 

It is important to emphasize here that the ANOVA was 
performed to provide more comparative information among 
Study units. Tne basic objective of the study was not 
comparative in nature, but one which required the develop- 
ment of a predictive model for condemnation rate. Thus 


multivariate causal models were utilized for such a task 


Geses 


Geer 


and these will be presented next. 

Multivariate analytic methods: The analytic methodology 
to develop the predictive models for poultry whole carcass 
condemnations due to diseases were based on linear 
multivariate models which involved multiple linear 
regression and discriminant analysis. The rationale for 
such an approach and the diagrammatic visualization of the 
static linear causal model which interrelated the 
predictor variabies for condemnation rates (Y) have been 
given earlier. (Figure 12, 13). 

For these two models, systematic variable reduction 
steps to select the best subset of predictor variables, 
were required. Following the decision-making criteria for 
variable reduction/deletion (5.10), the data base was 
examined critically, in steps, to arrive at the final 
conciusion. 

Variable reduction phases: 

Phase I - Preliminary — Based on criteria given under 
5.10, from the 286 original potential predictors, 162 
variables were deleted and only 124 independent variables 
were retained. These variables formed the initial data 
matrix upon which the various BMDP statistical programs 
were run. A listing of these initial variables (X) and 
“the dependent variables retained at the end Phase T, are 
given in appendix 9. Although, the size of the variable 
list was substantially reduced, this still left a large 


number making further reduction necessary. 


154 



































153 


Pnase II —- At this point a series of BMDP programs were 


used specifically to facilitate variable screening. The 
criteria for selecting and retaining predictor variables 
during this stage has been preserited under imnethodology. 

Unlike the description presented under 6,2, where the 
analysis was performed using four partitioned matrices, 
during this last phase of the project, the BMDP program 
was already modified to handle as many as 200 variables 
Simultaneously. Therefore, no partitioning was required 
and all the variables, whether for descriptive purposes 
(P2D, P6D) or regression analysis (PIR, P2R), were handled 
in one sweep. Specifically, in the case of PIR and P2R, 
this meant that all potential predictor variables with 
some infiuence on the dependent variables (Y) could be 
evaluated either individually and/or in combination with 
others at the same time. Utilizing this method, those 
With the least influence could then be weeded out with 
ease, efficiency and with minimal computing cost. 

The latter term is used purposely to emphasize that 
with 4 partitioned matrices, for example, P1% would have 
had to be run 4 times (one run per partition) while with 
the non-partitioned case, only one computer run was 
necessary. Tnis reduced the cost of computation and the 
time involved in performing the complex set of analvVScs. 

Briefly, during phase II, the RMDP analysis strategy 


consisted of two steps. 


step 1 COs eUeUR Ol 


a) Running P2D for all the variables, 


b) Performing P6D and obtaining scatter plots for 
each of the variables: Xeom* ay? x eae iit vs 


Xi09 (the disease rate based on necropsy), 
each of the variables Xo “Tone vs Yo (condem- 
nation rate). 

Only variables which were non—dichotomous 
were used in these scatter plots. 

¢) conducting PIR and P2R for all predictor 
variables (X 


- Y,) with X as the dependent 


a 109 


variable. When X409 Was in use, the disease 

specific rates corresponding to this variable 

(Xi 4 “Xi45? were omitted from the data matrix. 

Similarly, the above steps were performed with Y, 

as the dependent variable. 

Step 2: Subsequent to the results obtained from step 

1, another set of P2R were run for each dependent 
variable, to further reduce the number of variables. At 
this stage, the predictor variables selected for retention 
were based on the results of the stepwise forward and 
backward stepping sequence of P2R. Tne final potential 
predictor variables were further scrutinized at the end of 
step 2 and the following steps were taken. Although the 
hatchery source variables were useful predictors of 
condemnation rates in the study area, since these sources 
could not be generalized to other areas, hatchery source 


variables were deleted, Similary, the indicator variables 
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for air intake and exhaust directions in the hatchery or 
broiler nouse were also omitted, 

Therefore, at the end of of phase II, the variable 
reduction task was completed (table 16). Using the 
predictor variable provided in table 16, two linear 
multivariate models were developed, one for predictive 
purposes using BMDP2R, and the other using BMDP7M for 
Classification purposes. 
separation of cases into study and test data sets: In 
order to develop the predictor models and then test or 
validate the performance of the model, the flock testing 
Study data base was divided into two sets. One of the 
sets on which the models were developed was referred to as 
the "study" data base, and the other one as the "test" 
data set. The "study" data base was composed of all cases 
in the file, with a complete set of records for each of 
the variables used in the analysis. The "test" data set 
were those that had an incomplete set of records but could 

still be used to validate the final model when a very few 
set of predictor variables were in use. The point is 
that, when the number of predictor variables In use were 
large, the number of cases with a complete record were 
restricted while if the number of predictors were samll, 
then a larger number of cases were regained. For example, 
there were a total of 154 cases in the data base in the 
computer. ‘When developing the predictor model for Yos at 


the end of phase I, there were only 26 cases with a 


ey 


Table 16 
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Variables Selected For Use In Developing Final Predictor Model a,b 


Dependent Variable 


07, 


YQ 


wv 


Y11; 


2G) 


Y13 


Yy4 


X109 


Number of 
Predictor 


Variables 


14 


25 


14 


18 


16 


19 


Predictor Variables 


X60> X61> X62. X63» X64» X14> 
X109» X129> Xa0q4> X205> %223> 
TT's X99, X18 


X17, X22, TPOP, Xj99, X119> 
X205> Yor, Y10T>» X%237 


X17» X22, TPOP, Xj09, X119, 
X905> Yor. Y1i0T, X237 


X61> X62, X63> X65, Xg5» X88, 
X116> X119> X127> X129> Xj g9> 
X267> X208> ZR» Xo» Xy4> X15> 
X19, X49» X91, X130. X205, Y11T 


X64, X65, X66, X7, Xy4> X19; 
Xgg, X125, X129> X148> X133> 
ZR, X67» X140 


X625 X63 Xg5» X85» X88, X119> 
X1255 %202> %224> %235> X245> 
X133> ZR. Yau» Yi3p> %X117> %190> 
X206 

X61, Xg4, Xg, X19, X40. X88, 
X109, Xy43> X189> X267> 202; 
X905> X208> X224> X237, Yar 


X60, %63> %65> X222, X16, X17> 
X39> X40. Xgg> X1175 X148, X189> 


X190> X208> X206> X231> %133> ZR> 
TROP 


\ 








fable 16 


cont. 
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Variables Selected For Use In Developing Final Predictor Model a,b 


Dependent Variable 


Xii1 


X112 


X113 


X114 


X115 


Number of 
Predictor 


36 


13 


23 


16 


18 


Variables 


Predictor Variables 


X575 X58, X262, X63, X67> X265, 
Xg2, Xg4, X17> X27, X85» X116> 
X117, X169> %202, X203> X204> 
X205; X2355 X133> 29, Y2T, Y1oT> 
X59, X264> X9, X14, X22» X397> 
X91, X119> %130, %189, X206> 4R> 
TPOP 


X62, X63- X64, X16, X119> X189> 
X267+ 4R> X61 > X75 XQ, XA0> X148 


X260, X60, X61 X64, X65> %67, 
X264> X82, X7, X9, X19, X22, 
X23, X27, X28, X40, X85, X116> 
X129> X169, X205> 29, Y12T 


X65 479, X15, X23, 4125, 4129; 
Xj69> X172> X%186> X208, 223, 
MID ON950 20127 > OX 13T> TEOE 


Xg4> X14, %i5> X16, X23. gi; 
X148> X%369> %170> %267: %202> 
X205> %206> X%937> 245, 29, Y1aT» 
Seon 


@ The dependent variable X 110, disease rate due to leucosis based on necropsy 
had a value of 0 and therefore, not used in this analysis. 


b TPOP = 


“ne 7 Gy » where: 


X3 = population of poultry on each study farm 
X3 during the first week 


X92] = population of poultry of each study farm 
at the end of 7 weeks. 


Sees 5 


complete set of records. At the end of phase 2, there 
were 37 cases. After the final subset of predictor 
variables were obtained, there were 17 additional cases 
which were used to validate the model. 

Although, the original intention was to divide the 
data base into study and data sets randomly, this approach 
was ignored since the number of cases with a complete set 
of records for all the variables in use were limited, The 
lack of random allocation of cases should not be a 
drawback since in the system adopted here, there was no 
bias in allocating cases to "study" or "test" data sets. 
This simply was the result of restrictions imposed on the 
case records themselves. The validation task therefore 
was performed satisfactorily as planned as is discussed 


further below, 


Attempts to improve the model: Upon completion of the 


selection of the variables to be used in the final 
preaictor model, two techniques were considered to improve 
the accuracy of the model. These were: 

(a) transformation of predictor variables 

(b) weighted least squares procedures 
Transformation of variables: after examination of the 
scatter plots of independent variables vs the dependent 
Variables obtained from BMDP6D, as well as the correlation 
matrix from BMDPIR, selected predictor: variables were 
transformed into their respective square root and 


ia! e 4 faik 
exponentlale forums.  ihneseavar ables were X65 63° Xoo 
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6.3.4 


“201s Xin» X39+ Xuo» X95, X87, X109, X2au, and X35, 

The square root transformation of these variables did not 
substantially improve the absolute values of the 
correlation coefficient. In the case of transformation to 
the respective exponential forms, there was a substantial 
improvement of the correlation coefficient of X39 from 
0.0068 to -0.1559,. There was a slight increase for the 
transformed value of Xi09 too, from 0.4146 to 0.4848. In 


the case of X although, its exponential form had a much 


39” 
better correlation coefficient, the sign of the 
coefficient had changed from positive to negative, From a 


biological aspect, it was feit that the relationship 


between X39 and Y>5 would be direct. The above described 


cases were not encouraging enough so as to pursue this 
time consuming task on all other variables. Therefore, 
the variable data were used directly without 
transformations. The second aspect of attempting to 
improve the accuracy of the model was the use of weighted 
least squares procedures. This holds a significant 
position in the data analysis and is covered in detail 
further below under its own subheading. 

Predictor models: The final selection of the "best" 
subset of predictor variables was accomplished via a 
computerized stepwise selection process which utilized 
forward selection followed by backward stepping regression 
analysis techniques. This, as referred to earlier, was 


obtained via BMDP2R, with options for providing covariance 
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en gee ean ilewen Ne a8 3 c) H eee y ek iyo 
ana correlation matrices, residual analysis pilots ann various 


tabular cate useful for evaluating the adequacy of the fittec 


rhe results are presentec as cases for varions tynes 
of soncemnation rates viz. 2n overall condemnation rate 


due to diseases and disease specific condemnation rates, 
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The variables selected for such a final step w 
quantitative in nature or were such that the resuults would 
be generalizea to refer to other regions, Conversely, 

a more bocalizea result for the study rection a filam) es 
also previded. In the latter case, Such Vanieuleseas 
Bavenery sources, processing nlants | ana) ouner er pac ome 


only to Alabama were used in oeveloping the final medels, 


Preaicteor models for everaill wates:  Inverder iia eevo71e 


¢ 


reaundancy anda reduce the number of mathematical equations 


USeG “ih GALS Sections the fol lowiss = suiiels eeg wears 
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regression coefficients, will be presented only 
for selected cases, 
d) The description or definition of each variable 

used in the final analysis (X,Y) is in appendix 9. 
The strategy in the development of the predictor models, 
as well as the presentation of these models further below, 
revolved around: 

a) developing a predictor equation for condemnation 

rate (Y,) or the respective disease specific rates 

(YoY yy) using X 1099 (disease rate based on 

necropsy) as one of the predictor variables. 
One of the important phases of this project was to utilize 
the available data in order to evaluate the predictive 
Value of X109. After analysing the data both on a house 
basis as well as on a farm basis, it was found that the 
predictive vlaue of X109 was significant. The variable 
had a strong linear relationship with Y2 (the condemnation 
rate at processing), and its coefficients (0.109 on a 
farm basis and 0.147 on a house basis) appeared relatively 
stable. This indicated that this variable could be used 
as a predictor if a practical alternative of sampling 
birds on the farm 1-2 days prior to shipment to slaughter 
could be devised, Since the concept of flock testing 
partly relied on the result of the performance of this 
Variable. it was felt that this finding was very 
Significant. 


It is again vital to briefly touch on the idea of why 


Oe 


164 
flock testing and the use of X109 as a predictor is a 
rationai adjunct (conceptually) to poultry inspection by 
referring to figure 3. Other than processing plant 
errors, the major determinant of condemnations and the 
area of responsibility for inspection is that of 
condemnations due to diseases. In this case, since it is 
highly unlikely that diseases suddenly arise in transit 
or while birds are in processing, it is only logical to 
think that the disease conditions that led to 
condemnations during inspection must have had some time to 
develop and be recognized as lesions. In that case then, 
this must have taken place while the flock was still on 
the farm, Although the disease rate (variable X109) was 
determined at the time of slaughter, it could be inferred 
that X109 wouid be representative of the disease rate 1-2 
days earlier while the birds were still on the farm. 
Therefore, if a mechanism could be devised to obtain such 
information prior to slaughter of the flock, a reliable 
indication for the condemnation rate could be projected. 
b) Although X109 was a significant predictor, the 
question of how to devise an approach to the 
random sampling of birds while the birds are 
Still on the farm was raised. One alternative to 
this problem was to also devise a predictive 
model for X g as the dependent variable. This 


10 


way, one first predicts a value for X for a 


109 


broiler flock to be processed using available 




















data as prescribed in an equation for X1099° The 
predictor equation for Y5 which would have X 109 


as one of its predictor variables will then be 


replaced by the predicted X value, and an 


109 
estimate of Y> Will be computed. Based on Yo» 
decisions may then be made for an appropriate 
type of inspection. 

ce) Weighted least squares equations were used to 
improve the predictor equations for toy YoY oy 


and X409° 


The results of the final analysis provided the 
following predictor models presented under cases 1-3. 


Tnese equations therefore were the end products of the 


total effort of the project carried out during 1981-82. 


Case 1. Predictor models for study area: 


Since the study was centered around the poultry industry 
in Alabama, some of the variables studied were unique for 
the region, e.g. hatchery sources. Therefore, the first 
objective was to present the results of the predictor 
equations pertinent to the study area. The equation for 


predicting condemnation rate in the study area was: 


Same », 
0,45324X,ye+ 


0.03781X 544 


a 
Nee 16 .48972-0. 18201X_ 340, 22468%,, 


2 9 


(ommnorao PRG loi7o UPAPoti® 4030.83)e0 (A) 


Where: 
Can . 
Y. = predicted (%) condemnation rate 
X.. = average hatchability rate (%) for hatchery 
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Ps 
' 


= average disease rate (%) based on necropsy 
- 4 2 a 

X595 = total precipitation (in) 

Xony = number of inspectors in plant 


R = multiple correlation coefficient 


> 2] 
oe 


multiple coefficient of determination 
Correspondingly, the equation for predicting the disease 
based on necropsy was: 


A 
Xi09 = 39 54751-2..91267X, 9-2, 02945X -0,44252X 


261 
gt 129965 1x 


65 


+3.23479X._-0.55159X #1657413% 
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255326 


27 3 


-1.8 704K, , 


TE 


+1,04076X as} Cle 


189 8 20671 +9 
36482 2 
+9, 36 Beier (2) 


(Reen0293; Ro =EOROO Erte LLOR= Siocon 


X109 = average disease rate (4%) based on necropsy 
X-, = hatchery source no. 2 

X56 = hatchery source no. 5 

X65 = average hatchability rate (%) for brood 

X = location of air exhausts in broiler house 

(NNE) 
X29 = age specific mortality rate at 4 weeks 
Xx = type of feed (crumbles) 
169 = vaccine manufacturer # 1 
xX = litter type 1 
x = season (spring) 
X5o6 = Season (fall) 
Ze = location of air exhausts in broiler house 


(SSE) 











eooD = change in population size from first week to 

week of brood 

Two equations are presented above: one for predicting 
the overall condemnation rate (Y,), and the other for 
predicting disease rate (X199)> For any predictor 
equation, it is essential to examine the types of 
regression coefficients including their signs and 
magnitude. In this case, the signs of these coefficients 
indicated that X63 and X5Q5 had an inverse relationship 
with the condemnation rate (Y,), e.g. if the average hatch- 
ability rate for the hatchery increased, the condemnation 
rate decreased. This is consistent with the biological 
reasoning that a hatchery of a firm with good management 
would lead to a lower condemnation rate at processing. 
Conversely, X 109 and X,,, had a direct (positive) 
relationship with Y. For example, a high disease rate or 
X 109 would correspondingly iead to a high condemnation 
rate (Y,) at processing. Again, this is consistent with 
what one would expect with an understanding of poultry 
population health. A high morbidity rate (Ky 99) would 
necessarily ‘lead to a higher condemnation rate. In the 
case of X,,,, which represents the number of inspectors in 
a processing plant, the interpretation is two fold. In 
one case, this means that if there are more inspectors in 
a plant, they may tend to condemn more due to the fact 
that they could identify more diseased birds, or, they 


could make more errors in condemning birds which should 
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have passed inspection, In both cases, an inerease in 


Xony leads to an increase of Yo° 

To evaluate the importance or magnitude of each of 
the regression coeffieient in equation 1, it is useful to 
refer to the corresponding standardized regression 
coefficients wnich were: 
+0 ,267X (3) 


= “00415 .X 6240 .554X 44 =U5029% 


9 205 24] 


Y, = predicted (%) condemnation rate 
X63 = average hatchability rate (%) for hatchery 


X19 = average disease rate (%) based on necropsy 


os 
4 


= total precipitation (in) 

Xony = number of inspectors in plant 
In this case, direct comparison could be made amongst the 
predictor variables, and the one with the highest 
numberical value which would be considered as the most 
important; the one with the least numerical value, as the 
least important. Therefore, amongst the four predictor 
Variables, the rankings in terms of importance from 
highest to lowest were: Xi09° XoQ5 a5 and Kony? The 
effect of X and X on Y. are obvious; however, the 
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importance of X the total precipitation during the 


205" 
brood cycle having an inverse influence on Ys may not he 
easy to explain. One explanation may be that rain 

cleanses dust particles and the associated organisms in 


the air so that the risk of exposure to such agents is 


minimized. Additionally too, during the high rainfall 


months in Alabama (late spring and early summer) the 
combined effect of high precipitation and high temperature 
may have adverse cffect on viral and other agents of 
concern in poultry production. Therefore as the total 
precipitation increases, the condemnation rate decreases 
rindicating that the seasonal effect of rainfall on 
poultry production is an important consideration. 

Other pertinent statistics are the multiple 
correlation coefficient (R) and multiple coefficient of 
determination (R°). Although R measures the 
proportionate reduction of total variation in Y5 
associated with the use of the set of X variables, a large 
a does not necessarily imply that the fitted model is a 
useful one. In general, though, it is an acceptable 
measure of the performance of a model in that, the higher 
the value of R° or R, the more useful of the model. In 
equation 1, an ays 0.79 indicates that 79% of the 


variability in Y. could be explained using the four 


2 
predictor variables given. This is quite significant, and 
when examined the corresponding R value of 0.89, which 
indicates the magnitude of the correlation between Y, and 
the four multiple variables, the predictor model appears 
very satisfactory. Note that one attempts to obtain an R 
value close to 1.0 or one tries to explain 100% of the 
variability in the dependent variable. However, in 
biologically complex processes such as the case studied, 


an Ro of 0.79 may be considered very satisfactory. 
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In light of the above discussion, it was concluded 
that Yo. the condemnation rate of poultry at processing, 
could be predicted satisfactorily. Ina later section, 
the validity of this model was evaluated via residual 
analysis and using a test data set. 

In general, adding more independent variables in a 
model inflates Re, and it is sometimes useful to utilize a 
different measure which recognizes the number of 
independent variables and adjusts for it. Tnis is the 
adjusted coefficient of multiple determination, In 
equation 1, the predictor model for Yoo R° was 0.79, and 
the adjusted R” was 0.77. It should be noted that only 3 
independent variables were used for predicting Yo: In the 
case of equation 2, where 12 independent variables were 
usea to predict X109" Ro was 0.85, while the adjusted Ro 
wes 0.81. Although there is a slignt decrease in the R? 
Value, it had not decreased substantially. It also should 


be pointed out that X could be better predicted than 
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Yoo as seen from the Ro and R values, 

The variables useful for predicting Xiog were; 
hatchery source (Xegs ogy) average hatchability rate for 
brood (Xe) location of air exhausts in broiler house 
(Xo75 Ze). age specific mortality rate at 4 weeks (X39). 


type of feed-crumbles (X,._), vaccine source/type (X,-,,) 


EVE 169° 
litter type (Xa pq) season (Xogps Xoo6)s and change in 


population size CT oes With respect to the signs of the 


regression coefficients, the two hatchery sources out of 
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eight selected in this model were inversely related to the 
disease rate X409° This means that when the source of 
chicks for one of the study farms was from Keg Oars. 
the birds had a lower disease rate and would have a lower 
condemnation rate at the time of processing, Hatchery 
Sources were represented by indicator variables and 
therefore, only the presence/absence of specific 
hatcheries (out of the total of eight used in the study) 
were considered, This means that in equation 2, the 
presence of Xeg will decrease the disease rate X109 by 
2.913 units (%). Location of air exhausts in a broiler 
house were also handled as indicator variables. In this 
model, when the direction of air exhausts was to the 
north, northwest, it was positively correlated with the 


disease rate (X ); while if it was in a south, southeast 


109 
direction, it nad an inverse relationship with X409° 
Tnus, the location of exhaust fans apparently is 
Significant in affecting the morbidity or disease rate in 
broiler houses. Although the explanation for it could 
only be speculated here, in terms of its influence for 
spreading organisms within the broiler house, this may be 
an important structural and biological consideration in 
poultry housing. Tne other set of indicator variables 
selected in this model were seasonal factors, Xoo and 

X 503° The signs for these coefficients indiate that the 


fall or spring seasons are correlated with a lower disease 


rate. The average hatchability rate of the birds under 
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study (Xee) was inversely correlated witn X109° If the 
chicks at the time of hatching had a high hatchability 
rate, this was correlated with a lower morbidity rate. 
The explanation here is that, a high hatchability rate 
Signifies good hatchery management, good genetic stock, a 
nealthly and mature breeder flock that produces stronger 
baby chicks. In terms of population related factors, the 
age specific mortality rate at 4 weeks (X39) had an 


inverse relationship with X while the overall 


109° 
population change Srey was positively correlated with 


X 109° when the flock size decreases due to mortalities, 


€.8- those birds may have died due to some disease 


X39" 
and if so, they did not have a chance to stay on and be 
part of the infected population, which X 109 tndirectly 
measures. Conversely, the rate of change in the 
population size, represented by Asa Cvizese ind taad 
population size during the first week, minus the 
population size at the end of the brood, divided by the 
initial population size), was positively correlated with 
Xing: This was consistent in the sense that it accounts 


for the change in population health which also 


i109 
estimates. 

In terms of the comparative importance of the 
variables used in equation 2, above, an examination of the 
Standardized regression coefficients indicated the 
th 


following descending rank: x 


117° 26° 


Again, the 


Xoqe pop’ X65? Xogs 


X and finally X 


z 169° *39 189° 


X 


X508° %o61° Xoo: 
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importance of air exhaust direction, rate of population 
decline and the average hatchability rate as the top three 
preaictors of disease rate in poultry population are 
Significant points to remember, 

The other vital statistics pertaining to the 
predictor model were the R, Re andor ‘ratio; “all sof which 
have been given with the predictor equation, The R for 
predicting X109 was very high (0.93), and using the 12 
variables selected, 86% of the variability of X 109 could 
be explained. This also meant that the model was quite 
satisfactory in providing reliable values for X 109° which 
could then be used in predictor equation 1, 

In addition to predictor equations for t, and X09? 
ares other important equations which dealt with disease 
specific conaemnation rates for septicemia, airsacculitis 
and the "other" disease category were evaluated. These 
three disease specific rates form more than 80% of the 
condemnations in processing plants. In this regard all 
the steps described in the computer analysis of the data 
were followed in a comparable manner as far Y5> The 
resultant equations obtained at the termination of the 


bakeward stepping procedure were: 


A 
49 = 0, 11609-01467 2X, .40.07269% 5, .+0 .00363X,, 
=e ‘ 5 Y = Q 
+0,03661K, 4 0.11383Y.,, +0.83637 10T 0.09161X540 
(4) 
2 
GRaRO sap eh = 20860 bP cratio = 10829) 
Y = 5 =O? U 
Ya, = 12.81016-0. 13026%, ,+15. 138%, 4, 127032, 


A 
Tay 


X49 
The pred 


~0.16964X 56 (5) 
Gye tich) Ge = ieee etal Sag 
= 0.04641+0.03736X,,~0,.09023869X 5.40. 11217X, 
0 ..00062626X 4 9+0.000694% 11, (6) 


GS Wea FR 6 GGT, eae 6 Oe Ba 


= condemnation rate (4%) due to Septicemia 

= condemnation rate (%) due to airsacculitis 

= condemnation rate (%) due to "other" diseases 

= hatchery source no. 2 

= hatchery source no. 7 

= number of ventrillators 

= disease rate (%) based on necropsy 

= condemnation rate (%) at one previous brood 

= condemnation rate (%) due to septicemia at one 
previous brood 

= number of days involved in marketing 

= average hatchability rate (%) for hatchery 

= disease specific rate (%) based on necropsy 
(airsacculitis) 

= number of rainy days 

= total precipitation 

= location of air exhausts in hatchery (NNW) 

= date of placing of eggs 

= strain group .1 

= average distance between houses on a farm 


ictor equations for these three rates were not as 
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Boca as for i Although the disease rate, due to the 


"others" category, (Yo) could be predicted 
Satisfactorily. 74% of the variability in Yay could be 
explained by the variables given in equation 6, Fach of 
the variables used in the above equations could be 
identified using appendix 9. The signs of the regression 
coefficients, their magnitude or importance, as well as 
the significance of the R, R° and F ratios are comparable 
to what has been presented above under the discussion for 
Y, and Xi 09° Although the Ro values for Yo and Yay 


(septicemia and airsacculitis) were not as high as was 


hoped for, it still appears satisfacory for providing 


decision-making information. It is important to point out 


that, it is better to predict the overall average 

condemnation rate (Y,) rather than the disease specific 

rates. 

Case 2: Predictor models which could be generalized: In 
an attempt to develop predictor models for Yo5 X09 
and condemnations due to disease specific rates 
(Yo-Y 44), variables which appeared to be of local 
Significance were omitted from further analysis. 
Included in this category were hatchery sources, 
vaccine sources and location of air intakes or 
exhausts. Additionally, Xone the variable 


representing the number of inspectors working in a 


processing plant was also omitted since its predictive 


vlaue may be controversial. The data analysis was 


ry 


performed in the manner described under methodology. 


The predictor models derived in this section were: 


¥ 


Whereé 
& 


Y, = 


This 


19. 3583-0.21023X, 440, 10855X, 4+0.20511K 


14 


g Oe N357 1X55, -0-00702X553* 
2 


(R = 0.91, R° = 0,82, F ratio = 19.07) 
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9.3227X, 5 ‘0. 18356Y5, (7) 


predicted condemnation rate (4%) 

average hatchability rate (%) for hatenery 
total number of strains on farms 

disease rate (%) based on necropsy 

water spacing per bird (inches) 

total precipitation (in) 


distance from farm to plant (miles) 


= condemnation rate (%) at one previous brood 


multiple correlation coefficient 


multiple coefficient of determination 


equation was one which was obtained at the end of 


the forward stepping regression analysis program. It 


utilized 7 predictor variables which consisted of 


average hatchability rate for hatchery (Xe3), the 


total number of genetic strains of birds on the farm 


(X14), the disease rate variable (Xi 99) water spacing 


per bird (Xj) total precipitation during the brood 


eycie (X5o5)s the distance from the poultry farm to 


the processing plant in miles (X 


223) and the average 


overall condemnation rate at one previous brood cycle. 


For the variables postively correlated with Y 


>" the 
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id, 


finding is consistent with the underlying biological 
expectation that one may have. For example, if the 
previous condemnation rate (Yor) was high, one may 
expect the current condemnation rate or the response 
variable (Y,) to also be high. 

In the case of the positive relationship between 
Xi, and Y. which indicate that if the number of 
genetic strains in a poultry farm are large, this 
leads to high condemnation rates. This of course 
could be a reflection of limitations in the managment 
or that the chances of getting poor doers from among a 
larger genetic pool is greater than that from a few 
which would have been selected intentionally for good 
performance. A similar management related argument 
could be made for X 409° 

The variables with an inverse relationship, X63" 
X505 and X 593 are not so straightforward to explain 
except for the first case. If the hatchability rate 
(X63) is high, this certainly could be indicative of 
good management and therefore the returns in lower 
condemnation rates at processing for the firm may be 
expected. Tne model indicates that the total 
precipitation during the brood cycle (X555) is 
inversely related to condemnation rate. Although this 
could be questionable, one reason here is that 
significant precipitation is observed during the 


spring or early summer months in Alabama which may be 
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the most conducive time for raising poultry with less 


problems of infectious diseases. The variable for 
distance from the farm (X53) is somewhat oeculiar in 
the sense that one would in fact expect to see a 
positive influence of X53 on Y>5 1.e. as the distance 
increases, the birds are exposed to more stress due to 


travel which may then exascerbate existing 


abnormalities leading to higher condemnation rates, 
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However, the relationship in the model was inverse. 
Tnis can only be speculated here in that possibly due 
to the fact that transit stress is known to be 
harmful to animais in general the management may 
provide good care at the time of transporting of 
birds to the plants. For an examination of the 
Standardized regression coefficients, the following 


model is presented: 


A 
= ‘ v xX 
Y, Ssueihes + 0.121X,,, #0500 109 + 0.109% 1 49 


= x f 
-0.498X 06 0.204 223 * O.0145Y,, (&) 


2 
Che=eQ.0 eneeeroep2, f ratio = 19.07) 


Where the variable names are as described above under 
equation 7. 
Note that in terms of comparative importance, X409 ia 
the most valuable predictor followed by ons) X53) 
X53) Yor and finaily X,, and X199° 

For predictor equations 7 ana §, the R and R” 
values were 0.91 and°0.32 respectively. This is a 
very high correlation coefficient and the equation 
successfully explained 32% of the variability in Y,- 
VWiealsouhead ausienificant F ratio of 19.07 which 
Showed strong linearity between Y, and the 
independent variables. Table 17 provides other 
useful statistics such as the standard errors of the 
regression coefficients, the adjusted % squares and 


other such data. Based on these, confidence interval 


estimation and specific hypothesis testings could be 
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Table 17 


Regression Statistics For Equation 8 


Multiple R 0.9064 
Multiple R-Square 0.8215 
Adjusted R-Square 0.7784 
SID: Error of Est. 0.4538 


Analysis Of Variance 


Sum Of Squares DF Mean Square F Ratic 
Regression 27.493883 yi 3.927698 19.07 
Residual 289132733 an, 0.2059749 


Variables In Equation 


Variable Coefficient Std. Error Std. Reg. 
(Y-Intercept 19235430) of Coeff. Coeff. 
X63 =()o2 L026 0.0428 -0.480 
X14 0.10855 020922 OF121 
X109 0.20511 0.0357 0.506 
X129 0.32270 0.3197 0.109 
X905 ~0.43571 0.0885 -0.498 
X99 -0.00702 0.0034 ~0.204 


Yor 0.18386 0.1458 0.146 























conducted to further scrutinize the performance of 
Che model htathrs stagesthough, it} is*signif icant 
to note that this predictor equation for Ys was the 
best thus far. One reservation to be raised here is 
that with the rest of the predictor models, the 
selection of the equation was made at the end of the 
regression analysis. In this case, the equation was 
selected at the end of the forward stepping 
procedure, which meant that some of the variables 
entered into the equation could be removed during the 
backward stepping procedure allows that some of the 
variables entered into the equation be removed during 
the backward stepping phase of the computer analysis 


f necessary. That in fact was the case below where: 


i 
y O82X 0.2285 0.44S509X 
Y, — 1822799 = 0. 18993X63 + e SS DXaqq = {6) D) 9 205 
ws Q 9 
0.00858X,,,. (9) 
(R = 0.89, R° = 0.79, F ratio = 29.99) 
Where: 
A 
Y5 = predicted condemnation rate (4%) 
X53 = average hatchability rate (%) for hatchery 
X09 = disease rate (%) based on necropsy 
= ipitation.(im) 
Xoo5 = total precipitation (i 
X59 = distance from farm to plant (miles) 


Again this equation is quite good although its R and 
Ro values are slightly less than in equation 8. The 
F ratio shows much stronger linearity. In both 


cases, the number of predictor variables are of small 
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size ara therefore manageable and easy to explain. 


Predictor equitions Mforsiuns “six, caterories far 
disease specific rates (Y,-Y,,) were also evaliuaten., 
vy 


HO WeVer | Cyese yom met 


v 


SS She SiChory 15 wag. cas 
case for V5. The next step was to develop a 
Ppredi itor “oged stare Yas + the Iisease rate hased ‘or 
necropsy. Tnis variable was an indicator of 


morbidity on a farm. The predictor equation of “aad 


for the seneralizable model wast 


“4ag = OD «SOO 7 Oat 9.735 9X en - 0.63653% 75+ 
NG LOY mie inne ae HINRWDY 4 
0.954 quant 0.5492h 190 eel TG? ann 
| ee + OF2509K 5 =e Teese 


Ceo Sab (i?) 


This moael aia not aonear as satisfactory, althouszh 


that the equation was not deened satisfactory for 

tnis wort:, detailed discussions of “the var zables has 

rot been presented. Instead, an attempt was made to 
tuprove the predictor equation via weighted least 

squares procedures as deseribedd under Case 3. 

felentea least squares: A tho stage least 
squares procedure was usec to accomplish this task. 
Firstly, using e@quationg@e 7) “he prealictea valueceta one 
Y's for each of the study cases were obtained. Secondly, 
the estimates for the wetents tobe used for each 


study case were computed] “ese estimated sverguma 


were then used to obtain the weighted least squares 
regression model. It should be noted that this 
d-scussion applies to the generalized model. The 
weighted regression model was; 


A 
Yo = 4.8499 — 0.05 184X_., + 0.1508°X 4), + 0,05382Xx t 
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0. 14903K 159 


0. 10326Y57 — 0,34003X 


+ 0.03337% 45. - 0,00365X,,.,. + 


9 (11) 
(Rean.0503, Ro SOs Corer CON cell Oe.) 
This equation was obtained at the end of the forward 
stepping program of BMDP2R. The variables with 
positive correlations with Y5 were: the number of 
poultry genetic strains on the farm (Kade the 
estimated disease rate of flock (Xigq) water spacing 
per bird (Xing) the total precipitation for the 
brood (X59) and the condemnation rate for one 
previous brood (Yop). Those with a negative 
influence on Yo» were: the average fertility rate 
for the hatchery (Xe) s the distance between the farm 
ard the processing plant (X53) and type of strain-- 
(X 


). Exeept for X,, the interpretation of the 


3] Be 
variables selected above have been discussed earlier 
under Case 2. The weighted least squares model 
indicatea that the presence of strain type 1 birds is 
good in the sense that a fewer percent of these would 
be condemned at processing. 


The standardized regression coefficients of this 


model indicated that comparatively, Xg and Kay were 
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the most important variables followed by X09" X5931 


2 
4 and finally X The R and R- values 


at? *129* X60 205° 


were very good although the F ratio was not as 
Strong. With the aid of weighted least squares, the 
accuracy of the model as measured by ® and Ro had 
increased some but not very substantially. 

Therefore, in predicting Y5 for a generalized model, 
either the weighted least squares approach or the one 
given under Case 2 were satisfactory and comparable, 
In a similar fashion, a weighted least squares model 
for X09 Was derived where: 


& 


= a e 
X09 =e Ten Vor 0.37163%, g + 


6292) 1 
0.36292K,.4, (12) 


(R = 0.999, R° = 0.0997, F ratio = 117.07) 


4 0.1091X,,, 


Tne predictor equations selected were average 
fertility rates for the hatchery supplying chicks to 
the farm (X63) the duration in days of prophylactic 
medications (X,,,,), and the average feed utilization 
rate at the end of seven weeks (X433)+ The model 
indicated that the longer the duration of pro- 
phylactic medication, the nigher the estimated value 
for disease rate in the flock. This is somewhat 
understandable in that if the flock of birds is kept 
in a healthy, well-managed environment, prophylactic 
medications may be used in a planned and judicious 
manner so that only the least needed time span will 


be utilized. If the flock is generally unhealthy, 
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the duration of prophylactic medications will be 
longer. Therefore, this variable (XayQ) does not 
indicate a biologically reasonable influence on X 09° 


In the case of feed utilization rate at 7 weeks 


(X33). the positive correlation with X,,,. may be 


09 
indicative of mismanagement of feed utilization, The 
most plausible explanation is that if the flock of 
birds is heaithy, they will tend to consume more and 
gain more weight in the process, If the birds are 
healthy, one would expect to see a lower disease rate 
and an inverse influence on X99 In this case, 
X133° the feed utilization rate was positively 


correlated to 09 and one possible explanation may 
involve a mis-management of feed. Although the 
records of feed delivered during the brood cycle may 
be large and therefore the feea utilization rate 
coula be high; the effective utilization of the feed 
by the birds may not be optimal. The explanation for 
the inverse relationship with fertility rate of the 
hatchery (Xp) as pointed out earlier, relies on the 
premise that a well managed hatchery will show a 
higher fertility rate and hatchability rate in its 
chicks. Tnese chicks either come from good genetic 
stock or the hatching process is well managed leading 
to lower health risks and lower disease rate later; 
Since the birds get a good start during the most 


susceptible periods in their lives. 
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In the weighed least squares model for Xa09 the 
R and Ro values increased from 0.77 and 0.69 
respectively to 9.99 and 0.99 for each of the two 
narameters. The F ratio increased from a weak linear 
vaiue of 5.38 to a strong F value of 117.07. In 
snort, the least squares model was the best predictor 


EOvae and therefore recommended as an effective 
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predictor equation. 


Discriminant models — to perform the analysis for 


developing the linear discriminant functions. to classify a 
flock into high or low condemnation groups at the time of 
processing, various decision making options were tried. 
Tnese included criteria for nigh or low condemnations set 
at 5%, 2%, 1%, 0.5% and 0.1%. These demarcations were set 
somewnat arbitrarily but witn the idea of providing FSIS 
with a variety of options for consideration. 

Tne cepenaent variable used to develop the 
discriminant models was Yo The analysis was performed 
using BMDP743; a stepwise discriminant analysis program 
useful for such tasks. The F~to-enter and F-to-remove 
values were set at 4.0 and 3.995 respectively. For the 
two categories of condemnation rates (high or low), prior 
probabilities were set at 0.5. 

For the dependent variable Yo. the set of 
classification functions to be used at each demarcation 
point for high or low condemnation categories, as well as, 


the correct classification probabilities will be presented 
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in sequence. Then, one generalized discussion with 
interpretations directea at one of the classification 
functions will follow. This should suffice to present the 
reager witn the basic concepts and applications required 
to understand and infer what needs to be done with the 
Other functions. 


a) Demarcation between high and low condemnation rate set 





at 5%: 

LC = ~ 2061.35986 + 46,30182K_, - 6.5999X, 4. - 
17.21999K 159 + 67. 26443X 556 (13) 

HC = — 1639.04882 + 41.10843 — 2.87306 - 64 .26469X 459 
+ 56, 15462X,. (14) 


Percent correct classification overall — 1060 
Percent correct classification for cases under 5% — 100 
Percent correct classification for eases over 5% -— 100 


b) Demarcation between high and low condemnation rate set 


at 2h: 
| LC = - 2061.35986 + 46,30182X_, - 6.5989Xi 99 ~ 
| TT.21999K 19 + 67.26443Xo9¢ (15) 
HC = - 1639,04883 + 41, 109N3X%,. - 2.87306 — 
64.26469X,,, + 56. 1SH62X 5. (15) 


Percent correct classification overall = 190 
Percent correct classification for under 2% = 100 
Percent correct classification for over 2% = 190 

c) Demarcation between nigh and low condemnation rate set 


otal 





LC = — 4091,15234 + 114.77233X%, + NT.HHO3X,, + 


Ne a es a a eS ASS MA Re 


d) 


e) 


1s) 


HC = 
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WY .89483X,. + 4Y.969B9X,, — 18.59022x4, + 


63 64 
TV QNT52X 1g + 43.139 1X 556 - 4 TITS 2Y op + 
34. 62482X, Gl 


- 3969.55249 + 107 67 14X, + HY, 3B4HU5X, + 


44. 18753%,, + HH, 3899X - 17 .56498X,) - 


63 


12.645 1NX, = GSE NEN Gace 2 


09 5 bad 
31,05351Xq (18) 


+: 45 60967 1X, 


Percent correct classification overall = 91.9 


Percent correct classification for under 1% = 93.1 


Percent correct classification for over 1% = 87.5 


LC 


ae, 
OQ 
{ 


Demarcation set at 0.5%: 


- 4863.71973 + 78 £5039 1X6 + 17 -2642X,, + 


28.68534X,, - 85.92729X 159 — 8.7370RXK aoe + 

iB 32762X 553 (19) 
4699. 33447 + 76. 82549K 4 + 75. 84322X,, + 

28. 30284X,, - 83.62469X 155 - 7-79638Xa95 + 
1,21027X55. (20) 


Percent correct classification overall = 94.6 


Percent correct classification for under 0.5% = 88.2 


Percent correct classification for over 0.5% = 100 


EG = 


iG © 


Demareation set at 0.2%: 


— 5.14413 + 0,11869X (21) 
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=" 1055803740. 0b2g24 (22) 
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Pereent correct classification overall = 78.4 


Percent correct classification for under 0.2% = 75.8 


Percent correct classification for over 0.2% = 100 


Demareation set at 0.1%: 























Le 


HC 


= 6.00129 = 208138X¢ 4 + 1.42193X 159 a 


0.1236'X 53 (23) 
— 4.25845 + U.61612K%,, + 7.37256X 159 ~ 
0, 00404X,.,2 (24) 


Percent correct classification overall = 97.3 


Percent correct classification for over 0.1% = 97,2 


Percent correct classification for under 0,1% = 100 
where: 
LC = tne equation for classifying a farm into low 


HC 


tt 


ft 


condemnation category 

tne equation for classifying a farm into high 
condemnation category 

average hatchability rate (%) for hatchery 
disease rate (%) based on necropsy 

water spacing per bird (in) 

total precipitation (in) 

ancubator type 1 

incubator type 2 

average fertility rate (%) for brood 

total number of strains on farm 
condemnation rate (%) at one previous brood 
strain group 2 

average fertility rate (4) for hatchery 


distance from farm to plant (miles) 


In each of the six discriminant models given above, a 


classification function to specify whether a flock should 


be categorized into a nigh condemnation (HC) or low 
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condemnation (LC) group are provided. These equations 
were coupled with the respective correct classification 
probabilities. 
Using tne classification functions provided, a flock 
is assigned to the group (HC or LC) with the largest value 
of the classification fimections. For example, for Case 1, 


the following values are given: Xe = 4 =e0 ls. 


129 
X593 = 12. From the LC equation where the demarcation for 

igh/low conaemnation is at 0.1%, Case 7 will have a 
Computed value of -6.88. For the HC equation, its 
computed value will be 2.67. Therefore, since the HC 
equation nad the largest value, the case is assigned to a 
high condemnation category. All other cases were handled 
in a similar manner. 

In this classification function, the probability of 
correctly classifying a flock into high condemnation 
category was 97.2%, for classifying into low condemnation 
probability of 97.3%. This indicated that the accuracy of 
the modei in providing aecision making information was 
extremely high and very reliable. All the classification 
functions given under the discriminant inodels had very 
nigh probabilities of correct classification, except for 
One case viz. the probability of correct classification 
wnen the demarcation was 0.2% between high and low 
condemnation rates. This was only in the higner 70's 
percentage range, and it was considered less satisfactory 


than the other 5 cases. 
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6.4. Selection of the final models 
6.4.1. Analysis of resiquals 
) To detect model aeficiencies in regression analysis, 
a simple and effective method was to examine the residuals 


Ziven by: 


) The Yi are the observed values of the dependent variable: 

| the Yi are the predicted values obtained by using the 

| given predictor equations. Tne difference between the 

| observed and the predicted (vi-Yi) provided the vaiues for 
| the respective residuals (43). Corresponding, to each el, 
in the computer regression analysis, standardized 


residuais (eis) were computed using: 


eis =.c¢i/s 


wnere s is the standard deviation of residuals. The 
Stangardized residuals nave zero mean and unit standard 
deviation. Witn a moderately simple sample, these 


residuals should be distributed approximately as 





jnaependent normal deviates. An appropriate graph of the 
residuals will often expose gross model violations if 
these were present. In general, when the model is 
correct, the standardized residuais tend to fall between 
2 ana -2 stanaard deviations and will be randomly 
distributea about the zero mean. Tne residual plots 


should show no distinet pattern of variation which would 


inaicate an inherent deficiency in the model. This may 
indicate inadequacies in the underlying assumptions or 
errors in the specified equations. 

Since an essential part of any regression analysis 
includes a careful examination of residuals to ensure that 
the assumptions of least squares theory were fuifilled, 
residuals of two of the predictor models (equations 7 and 
12) were scrutinized with the aid of computer plots 
(figures 16, 17). Three plots were utilized in each case. 
In the first one, the objective was to examine the 
distribution of ei about the mean zero. In both figures, 
16 andl7, a plot vs the predicted Yi indicated that 
indeed the residuals appeard to be randomly distributed 
around the zero mean in a balanced manner in about half of 
the cases i.e. about half of the cases were positive 
(above 0), the rest were negative (below zero). Note also 
that the residuals were all within the acceptable standard 
deviation value of + 2.0. Tnerefore, this test criteria 

A ns 
was satisfactory for the predictor models of Y. and 109° 
In the second plot (figures 1lé., 19), the squared value of 
the residuals Bae were plotted vs the predicted Yi. 
This was expected to show a clustering of most of the ei 
Close to zero, the graph did indicate it to be so. ‘When 
there are outliers, when the residual term is squared, 
these will be exaggerated and could be picked up from the 
plots. In these two predictor equations, since the 


resigual deviated within the acceptable range of + 2, 
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NORMAL PROBABILITY PLOT OF RESIDUALS 
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there were no major problems with outliers and in terms of 
moaei inadequancy. 

Tne third part of the residual analysis involved the 
examination of a normal probability plot of the residuals 
(figures18, 19). In this case, the residual (yi-¥2) was 
plotted against the expected normal deviate corresponding 
to its rank. Tne result of the plot was expected to 
indicate a linear pattern. Such, of course, was the case 


aA 


a 
for both equations to predict Y., or X, 


2 09° 
In this case then, with the aid of the three 


graphical analysis of the residuals, the predictor 


A & 
equations, for Y, and X were found to be satisfactory 


109 
and acceptable. Since the model was developed using 
"study" data set, the second data set referred to as 
"test" aata set was utilized to validate the model. There 
were 37 cases in the study data set and 17 cases in the 
test data. Using test data, only 2/17 of the cases were 
predicted unsatisfactorily, (standard deviation greater or 
lesser than 2.0). On the other hand if Sees wes predicted 
first using equation 12, and then its value used in 


A t ne 
equstscn fate, orediet Y., 5/55 cases were*predicted to 


a 
have a standard deviation greater or lesser than 2.0. dhe 
rest of the predictea values were with reasonable ranges. 
In a form of brief comparative statistics, the mean value 
nx eo 
for observed Y, Was 0.677%, while for the predicted Yo» it 


was 0.768%. Tnerefore overall, both models performed 


satisfactorily and the validation process confirmed this. 
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Tne only other follow-up could be validation of these 
modeis using inaependent data from other regions. 

One important consideration in performance and 
Validation of the model was an examination of the 
regression coefficients to see within what limits their 
parameter values would lie. Therefore, construction of 
confidence intervals for the coefficients were essential. 
The 95 percent family of confidence intervals (CI) between 
two numbers for each regression coefficient where the 
population parameters could reasonably be expected to lie 
were computed. These were: 


rs 
(a) For the predictor equation of Y.,; the confidence 


> 

intervals for the regression coefficients: 
863 = -0.21, -0.123 
B14 = -0.08, 0.297 
BIO9T=s0e seen Oe cio 


B129 = 0.331, 0.977 


B2058=50-019 weno i 
B223 =) 020145302018 


BY 10e 0 ues Oeu 2 


eT 
(b) For the predictor equation of Xiag the 
confidence intervais for the regression 


coefficients were: 


B53 


0.372, -0.155 


B148 = 0.091, 0.127 


B133 = 0.2812 OF4kS 
Ones Validation of Models 


In summary both the analysis of the residuals and the 











6.5 


validation of tne model using "test" data set indicated 
that the predictor models would be satisfactory. 
Additionally, the 95% family of confidence intervals have 
been provided for each equation to provide information on 
the range of values within which one would reasonably be 
confident to find the true or population parameter of the 
regression cofficients. 
Appiication of results to poultry inspection 

It should be noted that the objective in this study 
was to develop a predictor modei which could serve as an 
adjunct to the existing poultry inspection procedures and 
Standards. Therefore, the results presented so far are 
visualized in terms of whether a2 decision pertaining to 
the current inspection procedure is to be made in light of 
information provided by models provided herein. The 
decision to be made may involve either a "scaled down" 
inspection, or to maintain the status quo as made 
operational in the various processing plants. 

In order to effectively utilize the modeis developed 
in this project, two alternatives are presented. Repard- 


\ 


less of the alternative though, two prerequisites for the 
successful use of the model(s) are: 
- First, information for the variables identified in the 
equations selected above should be made available by 
the firm/farm at the time of processing or immediately 


prior to processing of the birds. These data should 


be collectea and presented on standardized data 
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collection forms. 

- Second, the veterinarian must take the data for the 
variables selected and generate a predicted value for 
X199 which is the aisease rate for the particular 
flock to be processed. 


~- Finally, using the X Value computed above, the 


109 


veterinarian calculates a predicted value of Yo» the 


condemation rate expected for the particular flock, 


The conputation of Xi09 and Yy could be handled in one 


of two ways. 


Alternative 1 —- Using hand calculators-Tnis approach is 


rather simple and straight forward. The veterinarian only 
requires preexisting information on the predictor 
variables. The coefficients of these X variables is 
known from the model and once their values are provided, 
one then computes either X109 or Y, directly using a hand 
caiculator. With the availability of various hand 
calculators, such computatioins should be accomplished 
With ease. 

Alternative 2 — Using microcomputers-The manipulation of 
data for scientific computations and/or for record keeping 
are nowadays handled easily with the aid of micro- 
computers. Information processing which includes the 
Storage, analysis and display of information needed for 
decision making by individual, or organizations such as in 
poultry carcass inspection may have to now rely on the use 


of computers. In the case of poultry, it is obvious that 
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a large amount of data is being generated by the poultry 
inaustry and the Food Safety and Inspection Service 
(FSIS). Most of these data may best be stored and 
retrieved for various decision making tasks. In the case 
at nand, a predictor model has been developed which relies 
on the availability of data for selected variables, To 
facilitate this, a program could be prepared so as to 
enable: 

- the coding and storage of relevant poultry 

related data on a computer, 
- the analysis of the data including computations 


ee 


of predicted Y. values using the predictor models 


2 

developed. 
Based on current information, the veterinarian could then 
easily interact with a microcomputer and obtain a set of 
instructions to facilitate the decision making task for 
the type of inspection best suited to individual flocks. 
Such a computer program could be prepared in the form of a 
menu where one peruses available alternatives and responds 
to those questions until the terminal objective has been 
achieved. 

Additionally, via such microcomputers located in the 
various poultry processing plants throughout the country, 
FSIS could maintain a network of information ge2thering 
devices with connections to a central data base in some 
specified place. Since such a data base is inherently 


rich ana quite accurate and reliable tyne case OF 


poultry, FSIS may be in a position, not only to use micro= 
computers for implementing the modeis developed in this 
project in the field, but it would also be in an enviable 
position of utilizing its data base for various tasks with 
potential benefit to the poultry industry as well as to 
FSIS. Therefore, it is felt that one of the best 
aiternatives for implementing the use of the predictor 
modeis developed here in the field could be via the 
introduction of microcomputers to facilitate data handling 
ana decision making in poultry inspection. This would be 
of great benefit to the poultry industry, and to the 
veterinary profession itself, whereby, for example, some 
of the expertise of veterinarians and poultry inspectors 
may be diverted to expand the utilization of current 
information technology in their respective areas. 
RECOMMENDATIONS 
Which model to use for prediction: 

For predictive purposes, the best predictor model to 
use for predicting Yo: 


A 
Y5 = 79,3453 - 0.21023X_ 0. 10355X.y + 0.2051 1X ing + 


3 
0.3227X.., - 0.4357 1%, 


“= TO02k 5OY 
0.00702 992 + 0. 18366 or 


9 5 23 


For predicting disease specific condemnation rate 


(Yo): 

i, = 0.21854 + 0.028 16%, , - 9.090534% 9 + 0. 13305%, = 
0,.0090606X, , + 0.0274 EK. + 0,00451% 149 - 
0.00153X 4119 - 0, 01245X a6 - 9.0394 2X0 


For predicting the disease rate for the flock (X99)? 
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(SSF 


A 
Xi09 = 27.61194 —- 0.37163X¢2 + 0.1091% 44g + 0.36292X132 
or 
A 
= L i ; “7% ns sQNCY oe norc|; 
Xigg = 39-S4751 — 2.9:267% 4 ~ 2.02945%,., - O.NN252X, 
e234TIXn =— 0455.159X 9065X 574 13% 
+ 3.2347X,, O.55159X35 + 1.9065K 19 + 1.5741 169 
- af - 5 Hy! 
+ 1,0H076X4 94 1,8702X.40 1,.925532¢ + 9-3ON82T A, 


The model for classification purposes: 


For classifying a given flock into high or low 
condemnation groups, the best discriminant models for 
arbitrarily set demarcation points for condemnation rates 
have been provided under discriminant models. Those 
equations were very good in provicing decision making 
information at the indicated levels which varied from 0.1% 
to 5%. 

How to use the models: 
Tnis had beeen described above under the various 
sections. However, in surmary form, in order to use the 
models developed rerein, the steps required are: 
~ collect data for the selected variables identified by 
the model using standardized data collection forms. 
This way, each farm would be collecting and providing 
FSIS with records which would be comparable and 
hopefully reliable. Such records should be made 
available to the inspecter in charge before the birds 
from a given farm were to be processed. 

- using these records, the inspector in charge will first 
compute a predicted value for the estimated disease rate 


} in the piven Ghee he SO Malt 
(Li gq? in the piven flock o 
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- then, the value for X09 and the other selected variable 


data are utilized to compute a predicted value for Yo or 
Yay the average condemnation rate due to diseases or 
the disease specific conaemnation rate referred to as 
the others category. In a similar manner, vaiues for 
low condemnation rate (LC) or high condemnation rate 
(HC) are computed to facilitate the classification of 
flocks into nigh or low categories. 

~ based on the value of Y5. LC or HC, a decision will. be 
made by the inspector using guidelines set up by FSIS as 
to where to draw the line between a (predicted) 
condemnation rate which may require “scaled down" 


inspection vs another alternative. 


Areas for further studies in the future: 


The work involved in this project, overwhelming as it 
Was, haS opened up numerous areas for which further 
studies are required. The results from such studies will 
be useful in future decision making in poultry inspection, 
A few of the major areas for such a follow up are briefly 


presented below. 


Independent validation of the model developed in this 


project in different regions of the U.5.:; 


To test and verify the predictive value of the model, 
it is essential that it be utilized in the field in 
different regions of the U.S. on a pliot study basis? 

Only then could its validity be established before wider 


applications as an adjunct to existing poultry inspection 
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procedures are carried out. The pilot validation test 
implementation could involve randomly selecting a few 
processing plants in the country and then testing the 
results under those situations over a period of a few 
months to a year, Based on the performance of the model 
under these circumstances, a final decision could be made 


on its utility in future endeavours. 


Conducting comparable types of studies, as conducted in 
this project, in other regions of the U.S. - This is 


another option which could be beneficial in the long run, 
This point becomes significant in the sense that, rather 
than one general model which could be applied uniformly 
across the country, the preferable option may be to 
develop modeis which could be useful on a regional basis. 
Additionally, due to passage of time, the diseases of 
importance in current poultry condemnation or their 
determinants, whether on a national or regional basis, may 
change and therefore, the predictive models in use will 
nave to be reevaluated and adjusted to accomodate such 
changes. 


Condemnation of parts of poultry carcass due to diseases - 





This aspect of poultry condemnation, no doubt, will be an 
area which urgently needs further studies. Currently, it 
is not required by FSIS, to recora separately from other 
causes, the parts condemned aque to diseases. Yet, a size- 
able mount of inspection time is spent on determining 


which part(s) are to be trimmea by plant personnel. Poth 


(rer 


(cheer 
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of these tasks are not only time consuming, but also 
complex. There are numerous disease conditions that 
require trimming. To name a few, gangrenous dermatitis, 
hemorrnagic syndrome, scabby nip syndrome, staphylococci 
infection, and fowl pox. When whole bird condemnation due 
to Gisease decreases substantially, parts condemned may 
then become the important contributor to economic loss, 
Therefore, accurately recording the trimming due to 
different reasons will benefit the processors to correct 
these problems. 

Revision of recording method and terminology. The MPT 
form 514-1 (Poultry Condemnation Certificate) should be 
updated. For example, tuberculosis is no longer a common 
disease encountered in poultry processing plants. 
Tnerefore, this item should he deleted from the form. 
Numerous disease conditions which result in condemnation 
of parts of the bird are not recorded at the present time 


and shoula be aaded to the form. 


Serological profile studies — Although the value of 


serological profile as a predictor for condemnation has 
not been snown, there is a need to conduct more studies 
before its importance can be established. There are many 
factors affecting serological profiles. Tne limited 
resources and time spent in serological tests in this 
project, by no means were adequate to examine the validity 
of this factor. There are several disease agents proven 


to be directly or indirectly associated with 











eects O, 


airsacculitis. They are MS, "G, “Newcastle disease virus, 
anfectious bronchitis, infectious bursal disease and 
pernaps adenovirus, Among these agents, commercial flocks 
are routinely vaccinated with Newcastle disease and 
infectious bronchitis virus. An average vaccination titer 
of normal birds at a specific age should be predetermined 
before a diagnosis of disease condition can be made, 
(although in a disease condition the serum titer is 
usually higher). In order to thoroughly study the 
Serologicai profile of a flock, antibodies against 4G, MS, 
Newcastie disease, infectious bronchitis and infectious 
bursal disease should be monitored throughout the life of 
the flock, instead of testing only at 4 and 7 weeks of 
age. An increase in HI titer for Newcastle disease and 
virus neutralization titer in the sera of broilers about 
10 days before processing, may indicate a condition of 
active infection in the last few days, as a result of the 
infection, condemnation due to airsacculitis may increase 
since airsacculitis is believed to be related to ND and 
IB, However, the relationship of serologic titer wltn the 
rate of condemnation nas to be determined before one can 
Speculate its importance in predicting the rate of 
condemnation. To the investigators of this project, this 
factor (serologic titer) may still nave potential in 
preaicing condemnation. 


“Microbiological profile - Bacterial and fimgal flora in 


a 


c 


hatcheries and poultry houses were stidied in very few 
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flocks. Since the presence of large numbers of E. coli, 
Pasteurella multocida and Asnergillus fumigatus will 
influence the incidence of airsacculitis, more studies 
concerning the bacterial ana fungal flora in their 
environment are needed before definitive conclusions can 


be made. 


Specific guidelines for condemnation of septicemia and 


toxemia are needed = Septicemia and toxemia are conditions 


developed in animals as a result of infections caused by 
different agents. The ciinical manifestations are 
different, and guidelines for diagnosis are poorly defined. 

Birds condemned for septicemia-toxemia have become 
the single biggest category of economic loss for poultry 
companies. Apparently, due to ambiguous criteria for this 
disease condition, sometimes instead of condemning part of 
the bird, tne entire carcass had to be condemned, simply 
because the bird "looked septicemic". Therefore more 
clear cut guidelines for this condition is urgently 
needed. However, before the new guidelines can be 
established, research concerning the true picture of 
septicemia—toxemia in terms of signs and symptoms, and 
whether the meat from birds showing signs of septicemia 
and toxemia were detrimental to animal and human health, 
must be established. may be possible to salvage part of or 
the entire bird, provided that clear cut guidelines are 
established. 


In recent years, condemnation due to this category 

















(ees 


has increased, the increase may have been the results of 
mild forms of disease conditions such as ND and IR, or 
adverse effects of vaccination. There are several 
emerging disease conditions wnich have been reported in 
recent years, e.g. malabsorption syndrome, pale bird 
Syndrome, brittle bone disease, proventriculus adenovirus 
infections. In addition, infectious bursal disease has 
become wide spread in all parts of this country; and early 
infection of this disease in baby chicks resulted in 
immmosuppression, and nany emaciated birds. All these 
disease conditions could result in increased condemnation 
rate in the septicemia and toxemia category. Therefore, 
more research is needed in this area in tne near future, 
Others — It was consistently demonstreated in this study 
tnat there is a close relationship between weather, 
ventiiation, strain of birds, number of strains of birds, 
medication records, different vaccines used and 
condemnation rate. It is the authors opinion that long 
term studies concerning tnese aspects of poultry health 
Should be initiated. 

The relationship of feed consumption and feed 
conversion to condemnation has not been established, due 
to the fact that there is no accurate method for measuring 
the weekly or daily feed consumption. However, accurate 
measures of feed consumption during a given period will be 


neipful for future analysis. 


209 


[ae 


Concluding Remarks 
Finally, it is significant to reemphasize here that 


analytic models wnich interrelate such biological and 
mathematical complexities are relatively new in veterinary 
medicine and are now receiving long overdue attention (68, 
69, 75, 76). It is encouraging and pleasing that FSIS had 
the foresight to consider such an alternative to poultry 
inspection wnich provided us an opportunity to explore 
extensively, the utility of such hoiistiec approaches to 
problem solving in an industry where population health 
data is relatively valid and well maintained. 

The future for the application of systems analysis 
and the development of computer models to facilitate 
decision-making tasks in veterinary medicine appears 
bright; this project was one such example. At a time when 
information processing is erucial to decision-making, more 
areas of veterinary health care would have to be managed 
With the aia of computers. The inspection of poultry 
products appears to be one such area where an introduction 
of computer asided decision-making models could be of 
great value. Although, this project was but a beginning 
in such a direction, further research and testing of 
results will be essential before the full benefits of such 


probiem solving tools could be realized. 
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APPENDIX 1 Form Code 
POULTRY DISEASE RESEARCH 
School of Veterinary Medicine 


Tuskegee Institute 


Flock Testing Feasibility Study 














4. Breeder source 


A al Ee 


CODE/COL. 
_ Cross-Sectional Sanpling. - Farm Data 
NOTE:[ The farm data should be followed up by corresponding data 
PrOvsprocesciDceplLalt. awe wee Pe 
Meee ema /dy/ vt) kre sabe Visit No. (1.2,3,4) 0...» 
_ Data collected by (Code) 
_ County 
Firm Code No, of growers of this firm  -—s§s— ss s—s—sdg¥|____ 
frre Colcmmennns WP eeeBroodsNows(15255,4555)—— — = 
Distance to nearest poultry farm (miles) __ ics an, ate eee 
A. Population Data For Farm: 
1. Total number of birds (lst week) Tony) ae td. wien Sa | 3 es 
2. Age of birds (days) K4 
3. Average weight of birds at this age (lbs) 35 : 
4. Are birds of this age group of uniform size? X6 
l yes 
0 no 
5. If not, estimate proportion of .undersized birds (4%) a7 
| 6. Major strain of birds: 
1. Strain group l hae Oe ae BNE Ns oe 
| 2. Strain group 2 Anime es ee. ee eee 
B “ ceo 0 
SOL ali 2rOup 1 eee ee ers | eps 
| 4, Strain group 4 pes 5 nit 3 Cle] te eReader 
| 
; 7. Breeder Flock sources 
1. Breeder source 1 ae ees SS Sages. 
2. Breeder source ? rivet) ee eee xe 
3. Breeder source 2 E a lg a X15 
4 22 


8. Strain/breeder follow-up: 
House Id No. of Birds Strain Type Breeder Flock 


ome 











Totals — SS 


Page 2 


hea 


14. 
Le 


16% 
ie 


LB. 


Totals 


Total number of strains on farm _ 
Total number of breeder flock sources for farm _ 
Housing: 

1. House type 1 (conventional/fan) 

2. House type 2 (conventional/no fan) 

Si aHOUSes LYPens a. OLUCLS) aaa 


a a mn sae 


House type follow-up: 


House Id Type of House Area of House Total Birds in 
(sq. ft.) House 


oe a Ser 


Average population density per square feet 


Average distance between houses on farm (ft.) 


Type of ventilators: 


Le elype ly G CEH) 
2.) alype 2 (> CEM) 
3. Type 3 (Others) 


Number of ventilators (Total for all houses) 


Location of intakes: (geographic direction) 


ine Stefehe 


2. (hast. 2noGi he eactaus els et bale) Se eee 
3. SNorth, norch westaa | 
» West, south west 


5. South, south east 


Location of exhausts: (geogrpahic direction) 
Lat ROOL Rae eee 
2. sbast, northeast 
3. North, northwwest 
4, West, south west 
5. South, south east 
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Form Code 


_CODE/COL. 
X14. 
X16 


X16 
X17 
Z3 





__ X18 


X19 


CAS 
X26 
X27 
X28 
Z6 


B. 








**k Denotes values computed. 


Morbidity/Mortality data from farm records 


lise 


6. 


By housing: No. 


Totals 


Farm morbidity rate (prevalence) 


Farm mortality rate 


House type specific mortality rates: 


House type 1 _ 


House type 
House type 


House type 


House type specific morbitity rates: 


House type 
House type 
House type 
House type 


By age (week) Estimated % sick No. Of ebirds 





ne ae 


lst 
2nd 
3rd 
4th 
5th 
6th 
7th 


of birds sick 
House Id. so far 


2 


WT 


ih 


1 


Z 
3 
” 





No. of birds 
dead so far_ 


LotaieNOrs 
birds in 
house 
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Form Code 


CODE/COL. 


en ee eeerene fone 


Total No. 


Gtebirds 
on farm _ 


220 














Page 4 Form Code 
= CODE/COL. 
7. Age specific morbidity rates: 
1. end of week 4 _ ck es. X37 
2. end of week / 78 X38 
Kk 
8. Age specific mortality rates: 
lo end sol tweekta tee re eee ro teeeene aan e et ate ee X39 
228 endsoreweelm/s Chive cal X40 
9. Major causes of morbidity (diseases) on farm based on 
farm records (top 5 diseases at end of 7 weeks): 
House Id. Diseases No. affected List of Diseases 
fee infection)” © Oe ee eee eee 
i Leukosis =e 
2, Septicemia 
5. Alrsacculitis 
Ter ee RE rc We Aids NONE woth ES CRS OVA Sateen 
5. Tumors 
6. Others 
16 Disease specific morbidity rates for farm: 
4 es X41 
ae eS oA. et | ee X42 as 
a Areas nets ES ep oeaaS 
4, << a etn Ns X44 
ee ee. <2 Rs I nc cmnery io aertee | an X45 
6. 4 8 SS eee feet rhea RAPE 4 
11. Estimated (recorded) incidence of diseases on farm 
(2). Saar ereniet bee X47 
12. Estimated (recorded) prevalence of diseases on farm 
(Z) X48 


% 


* Denotes values computed. 


aa 


Form Code 

















‘ CODE/COL. 
13, Major causes of mortality on farm based on necropsy 
records of farm: 
House Id. Disease No. Affected |List of diseases 
(at Necropsy) 
1. Leukosis 
_ 2. Septicemia 
3. Airsacculitis 
4. Synovitis 
Jee LUMOLS 
(ens 6. Others —. 
re caacuara Cpa cadgead uaa 
14, Disease specific mortality rates from farm necropsy 
records (%): 
Tg ca Wea SS acta 
a X50 
3 . a X51 
4 X52 
5 X53 
6 X54 
ee eene snc ener ren See se ee 
15. Estimated (recorded) mortality rate on farm (4) 
(4 weeks) X55 
16. Estimated (recorded) mortality rate on farm (%) 
Cjaiceks) [eens Wea Sa 2 kOe 
CeeHatch Data: 
1. Hatchery source for this brood 
Te pHatchery Source 1 So $a tp tee RO] 
2. Hatchery Source 2 > X58 
3. Hatchery Source 3 ear weak’ Set ae X59 
4, Hatchery Source 4 Sie p X260 
Depmratchery isounce’ 2 , rrr X261 
6. Hatchery Source 6 X262 
7. VHatchery Source, / X263 
8. Hatchery Source 8 Ll 
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ity 


ce 
eo 
14. 
1% 
16. 


Type of incubators: 


Incubator type l 
Incubator type 2 
Incubator type 3 


Distance from hatchery to farm 
Average fertility rate for hatchery 


Week Breeder Flock Fertility Rate Hatchability Rate 
(Age) (%) (%) 


oo 





Average hatchability rate for hatchery (%) 
Average fertility rate for this brood of birds (4%) 
Average hatchability rate for this brood of birds (%) _ 


Number of days eggs held on farm 
Number of days eggs held in hatchery 


Are eggs from more than one flock hatched simutaneously 
in same hatchery? 


l yes 
0 no 


Number of breeder flocks supplying eggs for hatchery 


ee ee 


Average temperature (F) for holding eggs (farm) _ 


Average temperature (F) for holding in hatchery _ 


Average temperature (F) in hatchery 


Average relative humidity in hatchery (%) 


Sanitation assessment for hatchery area: (where ] = 
good, 2 = satisfactory, 3 = poor) 


Criteria Grader (Gaze) 
1. Offensive odor (unspecified source) 
2. Dusty room (and equipment) 
3. Cobwebs on walls/ ceilings 
4, Freedom from garbage and refuse eS Re. 
5. .Wast disposal containers availability 
6, Unwashed (unsanitary) equipment/areas) 


X60 
X61 


Z8 


X268 
X62 


X63 
X64 
X65 


X66 
X67 


X68 


X69 


X70 
X71 
X72 
X73 


ae ats 


Form Code 


CODE/COL. 
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Page 7 Form Code 
7 ——————--—-~-| CODE/COL. 
ih Hatchery sanitation index (average of grades for ak Ch EAR oc ae 
CEEcenial.0) ce X74 i 
18.. Type of ventilators: 
Om Cre ALO Ce LY Degas nCrM) 9 Ty Sen Al) Ye ah 7 RSE XT 
2. Ventilator type 2 (5 CFM) X76 
SvVentilator type 3" (Others) =" 9" Ee ed tea ae 
19<. Number of ventilators onstrate GU bes ss oe 
20. Location of intakes (direction): 
1. Roof — ees Sy X264 
Zee bast, anorth east ene Ih) 
3. North, north west a X79 a 
4. West, south west fans oan. oe OU ee 
5. South, south east . PE ee eT eh LLU ee 
21. Location of exhausts (direction): 
Te ROOf eR: X265% | 
Ze bast, vorch east AS 1 Stee 
3. North, north west at everwes £8). A ew 2X8 2% 
4. West, south west X83 
5. South, south east PAL 
D. Brooding Data: 
I. Date of setting of eggs for this brood (mo/dy/yr) X84 
2. Date of placing of chicks for this brood (mo/dy/yr) X85 
3, Number of chicks dead or arrival X86 
4. Length of time for preheating (brooder) (hrs) ss X8 (ee 
5. Brooding: ein 39 Sees 
DaaNnOLesS : 
Ome rartial 
6, Length of stay of birds in broiler house (days) sss], KBD 
JemeAce wunensbrooderseuards are removed (days) “| KIO 
8. Number of chicks per broiler house (average) ss iT CX 
Geestitares footage per, bird per broiler house — is). wt ee 5 h92 


10. Microclimate (Average temperature/humidity) in broiler 
house (if available) 


Week Average Temp. (F) Average Relative Humidity (%) 


Lec 
2nd 
ord 
4th 
5th 
6th 
7th 
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Page 8 Form Code 
oy ~ CODE/COL. 
ws 
ll. Average temperature (F) up to 4th week se X93 - 
Kx 
12. Average temperature (F) end of 7th week x94 
ke 
13. Average relative humidity up to 4th week (4) X95 
KK 
14. Average relative humidity end of 7th week (%) X96 

E, Serologic Profile: (Laboratory) 

l. Indicator Age of Birds No. of Samples No. of Pos. 
Disease Examined bad 
ND ___ 4 wks Ha 
MYC 4 wks Pate 
ND 7 wks 
MYC 7 wks ie av 
x ——— 
2. Average infection rate based on serological tests (%) 
Loisat 4 weeks __ MT og Sey X97 
2. at 7 weeks X98 
3. Commercial serological data: 
Date of Age Tested Serological Tests Total No. 
Testing (days) Performed (Type) Samples of 
(mo/dy/yr) Tested Pos. 
4. Average infection rate based on commercial tests (%) “00 X99 
F. Microbiological Profile: (at 4 weeks of age) 
1. Lecation Sampled Population estimate 
Hatchery Room X100 
Brooder House X101 
2. Number of microbial colonies counted: (Average) 
Hatchery Brooder House 
1, Media J (@iocod*base) ro? 


2. Media 2 (McKonkeys) 
3. Media 3 Cornmeal) 


** Denotes values computed. 
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4. Major types of genera identified: 
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Form Code 


SCODE/COL ee 










Hatcher Brooder Ho 
1. Genus ee ——ee X104 
2. Genus XLO5 
3. Genus _ peooee: Be: X106 
4, Genus Tat X107 
5. Genus atttet vod ly hih vee See XL035 7 
G. Necropsy Profile (Laboratory): 
1. Age at Necropsy Total No. of Birds Total No. of Birds 
Necropsied ____—_s—s«éDiseased (grossly) 
sm niatlatl Se 5) a a th ot 
ate ALES elas iam 
KK cere 
2. Average disease rate based on laboratory (necropsy) (%) 
1. at 4 weeks X108 
2. at 7 weeks X109 
3. Top 5 common diseases diagnosed (Laboratory): 
Total No. Birds Diagnoses No. Diagnosed |List of 
Necropsied (disease) Diseases” = 
1. Leukosis _ 
ee es nanos wie 2 tsephicemtas 
«| 3s _Airsacculitts 
| Lae oe) St he synoviae 
Ree ee oe Tumors 
Pere ey Cee 
x 2m Venda oc ph 
4. Disease specific rates based on necropsy (laboratory): 
1s RT oka eae 
ve Be nt X1il 
3, X112 
vo ee eet ence ee a Saws 
lee X113 
nn bere ee 8 ee beet AS 
5 X114 
4 OOO a pk Sil ete is ss Se 
6. -X115 


**k Denotes values computed. 
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Page 1l Form Code 
CODE/COL. 
15) Average feed utilization rate per bird for farm. (lbs./ 
bird) 
1. First 4 weeks Ki32 
2. End of 7 weeks _ X133 


a ee er 


wk 
16. Rate of (average) weight gain per (average amount lbs.) 
of feed consumed per bird for farm: (lbs./bird) 
1. First 4 weeks X134 


2. End of 7 weeks X135 





17. Information on feed analysis (if available): 


Contents Microbiological Chemical (Toxicological) 


18.. Microbiological contamination: X136 


l present 
O absent Aes 

19. Chemical (toxicological) contamination: ZL37 
L Spresenct 


QO absent 


I. Therapy/Prophylactic Drugs Data: 
1. Prophylactic drugs used: 


Week No. cf Birds “No. “of Birds) Noy of bircesleaa 
Medicated Became Sick 





ist 
2nd 
3rd 
4th 
5th 
6th 


7th 


*x*sDenotes values computed. 
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H. Feed Consuption/Utilization Data: 
l. Type of feed used 


de aeiash os 
Phe Crumbles _ r: 
3; -. Pellets 


2. Feeding schedule (times per day) 
leaad lib 


Oaeat intervals _ 


3. Feeder space per bird 


——— 


4. Location of feeders during brocding (height in in.) _ 


5. Feed spillage on ground: 
Les 


Q No 





a 


6. Feed source (company): (is same as firm) A,B,D 


Pee sources L. ( ) 
Zou oource, 2. ( 
SMES OUTCe TG aC 
4. Source 4 ( 


Spee pouLces a 


— we WY YS 


7. How often is feed delivered? (per month) 


8. Last date of feed delivery (mo/dy/yr) 
9. Number of waterers used during brooding 


10. Type of waterers 
1 manual 


Q automatic _ any! 
li. Water spacing per bird _ 
12. Age manual waterers removed 


13. Water spillage on ground 
1 Yes 
OTN 

14. Feed utilization by farm: 


House Id. No. of Birds Feed Consumed/ 


Week/house __ 


Estimated Weight 
Gain/visit 


area | ore 
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Estimated efficiency rate for prophylactics (4): 
1. End of 4 weeks 


2.7 End of -/oweeks 


Route of prophylactic medications: 


l via feed 
2 via water 
3 via aerosol 


4 others 
Type of prophylactic medications: 
1 antibiotics 


Z santicoccidials 


De OLCloonaae 
Types of specific diseases prophylactics applied for: 
1 Bacterial pay 
Ze Viral 
3° Parasitie 
4 Others 


Average duration (length of prophylactic medications) 
(days) 


Dosages of prophylactic medications are based upon: 


1. per ton of feed 


2. number of birds 


3. other guidelines 


Therapeutic drug use: 


Week No. of Birds 
Sick 


No. of Birds 
Treated 


No. of Birds 
Recovered __ 


No. Birds 


Nksit: 
2nd 


3rd 
4th 
Bye 4) 
6th 
7th 





Dead ae. 
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Form Code 


CODE/COL. 


X140 
X141 
we ALA 
214 


X143 
X144 
ZS 


J 8 Ae 


X146 


X147 
216 


X148 


oa 
cKO 
Z17 


ee eee ere ern seraenneripnenerneanenies 
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10. 


11. 


fifa 


Ske! 


14. 
15. 


16. 
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Form Code 
ee Arey CODE/COL. 
Estimated efficiency rate for therapeutics (%) 
iS End of 4 weeks XL51 
2. End of 7 weeks X152 
Types of diseases therapy applied for: 
l Bacterial X153 
2 Viral X1L54 
4) madaewysive, |) Me ites X155 
feOCH erst As ees No code 
Types of therapeutic drugs need: a " ey 
1 Antibiotics X156 
2 Sulfas X157 
ou Others) © No code 
Routes of therapy: 
1 via feed ean? Mh, . Sherr £ X158 
2 via water R59 
3 via aerosol X160 
4 Others No code 
Dosages based upon: 
Number of birds X161 
2 Pewiton. of, reed cil Ma Rae X162 
3 Others rer e No code 
Average length of therapy (days) X163 
Feed additives and growth promoters: 
UrecComsdd ttt Ves LOL yoo Foe ee ee X164 
2 growth promoters only __ X165 
3,9no additives ZL0 
Source company for additives and growth promoters 
LeeCOmpanye. X166 
JeaConpetye. 2 X167 
3 Company X168 
4 Others No code 


Se ee ee eae 
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Page 14 Form Code 
J. Vaccination Data: = CODE/COR ae 
l. Manufacturer(s) of vaccine : 
Vaccine Manufactuter Names of Manufactuters 
1. Vineland ms i XI6g 
2. Sterwin : X170 
3. Salisbury 4 XE 
4. Others r Z19 
2. Routes of vaccinations: 

17 Oral en tan eye 
injections X173 
others Z20 

3. Schedule of vaccination: 
Time Vaccine Type (l=live virus,|List of Vaccines 
(Week) given (used) 2=attenuated, 3= 
killed product) on 
lst ieeMarekssiDise 
2nd 2. Newcastle 
3rd 3. Bronchitis-Inf 
- Gumboro (inf. 
4th bursal 
Sth aA Pie es 
6th —— SL os Oe | ONL Ea eee 
(Chee eat Nel Ste eee 
4, Vaccines given before 4 weeks of age: 

isu Vaceine (g/d | 

Zoe Vaccine 2 ee ae XS 

3. Vaccine a ee X176 

4 Vaccine ee 221 

5. Vaccines given after 4 weeks of age: 

Leeeace le on X177 

2.) vaccine tans X178 

3. Vaccine X179 

4. Vaccine No Code 

6. Number of repeated vaccinations X180 
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Page 15 Form Code 
7. Type of vaccine(s) given: CODE/COL. 
1. Killed product _ Wea ee X181 
2. Live virus X182 
3. Attenuated Opa 


rr 


8. Previous types of vaccines used on farm: 


1. same as above X183 


O not same as above 





9. Number of birds showing reaction to vaccine X184 


10. Severity of reaction: wees ; 


1 severe 


2 moderate 
3 slight 
K. Sanitation Indicators: 
1. Assessment of broiler farm sanitation (where grade of 
l=good, 2-satisfactory, 3=poor) 
Criteria PEAdeomure can) 
1. Offensive odor (unspecified source) 


2. Dusty room (and equipment) 


3. Cobwebs on wall/ceilings 


- Freedom from garbage and refuse 4 


ZL 
5. Waste disposal containers availability 
6 


« Unwashed (unsanitary) equipment/areas __ 


2. Farm sanitation index (average of grades rOL criteria 
1-6) bona A BOYS 


eS EN NE SC I EAST RS RTE RS SLO NO 


3. How are dead birds disposed of? 


fi hae — X187 
2 incinerated X188 
3 others Z23 
4. Type of litter used? X189 


1 shavings ieee 
2 peanut hull ie ae Ay Saree 
3 sawdust 224 
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L. 


DM. 


107 


Depth of litter (estimate) (inches) 


Howsotten isi theslitter chanses aa 


Is the litter changed: 
eeePactie yee 
Q completely 

How is the litter removed: 
f tractors 
Q others 


How many batches in this litter? 


Procedure for house cleaning after litter is 


1 remove litter 

2 wash floor 

3 disinfect floor 
4 not any of above 


Health Care Data: 


1 
Pe 


af: 


4. 


Cost of medication ($ per brood) 


Cost of vaccination ($ per brood) 
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Form Code 


CODE/COL. 


ee end 


removed? 


No 


Number of times a veterinarian visited farm to provide 


health care for this brood 


Total number of hours veterinarian spent on farm for 


this brood 


Macroclimatic Data: 


Le 


‘Le 


oe 
4, 


County weather station where farm located 





Summary macroclimatic data at (two week) intervals; 


Dates Avge max. Avg. temp. 


(mo/dy/yr) temp. (F) (F) 


Avg. min. 


Total Preci 


temp. (F) pitation (in) 


ee ee 


ne ee ee 


Se —- — OSSS——-,—C NNN ---=- = 


ee 


DS eee Eee 


Average maximum temperature (F) 


Average temperature (F) 


ee x0 


X192 
X192 


X193 


X194 


XI95 
X196 
X197 


code 


X198 
X199 


X200 


X201 
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Page 17 Form Code 
CODE/COL. 
5. Average minimum temperature Sok" eee eee ROG 
6. Average total precipitation ee SA 
7. Number of rainy days ZR 
8. Season during which this brood of birds was raised: 
he ew X206 
oom inter | anal. a! un 
3. spring X203— 
4. Summer a Z2IUR 4d 


Supplement: 


Management Indicators: (Below, a grade of l=good, 
2=satisfactory, and 3=poor) 


1. Assessment of farm record keeping: 


Criteria Grade (1, 2,3) 


1. Records well organized 


2. Records easily retrievable 


3. %5Z of information easily available 
(grade=1) 


4. 50%Z-75% of information easily available 
(grade=2) 


5. Less than 50% of information found (grade=3) 








2. Record keeping eS : ee Pe 
3. Are records kept in a secure place? pete AN ee as 
i yes ie etd, 
0 no 


4. Cooperative attitude of farm personnel: 


Criteria Grade (1,2,3,) 


1. Answer questions courteously 
2. Volunteer information 


3. Provide assistance when asked 


ana 


5. Index of cooperation (average) has; X211 





6. Assessment of farm (organization) and management: 
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HE Average grade for farm management 


8. 


10. 


ET. 


dias 


Criteria Grade (1, 2, 3) 
1. Grade for recording keeping 

2. Grade for cooperation 

3. Grade for sanitation 
4 


« Others (if any) 


Security ~ Is there a disease surveillance program on 
this farm? 


1 yes 


ORR aes 


If yes, explain briefly 


Security - Is there an emergency plan to deal with a 
disease outbreak on this farm: 


leaves 
QO no 


If yes, explain briefly 





Physical security of farm ~- Is there a program to prevent 
theft and other related problems 
Loevesan 


O no 





Are there safety and other warning signs for employees 
to prevent accidents, zoonotic infections, etc.? 


yes 


QO no 


Rodent control program used: 
Le glraps 
Ze seLOLSOuS 


3. No program 
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Form Code 


CODE/COL. 


No 


X212 


X213 


XZID 


X216 


X217 
X218 


code 
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Ss 


14. 


If a rodent control program exists, how would the farmer 


rate its success? 


1. successful 


2. partially successful 


32. snot successful. __ 


a ree ct a en rm 


Other management indicators which may be pertinent _ 
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Form Code 


_CODE/COL. 


Hi) eee 
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Pore? APPENDIX 2 


Poultry Disease Research 


School of Veterinary Medicine 


Tuskegee Institute 


Flock Testing Feasibility Study 
CODE/COL 










Ii. Cross-Sectional Sampling ~ Processing Plant Data 


NOTE: Date should be collected from birds for which farm 
data has already been obtained 


Date (mo/dy/yr) Visit No. (1,2,3,4) 


Data collected by (code)_ 


Processing Plant Code Location (city/county 





Firm Code __ Liesl 42 Bae 
Source Farm code 
Distance from Tuskegee (miles) _ Se) ek 
A. Pre-Processing Data: 
1. Number of birds shipped from farm for processing X221 
2. Average slaughter weight on date of marketing (lbs) 
X222 
3. Distance from farm to plant (miles) K223 
4. Number of hours of driving from farm to processing 
plant _ Re X224 
5. Trucks and crews owned by plant? K225 a 
Le yecr—is= oer | 
2 no ews. 
6. Number of different drivers involved in flock 
delivery. Shs fe, = eee Sata | X226 
7. Number of trips to haul the whole flock to the plant 
eee X227 
8. Number of trucks used to transport birds | X228 
9. Type of coops used to transport birds _ X229 
1 wooden frame X266 
2 plastic A275 


3 steel X230 
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APPENDIX 2 







CODE/COL 






10. Approximate number of birds per coop (cage) X230 
ll. Approximate number of birds per truck 
12. Estimate time of day when birds from the farm were 
processed and inspected 
13. Method of cooling birds in transit: 
i ffansy 
O others.en 


14. Method of cooling birds at processing plant: _ 


1 fans ex 

0 others 
wseenumper Of birds dead On arrival-at plant... 4/ X235 
16. Total number of birds slaughtered _ X236 





17. Number of days involved in marketing this flock (from 


Piectelomiaseebirdeto leave the farm) ) X237 


B. Macroclimate data on date of processing: 





1. Average temperature (F) for the day X238 

2. Was it rainy or dry (rainy=l, dry=0 X239 

3. Average relative humidity (%) for the day X240 

C. Inspection 

1. Number of inspectors on line em X241 

2-2 Number of veterinarians(s) in plant. §. —s_—s—_ __ | X242 

3. Line speed (average per day) per min. ss X243 : 
4. Type of inspection: 


Pmerands: OL ia trte. foie 
2. semi-traditional 

Bee Craditional 

Leathers 


D. Processing Plant: 


be processing plant No. 





processing plant No. 





processing plant No. 


processing plant No. 





processing plant No. 
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CODE/COL 









Daily total capacity of processing plant ‘ X251 he. 


Number processed this day 
Total number of employees in plant _ er eee 


Number of hours per day that plant operates _ 


Number of shifts X255_ 


E. Condemnation data: 


1A. 
Farm ID 


Study 


farm 


Totals 
Ali Farms A 


got te) eh lee 


Totals 
All Earns 4 


B 
C 
D 
E 


Grand 
Moma 


Condemnation due to diseases for specific farms on this day: 


No. inspected No. condemned Cause of Con- Lbs. trimmed 
and passed (whole carcass) demnations due to diseases 
7 kok & ? YS 
Kk _. Dxee 
~ Ph Oe See ze KK bee a Y2u 
kK * ee 
kk x ; ¥2a 


ween te URS ea 2 ee ee ee | cove/co. 


Average condemnation rate due to diseases for process- 








ing plant (this day) (whole carcass) Yi 








Condemnation rate due to diseases for study farm 


(whole carcass) nme. x2 


Disease specific condemnation rates for all farms 
(whole carcass) 












** denotes values to be computed 


a ne 


Disease specific condemnation rates for study farm 
(whole carcass) 








Overall condemnation rate due to diseases for study 
farm (whole carcass + trimmed parts _ ¥15 


ae an nn a a her 








i) 


Te 


Farm ID 


Study 


Farm 


Totals 
All Farm 1A 
B 


ri 


Totals 


Grand Totals 


No. condemned 
(whole carcass) 





a 


a a ee 


a ee 


a re 


rt ee 


Je 
ray 


Non-disease condemnations: 


Lbs. of parts 


condemned (trimmed) 


(ibs) 





a Rn a 


a ee a 


en ce 


Hee oa EE a ie 


Denotes values to be 


Cause of Condem- 
tations 


me 
De re a ne re 
a 
re vane 


2 Soe 


me ee ee 


ce 


ooo 
re on ne 


computed 


ee 
be 
Be 


How often are condemned birds disposed? _ 
Method of disposal OF Wi © 
Do birds fall onto the floor? 

1 yes often 

2 occassionally __ 


3 none 


**k denotes values to be computed 


2% 





Average condemnation rate due to non-diseases for processing 


processing plant (this day) (whole carcass) __ 





Condemnation rate due to non-diseases for study farm 
(whole carcass) 


os 





Overall condemnation rate (whole carcass + parts trimmed) 
due to non-disease cause 


a a 





1 


Miscellaneous data: 


1. Unloading and slaughter: 
a. Number of mean used to take birds from coops ____ 
b. Number of stickers (killers) 
c. Birds per minute per man 
d. Number of back up killers employed 


a 


e. Average length of bleeding time (min, ) Sauk 
2. Scalding: 
a. Type of scalder . 
b. Proportion of birds still struggling when entering 
scalding vat 
c. Temperature of scald water(F) 


d. Estimate of time in scalder (min.) 


e. Kind of additives used in scald water (if any) _ 


eee nt 


3. Other information: 
a. Number of break downs in line during average working 


day 





b. Distance between shackles (in.) 
c. Are there separate cans for different types of 
condemned birds? 
1 yes 
0 no 
d. Are contaminated birds specially marked? 
1 yes 


O no 


Y16 


gly 


Y18 
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[cone /cot 


Zao 


Appendix 3 


Feed Utilization Procedure 


-Volume of truncated cone was determined by taking 
measurements of the larger and smaller diameters and the 
side(s) of the cone, which were then used in the general 
equation for a truncated cone: 


Ve aus + RORD + R.) 


h 
= eee ee 
3 


Where? RS is the larger diameter, 


Ro is the smaller diameter, and 


h is the vertical distance of the cone (See figure on 
the next nage 


The volume of each right circular section on top of the 
cone was determined by measuring the distance between 
that distance by the number of rings. To aerive the 
vertical distance of the section, the diameter of the 
section was taken from the larger diameter of the 
truncated cone. These measurements were then used in 
the general equation for a right circular cone, 


Water a) 
Te 


Wheye: D is the diameter and 

Metis the height of the section. 

The volume of the section was multiplied by the number 
of sections. The total volume of the bin was taken to 
be the sum of the volume of the cone and the sections. 


-To determine the weight of one cubic foot of feed, the 
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Corrugated 
Rings 


TYPICAL FEED BIN 





total weight of the feed in the bin after it was filled 


(obtained from feed receipts), was divided by tne total 
voiume of the bin. 


-To determine the amount of feed the birds consumed per 
week, the feed level in the bins was measured on a 
weekly basis. The level of feed was determined by 
hitting the side of the bins until the feed was leveled 
off. Based on the feed level, the volume of the ieft 
over feed was determined in the same manner as that of 
the method of eorec rae ae the volume of the hins, Tne 
weight of feed left in the bin was obtained by 
multiplying the volume of the ieft over feed by the 
weight of per cubic foot of feed. Therefore, the amount 
of feed consumed by a flock could be determined by 
subtracting tne weight of feed left in the bin ata 
Specific time from the total weight of feed placed in 
the bin. From the above procedure, the feed consumption 
of individual birds could also be calculated by dividing 
total feed utilized by the number of birds in the flock. 
This method provided the average amount of feed consumed 


Dermoid, oer bin, per fiil up. 
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Appendix 4 


This appendix refers to all the computer printouts turned in to 
FSIS at the termination of the project. 


Appendix 5 
Report of Bacteriologic Identification 
Poultry Disease Research 
® School of Veterinary Medicine 
Tuskegee Institute 
Source of Sampie: 


Hatchery ( ) Poultry House ( ) 


Name of Company 


Name of Hatchery ne ie Soe 

Name of Grower Mina ne 

Date of Original Sample Taken + De BAM oR 
Date of Identification Attemptea Boo) Nie Si. alee 
I. Bacteria 


Gram stain: Positive ( ) Negative ( ) 








) 
) 


) 


Morphology: Rod is) Cocel eon a, 
a. Blood agar 
Colony morphology: Smooth a) Rough ( ) 
| Size of colonies: Large ( Medium ( ) Smad let 
Hemolysis: Complete ( ) Viridan ( ) No ( 
! b. MaConkey agar 
Colony morphology: Smooth Cae) Rough (9) 
| Size of colonies: Large CR) Medium) Gi a Sma lige ( 
Color: Red eG oh deat Ge ecoloniless (=) 
e. Entero Tube II. Y GS Gane) No ( ) 
ieeeGlUp eee) 2 .GAS (i) .) Age) MNES) CoS) 
ee -ORN: (ons) His Hs Cie) IND (f3) 
Bee S DON Cees) Cae LAG a 
Tee A Gene) SSO n Pie) 
" Oe VE Ga) Oe DU Ie (eae) PA ao) 
| 11, Urea ( ) ioeeciTs 
d. Others 
Diagnosis: 
Pile Fungus 
a. Based on microscopic examination of morphology ODLCURUITEe. 
Description: 
Beescolony color: 


Diagnosis: 


i ici Date 
Diagnostician az tisge, deere coat lo 


Te EN 
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FLOW CHART 


Appendix 5a 


Gram Positive Cocc@ 
No growth on MacConkey 





Inoculate Trypticase broth and TSI if necessary 


Gata lage eee ee eee zone fp NE alpha or non-hemolytic 


complete hemolysis Streptococcus sp 


Pos 
Pos Beta hemolytic Streptococcus sp 


Neg 
Double-zone hemolysis ——-——-—_——-TSI: Acid butt 


Neg 
or zone of partial hemolysis : 3 


Micrococcus 


Coagulase Neg . Staphylococcus 
epidermidis 
Pos 


Staphylococcus aureus 5S. aureus 


a) A broth culture can be gram stained to confirm morphology including 
chaining of Strptococcus. 


b). TSI: The acid butt may be just turning yellow at 24 hours, 


*Source: John N. Berg. Diagnostic Veterinary Bacteriology Procedures and 
Species Identification. AAHA Annual Meeting Proceedings, 1979 (54). 


SO9TUOTOD T[BSUNT oj Jejei seoeds 1eay [Tee uz09 vu UE sisquny 
eds AsyUODIY pue Aese pooTq ey} UT sAequnyN » 
t 
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Appendix 5c 
Identification of Bacteria and Fungi Isolated 
From Hatcheries and Broiler Houses 
CODE ORGANISM FREQUENCY. 
A-4 Micrococci sp 2 
Bacillus sp ah 
Staphylococci sp Ys 
Cryptococci sp Z 
Moraxella sp. or Pasteurella sp. 
or Pseudomonas sp. or Group 2K-1 
Pseudomonas like or Alcaligene sp L 
Enterobacter agglomerans or 
Shigella boydii 
Enterobacter agglomerans or Citro- 
bacter sp ik 
BP-B55| Micrococci sp ae A 


Proteus jnealsesioyal J ahs 


Ep-B45| Pseudomonas sp. or Alcaligenes sp 
or Group M-4 Moraxella like Z 
Salmonella enteritidis bioses 
para A or Citrobacter sp.. 
Pasteurella hemoltica 
Listeria sp. 
2K-1Pseudomonas like 
Bacillus sp. 
Enterobacter agglomerans 
Non hemolytic streptococci 
Staphylococci sp 
Corynebacterium pyogenes 
Microcoeci sp; 4, 
Planococcus sp. L 
Aerococci sp. 2 
Group A Pseudomonas like or 
Pseudomonas aeruginosa L 
Escherichia coli or Klebsiella | 
" ozanae i | 3 
Pseudomonas stutzeri : 





BPN WeERP RPP 
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Appendix 5c Identification of Bacterial and Fungi Isolated 
From Hatcheries and Broiler Houses (cont.) 


CODE yx: |p oaen ene ORGANESHE 7 FREQUENCY 





Enterobacter sp. ox Citrobacter sp 1 
Enterobacter sp i 
Saccharomyces sp. 1 
Penicillium sp. a 


|Micrococci sp 

| Staphylococci sp. 
Bacillus sp. 
Aspergillus sp. 
Aspergillus sp. 
Penicillum sp. 


Meee ww 


NO 


Rhyzopus sp. 
Aspergillus sp. 
Penicillum sp. 


aed 


Pseudomonas cepacia 
Listeria sp. 
Corynebacteria sp. 
Pseudomonas sp. 
Micrococci sp. 
Staphylococci sp. 
Penicillum sp. 
Penicillium sp. 


oe ee) 


Aspergillus sp. 1 
Trichosporon sp. " 
Mucor sp. i 


Aspergillus sp. 2 
Penicillium sp. M2 


Penicillium sp. en 


Corynebacteria sp. 

Listeria sp. 

Bacillus sp 

Aerococci sp. 
Staphylococel isp) =. 


ln wo Ne 
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Appendix 5c Identification of Bacterial and Fungi Isolated 
From Hatcheries and Broiler Houses (cont.) 




















































_____ ORGANISM 


Aerococci sp. 

Citrobacter sp. 

Alcaligenes or Group 2K~1 
Pseudomonas like 

Citrobacter sp. or Enterobacter 
agglomerans 

Pseudomonas sp. 

‘Enterobacter agglomerans 
Shigella flexneri 
Non-hemolytic streptococci 

Pseudomonas sp. or Alcalicines sp. 
or Group 2K-1 Pseudomonas like or 
Group M-4 Moraxella like 





Staphylococci sp. 

Micrococe1 ssp. 

Pseudomonas sp. 

pseudomonas sp. 

Streptococci sp. (alpha hemolytic) 
Streptococci sp. (Non hemolytic) 

Aerococci sp. 

Escherichia coli 

Moraxella-like Group M-4 

Enterobacter aerogenes or 
Enterobacter cloacae 

\Enterobacter sp. 

Aspergillus sp. 

1Cryptococcus Sp. 

Kloeckera apiculata 

Aspergillus sp. 





OC eo 


ee 






Appendix 5c 


| Group 5A Pseudomonas like or 


}DL ~ D33/ 
D53/D63 | 


lp = $D147 
D24/D44 


Identification of Bacterial and Fungi Isolated 
From Hatcheries and Broiler Houses (cont.) 


ORGANISM 


Pseudomonas aeruginosa 


Aspergillus sp. 


B-hemolytic streptococci 
Pennicillium 


Aspergillus sp. 


|DL ~ D54 


E-1 


|Bacillus sp. 


Corynebacteria sp. 


|Aspergillus sp. 


Pullularia sp. 


Pseudomonas cepacia 
Listeria sp. 
Aspergillus sp. 


Pseudomonas cepacia 
Listeria sp. 
Aspergillus sp. 


Escherichia coli 
Pseudomonas sp. 
Bordetella bionguiseptica 
Pseudomonas cepacia 
i\Pseudomonas aeruginosa 
|Alcaligenes dinitrificans 
Listeria sp. — 
Cornebacteria sp. 
Klebsiella sp 

Providencia sp. 
Enterobacter cloacae 
Acinetobacter luoffi 
IStaphylococcus sp. 
Aerococci sp. 


Micrococci Sp. 





a Ce Or 


| 
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Appendix 5c Identification of Bacterial and Fungi Isolated 
From Hatcheries and Broiler Houses (cont.) 


CODES | _____ ORGANISM ” FREQUENCY 


DL — F45 | Micrococci sp 4, 
Shigella sp. 2 

| Klebsiella pneumoniae I 

| Alternaria sp. L 

| 


Penicillium sp. 


Dinar o> | Bacillus sp. 3) 
| Staphylococci sp. 5 
Micrococci sp. 6 
Alcaligenes sp. or Achromobacter 
icp. 
| Listeria sp. 

Proteus sp. 

}Non-hemolytic Streptococci sp. 
Corynebacteria sp. 

Enterobacter agglomerans 

| Alcaligenes dinitrificans 
Enterobacter cloacae 

Pseudomonas sp. or Achromobacter 
Sp. 

Aerococci sp. 

Klebsiella pneumoniae 

Acinetobacter luoffi 

Pseudomonas stutzeri 

Enterobacter sp. or Serratia sp. 





PRR RRR 


F-J] |‘Trichosporon sp. 
Micrococci sp. 
|\Staphylococci sp. 
Corynebacteria Sp. 
'Bacillus sp. 
Escherichia coli 
Serratia liquefaciens 


= BS te 1 





F-2 Trichosporon sp. a 
Bacillus sp. 2 
Staphylococci sp. 4 
\Micrococci sp. 4, 


fo a re 





nN 
sw! 
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Appendix 6 


BROILER SERA TITERS FOR NEWCASTLE DISEASE VIRUS AT FOUR WEEKS OF AGE 


Flock Testing Feasibility Study, Tuskegee Institute (1981/82) 
. 
CODE | HOUSE TOTAL 


H-1 
ae H-2 


1:80 | CONDEMNATION RATE 









16.67 
330] 933233 
16.67 
66.67 | 16.67 
50.0 33.33 


| 
| 





ApS 2ela 1 $e eee cen CO prs oot). 
lA-62 Hey 4 e100 
ae eye tae 100 cer Mitte UU 
| 
Verbal) ies | 4 | 100 
a oareo 4. ea 0 Eten Sp 383 (el ee 
one Hal ME 2) S100 Fee! 
| : | 
i | | | 
aos | wa | Ne Seek ae | 
| | \ 
FEES | Het | 4 100 
| be SH? a 100 my 
tee | 
(AZAOuRy Hod | 4 100 

| H-2 4 | 100 

Ho 3) | een 00 

eee he ae 


| 
| 


i ar 


‘copE | HOUSE | TOTAL 
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Appendix 6 (cont.) 


BROILER SERA TITERS FOR NEWCASTLE DISEASE VIRUS AT FOUR WEEKS OF AGE 


Flock Testing Feasibility Study, Tuskegee Institute (1981/82) 

















reo 


CONDEMNATION RATE 






















B-11 H-1 5 20.0 80.0 
6 5010 


B21 | H-l 16 62.50 

| Melee ee LG i+ iW ie) 

B-31 | H-1 | 10 | 40.0 NOC NO 

\; | 

B-51 | H-1 | a fiilely) ipa Damon 9 eo ta 
ecm lin | Wo 27327 18¢18. |:45.4 


CODE 


B=52 


260 
Appendix € (cont.) 


BROILER SERA TITERS FOR NEWCASTLE DISEASE VIRUS AT FOUR WEEKS OF AGE 


Flock Testing Feasibility Study, Tuskegee Institute (1981/82) 











1:80 | CONDEMNATION RATE 





OFZo 


0.0384 








0,231 


ee nes eee ef 


0.344 
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| Appendix 6 (cont. ) 

| BROILER SERA TITERS FOR NEWCASTLE DISEASE VIRUS AT FOUR WEEKS OF AGE 
| 


| Flock Testing Feasibility Study, Tuskegee Institute (1981/82) 


|CODE 


| 







b80 CONDEMNATION RATE 









a8) 
O 
Gq 
17) 
tr 
tH 
[oO 
ry 
> 
te 
7 
La 
fe 
= 
io) 


‘\B-65 | H-1 


\D-121 


J ieee een 
uw N 
ee i 


T 
i= 
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pee 
N 








jas WEI 


nieteasiene 
n 
Ww 
N 
T 
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a ee eee EIS ows 
i 
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N 
a 
ft 





0.477 
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Appendix 6 (cont.) 


BROILER SERA TITERS FOR NEWCASTLE DISEASE VIRUS AT FOUR WEEKS OF AGE 


Flock Testing Feasibility Study, Tuskegee Institute (1981/82) 











CODE | HOUSE | TOTAL  boeaat 1:10 | 1:20 | 1:40 | 1:80 | CONDEMNATION RATE 
Date) MeneT 4 
H-2 4 
H-3 4 _... |b 3 800R ae 
902 
023 
002 
865 
982 
960 
fuer Os 6/4 
0.579 
Does ek 4 100 
H-2 4 100 
H-3 4 75.0 25.0 
H-4 4 we TAOA® 2539 0.449 
eal pmeecenper base aa | 
| H-2 4 100 0.868 
D-55 | Hel 4 | 100 | 
| H-2 i544 R100 | 0.478 







Appendix 6 (cont. ) 


BROILER SERA TITERS FOR NEWCASTLE DISEASE VIRUS AT FOUR WEEKS OF ACE 


Flock Testing Feasibility Study, Tuskegee Institute (1981/82) 


NAY 8A Pe A ee oot 


HOUSE | TOTAL NT_ 1:80 | CONDEMNATION RATE 


50.0 
po 08905 Dene 





20.4940 





ZoeU 0.380 





er shake 





ie es a i Sl 
Teas 


0.296 


0.754 











264 
Appendix 6 (cont.) 


BROILER SERA TITERS FOR NEWCASTLE DISEASE VIRUS AT FOUK WEEKS OF ACE 


Flock Testing Feasibility Study, Tuskegee Institute (1981/82) 





eee ee rr a RR NN A oR 


HOUSE | TOTAL eeu CONDEMNATION RATE 





1:40 
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INITIAL VARIABLES USED FOR DEVELOPING A PREDICTIVE MODEL FOR 
POULTRY CARCASS CONDEMNATION RATES IN PROCESSING PLANTS. 


VARIABLE DEFINITION 


Firm/Farm/Brood Code 

Hatchery Source #1 

Hatchery Source #2 

Hatchery Source #5 

Hatchery Source #4 

Hatchery Source #5 

Hatchery Source #6 

Hatchery Soruce #7 

Incubator Type #1 

Incubator Type #2 

Average Fertility Rate (%) For Hatchery 

Average Hatchability Rate (%) For Hatchery 

Average Fertility Rate (%) For Brood of Birds 

Average Hatchability Rate (%) For Brood of Birds 

Number Of Days Eggs Held On Farm 

Number Of Days Eggs Held In Hatchery 

Number Of Breeder Flocks Supplying Eggs for Hatchery 
Location Of Air Intakes In Hatchery (Roof) 

Location Of Air Intakes In Hatchery: East, North East 
Location Of Air Intakes In Hatchery: North, North West 
Location Of Air Intakes In Hatchery: West, South West 
Location Of Air Exhaust From Hatchery (Roof) 

Location Of Air Exhaust From Hatchery: East, North East 
Location Of Air Exhaust From Hatchery: North, North West 
Location Of Air Exhaust From Hatchery: West, South West 


Date Of Setting Of Eggs For Brood (Calendar Day Starting at Jan. 


Total Number Of Birds In (Farm/House) 

Number of Birds Shipped From Farm For Processing 

Average Slaughter Weight On Date Of Marketing (lbs.) 
Proportion (Z) Of Undersized Birds 

Strain Group l 

Strain Group 2 

Strain Group 3 

Total Number Of Strains On Farm 

Total Number Of Breeder Flock Sources For Farm 

House Type l 

House Type 2 

Average Population Density Per Square Feet 

Average Distance Between Houses On Farm 

Ventilator Type 1 

Number Of Ventilators 

Location Of Air Intakes In Broiler House: East, North East 
Location Of Air Intakes In Broiler House: North, North West 
Location Of Air Intakes In Broiler House: West, South West 
Location of Air Intakes In Broiler House: South, South East 


1) 


a, 


APPENDIX 9 (cont.) 


Location 
Location 
Location 
Location 


Age Specific Mortality Rate 
Age Specific Mortality Rate 


Of Air 
OiPALE 
Of Air 
of Air 


VARTABLE DEFINITION 


Exhausts 
Exhausts 
Exhausts 
Exhausts 


In 
In 
In 
in 


Broiler 
Broiler 
Broiler 
Broiler 
(%) At 4 Weeks 
(%) At 7 Weeks 


House; 
House; 
House: 
House; 
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East, North East 
North, North West 
West, South West 
South, South East 


Date Of Placing Of Chicks For Brood Starting Jan. l 
Length Of Time For Preheating (Hours) 


Brooding 


Type 


Length of Stay Of Birds In Broiler House (Day) 


Age Broiler Guards Removed (Day) 
Number Of Chicks Per Broiler Hours (Average) 


Average 
Disease 
Disease 
Disease 
Disease 
Disease 
Disease 


Disease Rate (%) Based On Necropay At 7 Weeks 


Specific 
Specific 
Specific 
Specific 
Specific 
specific 


Rate 
Rate 
Rate 
Rate 
Rate 
Rate 


(%) 
(Z) 
(%) 
(%) 
(4) 
(7%) 


Mash (Type Of Feed) 
Crumbles (Type Of Feed) 
Feeder Space Per Bird 


Feed 
Feed 
Feed 
Feed 
Feed 


Sou 
Sou 
Sou 
Sou 


Number of Feed Deliveries Per Month 


ree #1 
ree #2 
rce #3 
ree #4 


Spillage (Present Or 


Based 
Based 
Based 
Based 
Based 
Based 


Not) 


on 
On 
On 
On 
On 
On 


Necropsy 
Necropsy 
Necropsy 
Necropsy 
Necropsy 
Necropsy 


Number of Waterers Used During Brooding 
Water Spacing Per Bird 
Age Manual Waterers Removed (Days) 
Water Spillage (Present Or Not) 
Average Feed Utilization Rate For (Farm/House) At 4 Weeks 
Average Feed Utilization Rate At 7 Weeks 

Duration Of Prophylactic Medications (Days) 
Growth Promoters Only 


Vaccine Manufacturers 
Vaccine Manufacturers 
Vaccine Manufacturers 


#1 
#2 
#3 


Farm Sanitation Index 


Litter T 
LLecer . 


ype l 
ype 2 


Depth Of Litter (Inches) 
Nubmer of Batches In Litter 
Average Maximum Temperature For Study Blood (F) 
Average Temperature For Study Brood (F) 


(1 = Leucosis) 

(2 = Septicemia-Toxemia) 
(3°-= Airsacculitiep 

(4 = Synovitis) 

(5 = Tumors) 

(6 = Others) 
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Appendix 9 (cont.) 


VARIABLE DEFINITION 


“3 


Average Minimum Temperature For Study Brood (F) 
Average Total eee For Study Brood (inches) 
Number Of Rainy Days 

Fall (Season Brood Of Birds Raised) 

Winter 

Spring 

Distance From Farm To Plant (Miles) 

Number Of Hours of Driving From Farm To Plant 

Number Of Different Drivers Involved Ia Flock Delivery 
Number Of Trips To Haul Whole Flock To Plant 

Number Of Trucks Used To Transport Birds 

Number Of Birds Per Coop 

Number Of Birds Per Truck 

Number Of Birds Dead On Arrival At Plant 

Number Of Days Involved In Marketing Flock 

Number Of Inspectors On Line 

Semi--Traditional Inspection 

Processing Plant #1 

Processing Plant #2 

Processing Plant #3 

Processing Plant #4 

Daily Total Capacity Of Processing Plan‘ 

Number Processed For Entire Day 

Total Number Of Employees In Plant 

Number Of Hours Per Day Plant Operates 

Number Of Shifts 

Average Condemnation Rate (%) Due To Diseases For Processing Plant 
Condemnation Rate (%) Due To Diseases For Study Farm 


Disease Specific Condemnation Rate (%) For All Farms (1 = Leucosis) 

Disease specific Condemnation Rate (%) For All Farms (2 = Septicemia- 
Toxemia) 

Disease Specific Condemnation Rate (4%) For All Farms (3 = Airsaceculitis) 

Disease Specific Rate (%} For All Farms ‘4 = Synovitis) 

Disease Specific Rate (%} For Al] Farms (2 = Tumors) 

Disease Specific Rate (%) For All Farms (6 = Others) 

Disease Specific Condemnation Rate (4%) for Study Farm (1 = Leucosis) 

Disease Specific Condemnation Rate (Z) Yor Study Farm (2 » Septicemia- 


Toxemia ) 
Disease Specific Condemnation Rate (%) For Study Farm (3 = 
Aupsacculitis) 
Disease Specific Condemnation Rate (4%) For Study Farm (4 = Synovitis) 
Disease Specific Condemnation Rate (%) For Study Farm (5 = Tumors) 
Disease Specific Condemnation Rate (%) For Study Farm (6 = Others) 
Overall Condemnation Rate (4%) Due To Diseases For Study Farm 
Average Condemnation Rate (%) Due To Non-Diseases For Processing 
Plant 
Condemnation Rate (%) Due To Non-Diseases For Study Farm 
Previous Condemnation Rate (%) Due To Diseases For Study Farm 


Appendix © {cont.) 


VARIABLE, CODE VARIABLE DEFINITION 

Nov isease Specitic Condemnation Rate (4) ot previeus 
For Study Farm (1 = Leucosis) 

YLT Disease Specific Condemnaiton Rate (%) of previous 
For Study Farm (2 = Septicemia-Toxemia) 

deta Disease Specific Condemnation Rate (4) ot previous 
For study Farm (3 = Airsaculitis) 

VEZ Disease Specific Condemnation Rate (%) of previous | 
For study Farm (4 = Synovitis) 

YI32 Disease Specific Condemnation Rate (%) of previous 
For Study Farm (5 = Tumors) 

C147 Disease Specific Condemnation RAte (%) of pervious 


For Study Farm (6 = Others) 
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